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Almost Fifty Years of URSI’ 


Address at banquet, December 13, 1960, of URSI-IRE Meeting, Boulder, Colo. 


J. Howard Dellinger 


(January 19, 1961) 


“URSI’—If I were Professor Kennelly, I would 
automatically say “URSI, the bears.’ Although 
a professor of electrical engineering, being a Harvard 
professor, Kennelly could not forget his Latin so to 
him URSI was always the bears. That remarkable 
man was one of the early presidents of URSI; was 
the discoverer of the ionosphere (some of you have 
heard of the ‘“‘Kennelly-Heaviside layer’), and was 
the father of the coded Ursigrams. 

All of us are in on a good thing. We are in on the 
ground floor because “URSI’s almost 50 years”’ in- 
clude three-fourths of the entire history of radio and 
practically all of the history of radio science. This 
is true even though the greatest single event in the 
history of radio science was its beginning, just a 
hundred years ago. That was Maxwell’s discovery 
by means of pure theoretical science, that electro- 
magnetic waves could be produced by electrical 
apparatus and that their velocity would be the 
velocity of light. Maxwell published this in 1864. 
Once this giant stride had been taken, however, 
science had little to do with radio for the next 50 
years. The world indeed took Maxwell’s hint and 
went about the business of producing and using those 
electromagnetic waves from electrical apparatus, in 
other words,radio. Therewasofcoursesome research 
on the means of generating and detecting the waves, 
but little in the way of scientific study either of 
the devices or the behavior of the radio waves or 
the medium through which they are propagated. 

An awakening consciousness that science and radio 
should get together was what led to the formation of 
URSI. The first step was taken in 1913, about the 
same time, incidentally, as the beginning of the 
Institute of Radio Engineers. In that year a group 
of nine scientists from six European countries met 
in Brussels and determined to establish an “Inter- 
national Commission on Scientific Wireless Teleg- 
raphy,” to undertake research on wave propagation. 
World War I prevented the organization from func- 
tioning. After the war, in 1919, the members of the 
group saw that the need for organizing radio 
science had greatly increased and they did establish 
the “International Union of Scientific Radio Teleg- 
raphy” with approximately the present objectives 
and form of the URSI. The first General Assembly 
was held in Brussels in 1922. We radio scientists 


1 Editor’s note: Dr. Dellinger is an honorary president of the International 
Scientific Radio Union (URSI), was the first chief of the Central Radio Propa- 
gation Laboratory, and was president of the Institute of Radio Engineers in 1925. 





have been having fun ever since, talking to and 
learning from one another. 

The 1963 General Assembly in Japan will celebrate 
the 50th anniversary of the 1913 beginning. It 
could well recognize at the same time the centenary 
of radio’s origin, Maxwell’s theory of electromagnet- 
ism. That supreme achievement is to be celebrated 
also at the 1965 Plenipotentiary Conference of the 
International Telecommunication Union. 

Radio today talks an entirely different language 
from that of 1913. You note that the organization 
formed that year dealt with ‘wireless,’ not “radio.” 
Also in those days we always described the frequen- 
cies we used by means of wavelength in meters; 
nobody dreamed of a megacycle, and what we know 
as a frequency meter was called a wavemeter. Also, 
the “ionosphere” had not been heard of, what we 
knew was the “Kennelly-Heaviside layer.’ To 
intrude a personal note, I had considerable to do with 
the change to all three modern expressions: radio, 
kilocycle, ionosphere. This I was able to do by 
virtue of my position as head of the radio work of 
the National Bureau of Standards. In the U.S. the 
change in common practice from ‘wireless’ to 
“radio”? came in about 1912 to 1916, from wave- 
length to frequency about 1916 to 1920, and from 
‘“Kennelly-Heaviside layer’ to “ionosphere” in 
about 1932. In Bureau papers published in the 
early thirties, it became more and more unsatisfactory 
to talk about the ‘“Kennelly-Heaviside layer’ 
when we regularly got reflections from two or more 
layers. These changes in terminology came some- 
what later in other countries. It is interesting to 
note that the changes from ‘‘wireless’”’ to “radio,” 
and from wavelength to frequency, were resisted by 
the British and are to this day. The British just 
don’t like change. 

Not only the language but the ideas were different 
in the days of URSI’s beginnings. Radio telegraphy 
was carried on by means of damped waves, and radio 
telephony which requires continuous waves was just 
beginning to be experimented upon. The damped 
waves were in use because only they could be gener- 
ated with enough power. Damped waves are a 
terrible nuisance, not only in practice because of 
their great interference production (like pulses) but 
also in all calculations and theory. The principal 
thing you had to deal with was a quantity called the 
logarithmic decrement, a measure of the rate at 
which each train of waves or current oscillations fell 
off to zero. Well, the decrement of a circuit was 
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proportional to the ratio of resistance to reactance 
or the reciprocal of Q, the sharpness of resonance; 
the trouble was you had to take account of this 
quantity for the damped current oscillations as well 
as for the circuits. The principal instrument a 
radio inspector had to carry around was a decremeter, 
a device for measuring decrement as well as wave- 
length. 

When URSI was beginning we were just develop- 
ing reliable sources of continuous waves. Not only 
did this make radio telephony possible but there was 
dawning hope of that bugbear, decrement, disappear- 
ing from radio calculations. Anticipating the demise 
of damped waves, the National Bureau of Standards 
during World War I prepared the first radio text- 
book which based radio theory on straight alter- 
nating-current theory, giving damped waves only 
minor and separate treatment. This was Bureau 
Circular 74, Radio Instruments and Measure- 
ments.” Issued 42 years ago, it is still being sold by 
the Government Printing Office. I suppose it is the 
oldest of the Government’s best sellers. 

We tried various ways of producing continuous 
waves (Alexanderson’s high-frequency alternator, 
arc generator, and thermionic tube). J. A. Fleming 
of England had invented the two-electrode thermi- 
onic tube, and used it as a radio detector, in 1899. 
In 1906 De Forest put in the third electrode and the 
development of the electron tube as oscillator and 
amplifier proceeded steadily from then on. The 
way was clear for radio telephony, for precise radio 
measurements, for penetrating the mysteries of 
radio wave propagation, and for the new world of 
electronics. The advances have snowballed ever 
since. 

Two of the thirteen General Assemblies of the 
URSI were held in the U.S. The second was in 
Washington in 1927, and the twelfth here in Boulder 
in 1957. I think E. V. Appleton and I are the only 
participants in the 1927 meeting who are still active 
in URSI. Another participant was the late Samuel 
S. Kirby (father of Bob Kirby who is present this 
evening); he was a very active URSI man and was 
secretary of the USA National Committee, which 
we then called the American Section. Also among 
those present was the late Dr. B. van der Pol, 
eminent mathematician who at the 1927 meeting 
initiated the Commission on Radio Waves and 
Circuits (called at first Radio Physics) and was its 
chairman for 25 years. Distinguished participants 
who are still living included Dr. A. H. Taylor, 
E. F. W. Alexanderson, H. Pratt, Dr. P. Le Corbeil- 
ler, Dr. G. C. Southworth, Dr. G. Breit, Dr. M. 
Tuve, Dr. C. B. Jolliffe. Also present were Dr. 
L. W. Austin and Professor A. E. Kennelly, both of 
whom served terms as chairman of the USA National 
Committee and president of the URSI itself. The 
president in 1927 was the French General G. Ferrié 
and the secretary general was professor R. B. 
Goldschmidt of Belgium. Professor Goldschmidt 
was the father of URSI. When the predecessor 
organization was founded in 1913 he offered to build 
a special radio transmitting station for research 








purposes, and he financed the URSI in its early 
vears. General Ferrié was a remarkable leader. 
He was the first president and provided the guide 
lines of growth for the first 10 years. He was the 
finest type of gentleman, a lovable character and an 
inspiring scientist. Dr. Austin was the second 
president of URSI. He was the first chairman of 
the USA National Committee and was the American 
pioneer in radio research. He was doing this work, 
on the staff of the U.S. Navy at the National Bureau 
of Standards, as early as 1905, and was active in it 
until his death in 1932. He was most famous for 
his semiempirical formula for long-distance received 
field strength for low frequencies (LF and VLF). 

Don’t start to worry. I am not going to recite 
the history of the URSI. That is largely available 
in the URSI Information Bulletins, and a compre- 
hensive publication on it will be issued as part of 
the golden anniversary celebration of 1963. 

The list of places where the General Assemblies 
have been held will convince you that we who are 
designated to attend these affairs are fortunate. 
Besides Brussels and Washington they have included 
Copenhagen, London, Venice, Paris, Stockholm, 
Zurich, Sydney (Australia), The Hague, and Boulder. 
Besides giving the participants a good time, they 
have had a mighty influence in coordinating and ex- 
pediting the extraordinary advances of radio; and, 
just by the way, also contributing to that good inter- 
national understanding which is the prerequisite of 
world peace. 

While recounting the privileges that URSI has 
brought to me, I might add that I had the honor of 
serving in the USA National Committee as Secretary 
from 1920 to 1933, Vice Chairman from 1933 to 1940, 
and as Chairman from 1940 to 1949; and in the 
international URSI as Chairman of the Commission 
on Radio Wave Propagation from 1934 to 1946, as 
Chairman of the Commission on Radio Measurement 
Methods and Standards from 1946 to 1952, as Vice 
President from 1934 to 1952, and as Honorary 
President since 1952. 

Since my service as URSI Vice President spanned 
the period of World War II, I had the opportunity to 
help in keeping the organization alive during that 
distressful period. For four years after Germany 
invaded Belgium in May, 1940, mail service between 
England and Belgium was cut off, so President 
Appleton in London and the Secretariat in Brussels 
had no communication. URSI work in the Na- 
tional Committees in the various countries also 
declined to nearly zero. I wrote to the Brussels 
Secretariat in August 1940. I received a reply two 
months later from Miss Straetmans, the office 
assistant, saying that she was trying to keep the 
office going but that Commander Dorsimont (Acting 
Secretary General) and Major Herbays (Editor) 
were prisoners of war in Germany. She gave me 
their addresses; I wrote to them and also tried to 
help Miss Straetmans’ efforts to get them released. 
Mrs. Dellinger and I sent them food parcels, a com- 
mon practice in'the war years. I also wrote to the 
chairmen of various National Committees urging 
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them to keep National Committee work going if 
possible. The spark of life was not extinguished. 
A General Assembly was held in Paris in September 
1946, a year after the end of the war. 

I wish there were time to tell you of many events 
and many people that have been in the forefront 
these 50 years. I will give you merely the titles of 
two great events, the Second Polar Year (1932 to 
1933) and the International Geophysical Year 
(1957 to 1958). The URSI is now recognized as the 
natural leader of great enterprises like the IGY. 
The fine work done by United States research people, 
and may I say especially you folks here in Boulder, 
is a major ingredient in our present prestige. 
Another ingredient is the man who served as presi- 
dent of the URSI the past three years and as presi- 
dent of the International Council of Scientific Unions 
just before that, Lloyd V. Berkner. We who belong, 
or have belonged, to the National Bureau of Stand- 
ards radio work, can take pride in the fact that he 
was one of us for several years (this was some 30 
years ago). 

The prestige and the glory of URSI mainly derive 
from the fact that it really accomplishes its objec- 
tive: the worldwide coordination of radio research. 
It is the place where radio science presents its find- 
ings. If you include the meetings of the national 
sections as well as the General Assemblies, it can be 
said that in the URSI meetings you hear the first 
announcement of everything new in radio science. 
Sometimes a great advance is not recognized as such 
when first presented. To mention one example: 
Karl Jansky described at a U.S—URSI meeting in 
1932 how he had received extraterrestrial radio waves 
and from their diurnal variation in direction identi- 
fied their source in the Milky Way. We did not 
know we were witnesses to the start of an extraor- 
dinary new branch of human knowledge, radio 
astronomy. 

I said that all of us are in on a good thing. It’s 
not only because we are privileged to be the pioneers 
but because of the tremendous scope of our enter- 
prise. The word “radio” does not tell the whole 
story. I am one of those who have kept the IRE 
from changing its name as more and more fields were 
added. The new fields do use radio techniques, so 
there is justification, even wisdom, in not changing 
the name. As a result, in IRE as well as in URSI, 
we have electronics, and what a field that is! Did 
you ever go through a modern airplane factory and 
see what there is between the gleaming outside and 
the glistening interior of the cabin? What you see 
is an amazing complex of the wires and fittings and 
devices of electronics. W ell, URSI also has cyber- 
netics and information theory and computers, and 
space communication and telemetry. Not only will 
radio with its branches, electronics and all, become 
the greatest of the industries, it will absorb most 
industry. I have been informed that electronics in 
the U.S. is already a bigger industry than agriculture. 

The scientific side of all this belongs to URSI, 
and everybody knows that without the science these 
things would not exist. I am not implying that 





URSI is a perfect organization. It is an evolving 
organization, changing as the needs of radio science 
change and as more men enter it. There is ample 
opportunity ahead, immediately and also later, for 
creative work in adapting URSI to this changing 
world. 

It seems likely that the next great step forward in 
science as a whole will come from our field, radio. 
I am thinking of what radio astronomy is doing and 
is leading to. I do not think men are ever going 
to go out among the stars in spaceships. (We shall 
probably go all around the solar system; not out 
among the stars: other creatures, not men, may do 
even that.) But we are going out on radio waves; 
and with our radio telescopes farther than with 
optical telescopes. 

At least one major and respectable radio astronomy 
observatory has a regular listening schedule for 
possible signals from intelligent beings somewhere 
in the universe. Many scientists (I among them) 
believe that there are myriads of planets of stars 

other than our sun on which there is life. Probably 
very many have life of more advanced character 
than mankind. Probably many are sending out 
signals with the idea of contacting life on other 
worlds. I have heard them credited with the idea 
of trying to contact us—No! Although intelligent, 
they never heard of us. Until very recently most 
of us had only an amused tolerance for any such 
ideas. By the way, the National Bureau of 
Standards continues to receive every once in a while 
at) inqury from somebody (a different person each 
tine) about alleged signals, described in meticulous 
detail, which were supposed to have been recorded 
at the Bureau in 1924 from Mars. That was a 
false rumor. But scientifically organized efforts to 
receive intelligent signals from the depths of space 


may very well be rewarded some time. It may 
happen some centuries hence, it could happen 
tomorrow. 

Another profound effort of the largest radio 


observatories is to reach the farthest observable 
parts of the universe. The universe is commonly 
thought of as expanding, constantly carrying galaxies 
out beyond the distance where their radiations can 
ever reach us, so no improvements in radio or other 
telescopes can ever observe them. Improvements 
have already taken us more than half way out to 
that limit. We are on the threshold of finding out 
whether the universe had a beginning or is being 
continuously created. If it had a beginning, the 
number of galaxies in a given volume of space was 
(according to the expansion theory) much greater 
near the beginning than now. Well, when we look 
out upon very distant galaxies, because of the finite 
velocity of light we see them as they were nearer 
the beginning. At a distance of several billion 
light-years, then, there should be more galaxies than 
in the same volume of space nearer us. But if, on 
the other hand, creation is and always has been 
going on, new galaxies are always being formed and 
their average number in a given space is presumably 
the same throughout the universe. Evidence bearing 
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on this is obtained slowly and with difficulty. It is 
one of the direct concerns of URSI and was one of 
the subjects of the General Assembly at London 
this year. 

The cosmological problem can be attacked by 
radio in another way. If creation is continuous, 
the first step in the process is the emergence (or 
creation) of hydrogen atoms at a slow rate through- 
out space. Since these radiate on a frequency of 
1,420 Me/s, it is quite possible that the presence of 
such hydrogen in space between the galaxies will be 
observed by radio telescopes. 

The implications of an imminent solution of the 
cosmological problem stagger the imagination. They 
could lead to a step in the understanding of the 
universe comparable to those taken by Maxwell 
and by Einstein. 

In conclusion, we tread a giant stage, don’t we? 
And don’t we have great fun doing it? The URSI 
meetings let us come out from our laboratories, and 





talk and listen. Lest we take ourselves too seriously 
in all this, I’d like to close with a quotation from a 
great astronomer who had some ideas that should 
help. He was the most noted of Persian astronomers. 
Here are three of the rubais of Omar Khayyam: 


““Myself when young did eagerly frequent 

Doctor and Saint, and heard great argument 
About it and about; but evermore 

Came out by the same Door where in I went. 


‘‘With them the seed of Wisdom did I sow, 

And with mine own hand wrought to make it grow; 
And this was all the Harvest that I reap’d— 

‘I came like Water, and like Wind I go.’ 


“The Revelations of Devout and Learned 
Who rose before us, and as Prophets burn’d, 

Are all, but Stories which, awoke from Sleep 
They told their comrades, and to Sleep returned.” 


(Paper 65D4-134) 
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Power Density Requirements for Airglow Excitation 
by Gyrowaves ' 
V. A. Bailey 


(December 9, 1960; revised December 27, 1960) 


The original proposal in 1938 for generating an artificial airglow by means of concen- 
trated and powerful gyrowaves is reconsidered and new estimates are made of the powers 
P radiated from an aerial array of D ideal dipoles which would suffice to enhance the normal 


brightness of the airglow by specified factors. 


Other interesting phenomena are mentioned 


which may be produced in the nocturnal E-region by means of such powerful beams of 


gyrowaves. 


Some 22 years ago [Bailey, 1938a and b] it was 
shown that it is possible to generate in the lower 
E-region of the upper atmosphere at night an arti- 
ficial airglow (or “aurora”’) by means of a concen- 
trated beam of powerful gyrowaves, i.e., waves of the 
same frequency as the frequency of gyration of a free 
electron in the local terrestrial magnetic field H, 
namely He/2rm. This frequency is about 1,300 
ke/s above London and 1,530 ke/s above Armidale 
(New South Wales) where experiments which bear 
on this subject are in progress. 

With an appropriate array of 800 horizontal dipoles 
the mean power required to generate an airglow 50 
times as bright as the night sky was estimated to be 
500 kw. Also, in a part of the region where the glow 
is generated the number-density of the electrons was 
estimated to exceed 50 times the normal density. 

The estimated brightness was based on the value 
16 v/em/mm Hg for the reduced electric force X/p 
found experimentally in the uniform column of an 
electrical discharge in air at a low pressure p, on the 
assumed pressure 2107? mm Hg at the 92-km level 
of the ionosphere and on an assumed luminous output 
of 1 candle/watt (c/w). 

New estimates have now been made which are 
based on the following facts: 

(1) In certain glow discharges in nitrogen the 
reduced “normal” electric force X,,/p is equal to, or 
less than, 3.25 v/em/mm Hg [Giintherschulze and 
Betz, 1933]. 

(2) In the uniform column of a discharge in nitro- 
gen the value of X/p approaches 4 as the tube radius 
or pressure is increased [Holm, 1923]. 

(3) From the uniform column of a discharge in 
nitrogen at the pressure 0.5 mm and with X/p=10 
the luminous output L is 0.5 ¢/w [Bailey, 1938a]. 

(4) In nitrogen the mean energies V,, of electrons 
with X/p=3.25 and 10 are respectively 1.37 and 1.80 
v [Townsend and Bailey, 1921]. 

(5) The resonance potential V, of the nitrogen 
molecule is 6.1 v. 

(6) Electrons in air have a mean energy of 1.37 v 
when X/p=4 [Crompton, Huxley, and Sutton, 1953]. 


! Contribution from Department of Physics, University of Sydney, Sydney, 
Australia. 





From the facts (3), (4), and (5) and the assump- 
tions that when X/p>3.5 the energies V of the elec- 
trons in nitrogen have a Maxwellian distribution? 
and the probability of excitation of light by a collid- 
ing electron is proportional to V—V,, we can deduce 
that with V,,=1.37 v the luminous output Z of a 
discharge in nitrogen is 0.08 ¢/w. 

Since 80 percent of air consists of nitrogen we may 
assume, with little error, that when the electrons in 
air have the same mean energy V,,=1.37 the 
luminous output is the same as in nitrogen. Hence, 
by (6), when X/p=4 in air containing free electrons 
then L=0.08 e/w. 

Since (1), (2), and (4) show that a glow discharge 
in nitrogen can be maintained when V,,=1.37 it 
follows from (6) that a glow discharge can be main- 
tained in air when X/p=4 and free electrons are 
supplied, by ionization by collision or otherwise, at 
a rate equal to the rate of loss of electrons by attach- 
ment, diffusion, and recombination.* 

We may therefore conclude that a glow discharge 
in the ionized air of the nocturnal lower E-region 
can be maintained by a circularly polarized gyro- 
wave which sets up a local rotating electric field X 
such that X/p=4, where p is the air pressure at 
the 92-km level corrected to 15°C. From rocket 
data we deduce that this value of X is equal to 
6.81073 v/em. 

In order to generate this gyrowave by means of a 
radio transmitter on the ground connected to an 
aerial array of D ideal dipoles the required mean 
radiated power P is given by 

P=8.64X10° X?/D kw, 
where _X is in v/em, i.e., by 
P=40,000/D kw. (1) 


2 Maxwell’s distribution is adopted here on account of the discussions of the 
energy distribution which are given in von Engel [1955], Bayet [1958], and Francis 
1960]. Since it relates only to perfectly elastic collisions, Druyvesteyn’s distri- 
ution is strictly not applicable to electrons in gases like nitrogen and air under 
the influence of the high-frequency electric forces considered here; for these forces 
are such that the average loss of energy by an electron at a collision is more than 
50 times as large as that lost in an elastic collision. Also several publications by 
M. Bayet et al., stress the fact that in many experiments involving relatively 
large electron temperatures (of the order of 30,000 °K) the energy distribution is 
found to be Maxwellian. For all these reasons it does not seem worthwhile to 
widen the present discussion to include Druyvesteyn’s distribution, 

3 This conclusion is in agreement with some observations made by M. Cutolo 
on low pressure discharges produced in air in a magnetic field by means o/ the 
corresponding gyrowaves (Private communication). 
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When D>80 the power flux density, cX?/4z, into 
the H-region overhead is uniform within 10 percent 
over an area which subtends at the ground a solid 
angle of 2/4D steradians. Within this area the 
power flux density for X=6.8X10-* v/em is 1.22 
1077 w/em?, and so the apparent brightness of this 
area would be about 10~* candles/cm’; this is of the 
same order of magnitude as the normal brightness 
of the night sky. 

Thus the gyrowave radiated with the power P 
given by eq (1) would enhance the normal brightness 
by about 100 percent; this is more than 10 times 
as large as the least enhancement which can be 
observed visually or by means of a photoelectric 
device. In particular we may conclude that an 
array of 80 dipoles radiating a gyrowave with a 
mean power of 500 kw would produce an easily 
observable enhancement of the nocturnal airglow 
overhead within a solid angle of 32 square degrees. 

If the transmitter could be made to generate 
pulses of gyrowaves each of power 3,500 kw and of 
length about 200 usec the corresponding value of 
X/p would be 10.5 and then each resulting flash of 
airglow would be about 40 times as bright as the 
night sky. 

A constant mean power of 500 kw radiated from 
an array of 500 dipoles would produce an equally 
bright steady glow within a solid angle of about 5 
square degrees. 

Since the existence of the normal airglow and the 
probable existence of notable ionizing streams of 
meteors in the #-region suggest that an appreciable 
fraction of the molecules present may be in excited 
states, the values of brightness given above may be 
notably exceeded or else attained with lower values 
of the product PD. 

Other interesting phenomena which may be pro- 
duced in the nocturnal E-region by such beams of 
gyrowaves, with similar or larger powers [Bailey, 
1959], are as follows: 

(1) Generation or enhancement of infrared radi- 
ation and of the radiation from sodium atoms. 

(2) Making visible parts of overhead meteor trails 
or the enhancement of such trails when visible. 

(3) Enhancement of radar echoes from such trails. 

(4) Making visible parts of overhead tracks of 
solar corpuscles or other extraterrestrial particles. 

(5) Production of transient local variations of the 
terrestrial magnetic field. 

The magnitudes of most of these phenomena are 
not easy to predict, at least because some of the neces- 
sary physical data are unavailable. But some sup- 
port for the conclusions (2) and (4) is provided by the 
experiments of A. R. Bevan [1949] in which the 
tracks of a- and B- particles through neon at 200 
mm pressure were revealed by the localized glows 
caused by strong microwave pulses. 

Even with less powerful beams of gyrowaves it is 
possible to make important, though invisible, changes 
in the nocturnal #-region [Bailey and Goldstein, 
1958; Bailey, 1959]. For example, with an aerial 





array of 40 dipoles and a power of 500 kw the corre- 
sponding value of X/p at 92-km height would be 
about 2.5 and so a pulse of length 2 sec. would nota- 
bly reduce the rate of attachment of local electrons 
to molecules and thereby increase the local electron 


density by a factor of about 5. The theory of this 
phenomenon also shows that such a pulse can be used 
to determine experimentally the rates of attachment 
of electrons near the 92-km level as well as other in- 
formation about the collisions of these electrons with 
neutral particles. Preparations are under way at the 
University of New England in Armidale to carry out 
experiments with such strong pulses and the use of 
auxiliary pulses for probing the ionospheric regions 
modified by the strong pulses. 

All these conclusions, together with the work pre- 
viously accomplished on gyrointeraction, show that 
experiments with beams of powerful gyrowaves are 
likely to throw light on several obscure aspects of the 
ionosphere such as the nature of the neutral and 
charged constituents of the H-region, the processes 
which determine the density of ionospheric electrons 
and some of the processes which cause the natural 
airglow. Moreover, the possibility of artificially in- 
creasing the local electron density in this way may 
have applications to radio communication. 


The author is indebted to Dr. R. A. Smith of the 
University of New England for his help with several 
aspects of this work. 
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On the Validity of Some Approximations to the 
Appleton-Hartree Formula’ 


Kenneth Davies and G. A. M. King ” 


(December 9, 1960) 


The validity of some commonly used quasi-transverse and quasi-longitudinal approxi- 
mations to the Appleton magneto-ionic formula is considered. Using the dipole approxi- 
mation for the earth’s magnetic field the various approximations for the refractive index 
are compared with the values computed from the complete formula for various geomagnetic 
latitudes and a frequency of 2.0 megacycles per second. It is found that certain approxi- 
mations become very poor only a short distance from where they are exact and so care must 
be taken in their use. It is shown that a choice of two suitable approximations yields refrac- 
tive indices of sufficient accuracy for all geomagnetic latitudes. Certain approximations to 


the group refractive indices are also considered. 


1. Intreduction 


Approximations to the Appleton-Hartree formula 
have been considered by a number of workers [Booker, 
1935; Rydbeck, 1940; Westfold, 1951; Ratcliffe, 
1959; Titheridge, 1959]. In most cases, however, 
the range of validity of the expressions used is dif- 
ficult to assess. Although, with the advent of 
electronic computers, the necessity for using approx- 
imations in numerical work has decreased; there is 
still some need for them in practical work and also 
they are widely used to simplify theoretical discus- 
sions. We shall consider here the applicability of 
some approximate formulas for the case of the 
earth’s dipole magnetic field and for frequencies 
normally used in ionospheric sounding at vertical 
incidence. Calculations involving the time of flight 
of radio pulses require the group refractive index 
and so we have included a consideration of some 
approximations to the formula for this quantity. 
(The effect of electronic collisions on the phase and 
group indices has been neglected.) 


2. Definitions 


N=electron density 
f= frequency 
fu=gyrofrequency 
fy=plasmafrequency 
6=angle between direction of propagation and the 
magnetic field 





IwY 
P= f— ffu 





! Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 

2 Present address, Geophysical Observatory, P.O. Box 2111, Christchurch, 
New Zealand. 





3. Approximations 


The expression for the phase refractive index in 
the ordinary mode, n+, is 
/y* sin* 6 ‘ 


a aaaling 


zs. . y* sin? 0 
2(1—2) 


1—n 





The quasi-transverse approximations, n; and m2, are 
given by 
ar 


1—n}? 





=1, 





ian t+a-2) cot? 6, 


and the quasi-longitudinal approximations are given 
by n3, 24, and ns; 
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The expression for the phase refractive index in 
the extraordinary mode, 7_, is 


e. ., are Ve ae 
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The quasi-transverse approximations are given by ng, 
No, and 7, and the quasi-longitudinal approximations 
are given by ng, 27, and Np. 


ing 


323 




















£ 
=1— 
ia 1—|y cos 6}, 
z _, y sin’ 6 
1—nj l—z 
OO. 
ia 1—2’ 
ane )_y? sin? 6 
a a ly cos | 2(1—2)’ 
x y? sin? 6 
neste aoe reno t= gf) COREG. 
ee [eee (1—z) cot? 6 


The group refractive index (n’ =n+f *) is given 
by the lengthy formula, : 


y’ sin? 6 


sited ina" meta (1—2)? 





2(1—2)* y? cos? ¢@— xy! sin‘ 6 } 
as ; . , 
2(nz)a(1—z)?Vy* sin’ 6+4y?(1— x)? cos? 6 








and we discuss these approximations to it: 
Ordinary Mode, Quasi-transverse Approximation 


(n2)’ = (ne) 


4 x(1+ cot? @) 
(nz) (1+ cot? 6) {1+ (1—2z) cot? 6}+ (n2)x cot* 6 





Extraordinary Mode, Quasi-longitudinal Approxima- 
tion 
22—xy 


(1)’=("s) +356) G—? 


Extraordinary Mode, Quasi-transverse Approximation 


a(1—z)?+2%y? 
Ng) (l—z—y?)? 





(7) ‘= (Ns) +7 


4. Discussion 


For a dipole field, the propagation angle, 6, and 
the geomagnetic latitude, , are related for vertical 
incidence propagation by 

cot @=2 tan dA. (See table 1) 
7 herefore, the discussion will be carried on in terms 
of i. 

(a) Phase index in the ordinary mode. In figures 
la, b, c, d, e, and f we have plotted the variations of 
n, and the approximations n, to n; as functions of 
geomagnetic latitude (the propagation angle, 0, is 
also indicated). The graphs are for a frequency, f, 
of 2 Me/s, and each one treats a different plasma- 
frequency between 0.5 f and 0.998 f. We also exam- 








ined the cases of f=4 Mc/s and f=8 Me/s, and 
found that all the approximations held more closely 
than at 2 Me/s, but they maintained their relative 
order. 

The best approximation to n, appears to be a 
combination of nm. and n;. At lower latitudes nz. is 
the more accurate with a positive error, and at 
higher latitudes n, is the more accurate with a nega- 
tive error. The latitudes at which these errors are 
equal in magnitude together with the percent errors 
are given in table 2. It is of interest to note that, 
where the errors are largest (fy/f=0.998) they only 
apply to a very narrow range of X. Figures 2a and 
2b which show n as a function of plasma frequency 
for f=2.00 Me/s and latitudes of 60° and 30° respec- 
tively confirm the above conclusions. Table 3 
gives the latitude ranges over which the various 
approximations are within 5 percent of the true 
value. 

The most commonly used approximation in the 
ordinary mode is 3. It is interesting to note that 
it holds well only at very high latitudes and, almost 
fortuitously, at the equator. (One should not use 


TABLE 1. Variation of propagation angle, 0, and gyro- 
frequency with geomagnetic latitude for a dipole field 


cot 9=2 tand 


Surv=fu0 V1+3 sin? 











n° 6 Su 

Me/s 
0 90° 00 0. 8000 
5 80° 04 . 8080 
10 70° 35 . 8354 
15 61° 49 . 8763 
20 54° 57 . 9295 
25 47° 00 . 9895 
30 40° 54 1. 0583 
35 35° 32 1. 1266 
40 30° 47 1. 1974 
45 26° 34 1. 2649 
50 22° 46 1. 3290 
55 19° 19 1. 3888 
60 16° 06 1. 4422 
65 14° 08 1. 4890 
70 10° 19 1. 5284 
75 07° 38 1, 5594 
80 05° 02 1. 5819 
85 02° 30 1. 5952 
90 00° 00 1. 























TaBLE 2. Point of equal error in ny and ny from ny 
in nN Error 
i 

deg % 
0. 500 12 0.6 
. 600 13 8 
- 700 15 1.5 
. 800 22 37 
. 900 33 es 
. 930 41 9.1 
. 960 52 11.7 
. 980 62 11.8 
. 990 7 12.0 
995 76 14.7 
.998 | 82 8.5 








For 0.998, the data are insufficient to accurately plot m4; the figures shown 
are the minimum values. 
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TaBLE 3. Areas of better than 5 percent accuracy in TaBLE 4. Areas of closer than 5 percent accuracy in 
approximations to n+ approximations to n- 
\=Geomagnetic latitude (deg) A= Geomagnetic lz atitude (deg) 
2 | | 4 fy | 
fn ; ™ ma 3 m4 ms — ne nz ns ng n10 
J s 
_ + + _ 0. 500 05590 0S\ S90 A<29. | A<56 0S\<90 
0. 500 A<41 o 37 0S 590 0OSAS90 | 0SA S90 . 600 12<A 0SA S90 A<19 . 0SdA3890 
600 A<25 <30 05590 0S =90 | 13<d . 700 23<A 45<A A<14 A<19 0S\A 390 
. 700 A<19 Noo 0S\ S90 04<X ; 2a - 800 31<A 61<A A<08 08<A 
800 | A<I3 A<30 30<A 16<d | 41<d : 
<10 900 | 38<r 75<A <06 A<12 23>d 
- 900 A<12 A<32 52<A 35<A 55<A -930 | 388<A 79<A A<04 31 
<03 960 | 40<d 83< A<O1 A<08 41<d 
980 | 42<d 86<A <1 48< 
. 930 A<12 <34 61<A 45<X 63<A - 990 38<A 87<A | A<OL A<10 | = 55<A 
A<02 | . 
. 960 A\<12 A<43 69<X 56<A | 70<d .995 40<X 88<A | »<O1 | 61<d 
A<O1 | .998 | 37<A 89<A A<Ol | A<O 66<A 
980 | r<14 A<50 77<X | 66<A | 77<A 
, AOL | 
990 | A<15 A<66 82<A 75<A | 82<A Likewi : 1 i 
A<O1_ | ew 5 , > eGus there it 
oe ore ae dike | s2<a beads akewise, mg is good only at the equator where it 
= | degenerates into m9, a simpler formula. 
998 | x<22 | x<s0 88<d 84<r 88<) Plots of n versus fy/f are given in figures 5a and 5b 
| A=0 so that the above discussion can be considered in 





nz near the equator if there is the slightest doubt that 
the propagation angle is exactly zero.) In the lati- 
tude range where ns is satisfactory, the slightly 
simpler formula, Ns, is almost as good. 
(b) Group index m the ordinary “mode. 


The group 
refractive index, n’, is given by 


nant fG 


This equation can be applied to the approximate 
formulas for n to give expressions for n’. Clearly, 
it is not worth doing this for n; and n; since these 
are poor approximations to start with. Also nj is 
almost as complicated as n‘,, and is, therefore, of 
little interest. So, figures 3a, b, c, d, e, and f give 
a comparison of only mn’, and n3. 

As the region making the greatest contribution to 
the group path is where fy is almost equal to f, it 
seems wise to restrict the use of n5 to latitudes less 
than A=30°. 

(c) Phase index in the extraordinary mode. Curves 
showing the variations with latitude of n_ and the 
approximations ng to ny, are given in figures 4a, b, 
c, d, e, and f; the latitude ranges over which the 
various approximations hold to within 5 percent 
are set out in table 4. 

mo is the extraordinary analogue of nm, and n; is 
the analogue of no. 40 and n, are satisfactory only 
at high and at low latitudes respectively for small 
values of fy/é. A combination of nj and 1, is not 
satisfactory for values of fy/£ approaching unity as is 
the case for nz and nz in the ordinary mode. For 
larger values of fy/£, the choice seems to lie between 
nm, at high latitudes and no at low latitudes. It 
appears that approximations in the extraordinary 
mode are not really satisfactory at low-middle 
latitudes (10° to 30°). 

The approximation most commonly used for the 
extraordinary mode, nz, is noticeably poor except at 
the pole (and at very low values of twit). A discus- 
sion of ns and n; has been given by King [1960]. 


588261— 








another form. 

(d) Group index in the extraordinary mode. Fig- 
ures 6a, b, c, d, e, and f compares n_, ng, and no. 
ng is surprisingly good at least for latitudes above 
45°. mg holds exactly at the equator, but its use 
elsewher e is questionable; of course, at the equator, 
my degenerates into the transverse case. 


5. Conclusions 


Two formulas n, and n, are preferred for the phase 
index in the ordinary mode at high and low latitudes 
respectively. The latter gives rise to a suitable 
formula for the group index at low latitudes. 

The phase refractive index in the extraordinary 
mode is not so easily treated. However, here again 
two formulas, ns and m9, are probably adequate for 
many purposes. The most serious difficulties occur 
at low-middle latitudes. ng gives rise to a suitable 
formula for the group index at high latitudes. 

We wish to point out that there is a tendency in 
the literature to stress certain approximations; e.g., 
Mm, N3, and n7, and to use them outside their proper 
range of applicability. 


The authors acknowledge the help of Irving Litt- 
man who carried out many of the computations 
involved. 
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Figure 1. Comparisons of the geomagnetic latitude variation 


of the phase refractive index in the ordinary mode (n+) and 
various approximations to it. 
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Figure 2. Comparisons of the variations, with plasma frequency, of the phase refractive index in the ordinary mode (n+) 
various approximations to it at geomagnetic latitudes of 60° and 30°. 
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The latitudinal variation of the group refractive 
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JOURNAL OF RESEARCH of the National Bureau of Standards—D. Radio Propagation 
Vol. 65D, No. 4, July-August 1961 


Amplitude and Angular Scintillations of the Radio Source 
Cygnus-A Observed at Boulder, Colorado’ 


R. S. Lawrence, J. L. Jespersen, and R. C. Lamb 


(November 3, 1960) 


Variations in the apparent flux and position of the radio source Cygnus-A were recorded 
at 53 and 108 megacycles per second using a two-element, phase-sweeping interferometer 


located at Boulder, Colo. 


An ionospheric sounder operating at Ellsworth, Nebraska, 


provided, for a few hours each day, simultaneous vertical-incidence measurements on the 


ionosphere at its intersection with the line of sight from Boulder to the radio star. 


Amplitude 


scintillations observed at Boulder over a twelve-month period are compared with ionograms 


taken at Ellsworth. 


Positive correlation is found between amplitude scintillations and 


spread F, while no significant correlation is found with sporadic E. 
Detailed analysis of the scintillations indicates that the probability distribution of the 


amplitude can be represented by the Rice probability distribution function. 


The zenith- 


angle dependence of the amplitude scintillations does not agree with a theory based upon 


isotropic ionospheric inhomogeneities. 


The root-mean-square value of angular scintillations is proportional to the square of the 


wavelength, in accord with a theory of diffraction by ionospheric irregularities. 


Comparison 


of angular scintillations with amplitude scintillations indicates that, for elevation angles of 
15° to 50°, the region of the ionosphere responsible for scintillation lies near the border be- 
tween the Fresnel and Fraunhofer diffraction regions for both frequencies. 

Slow, irregular angular variations are commonly observed in the daytime at both 


frequencies. 
dimensions as large as 200 kilometers. 


1. Introduction 


It has been recognized since 1950 that the apparent 
fluctuations in intensity of discrete radio sources 
arise in the ionosphere. In a manner analogous to 
the twinkling and dancing of visible stars, these 
“amplitude scintillations” are accompanied by ‘an- 
gular scintillations”’, which are fluctuations in the 
apparent position of the radio source. 

Radio-star scintillations offer a convenient means 
for studying the irregular structure of the ionosphere. 
The F region is generally believed to play an impor- 
tant role in the production of scintillations, but our 
knowledge of the effect of the EF regio 1 is confused by 
incomplete and apparently conflicting evidence. One 
of the purposes of the experiment reported here was 
to improve our understanding of the ionospheric 
conditions and levels responsible for scintillations. 
The experiment was designed to be particularly 
sensitive to any relationship between scintillations 
and the E region. The results have indicated little 
or no such relationship, but do show a connection 


between scintillations and spread echoes in the F 


region. In addition to the comparison of scintilla- 
tions and ionospheric phenomena, detailed analysis 
has been made of the scintillations themselves. In 
connection with this, some statistical aspects of 
amplitude scintillations are found to disagree with a 
theory based upon plane-earth geometry and spher- 


1 Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 





These variations are attributed to lens-like ionospheric irregularities having 


ical irregularities whereas the statistical ,variations 
in angular scintillations agree with this theory. 

A third aspect of the experiment has been the 
observation and study of slow angular variations in 
the apparent position of the source. These varia- 
tions are due to refraction by individual large irregu- 
larities, in contrast to amplitude and angular scintil- 
lations which are due to the diffraction effects of a 
number of smaller irregularities. 

A description of the planning and instrumentation 
for the experiment, along with a brief summary of the 
pertinent literature, has already appeared [Lawrence, 
1958]. In addition, two review papers [Booker, 1958; 
Little et al., 1956] are available which outline present 
knowledge concerning scintillations and which de- 
scribe the usefulness of scintillation observations as a 
tool for studying the irregular structure of the 
ionosphere. 


2. Description of the Observations 


The radio source Cygnus-A has been observed at 
Boulder (latitude 40°7’.7 N, longitude 105°14’.3 W) 
at frequencies of 53 and 108 Mc/s, using a two-ele- 
ment, phase-sweeping interferometer. The phase- 
sweeping system permits simultaneous observation of 
the apparent flux from the radio source and of the 
relative phase of the radio waves at the two antennas. 
Amplitude scintillations are simply variations in the 
flux, while angular scintillations may be deduced from 
the known spacing of the antennas, the general direc- 
tion of the radio source, and the variations in relative 
phase of the radio waves. 
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Frem February 1958 through February 1959 ob- 
servations were obtained regularly during the three or 
four hours each day following the rising of Cygnus-A. 
Figure 1 shows how the time of such observations 
varies throughout the vear and indicates the seasons 
during which the effective ionosphere was sunlit. 

The two antennas are separated by 475 m on a 
horizontal, east-west baseline and the resulting inter- 
ferometer lobes are of nearly constant spacing over 
the region of the sky in which the observations were 
made. At a radiofrequency of 108 Me/s the lobe 
spacing is about 0.4 deg; at 53 Me/s it is twice as 
great. The measured component of angular scintil- 
lations is, of course, perpendicular to the lobes. This 
means that, near the horizon, angular scintillations 
can be observed only in the horizontal direction, 
while at an elevation angle of about 35° the horizontal 
and vertical components of angular scintillations are 
observed with equal sensitivity. 

At the time when Cygnus-A is about 15° above 
the north-east horizon at Boulder, the line of sight 
intersects the EF region over Ellsworth, Nebraska. 
During the appropriate portion of each day, an iono- 
spheric sounder was operated at Ellsworth at 5-min 
intervals. As the radio source rises, the intersection 
point of the line of sight with the region moves to- 
ward Boulder, and the intersection with the F region 
passes within about 100 km of Ellsworth. Figure 2 
is a map of the Boulder-Ellsworth area showing the 
plan view of the intersection point for various 1iono- 
spheric heights. The numbers along the curves indi- 
cate the local sidereal time (LST) at Boulder. 
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3. Instrumentation 


The instrumentation for the scintillation measure- 
ments has already been described in some detail 
[Lawrence, 1958], but is briefly reviewed here. The 
two frequencies, 53 and 108 Me/s, are observed with 
the same antennas and with similar receiving and 
recording systems. 

Antennas. The two 40-ft paraboloidial reflecting 
antennas are modifications of standard 28-ft parabo- 
loids manufactured by the D. S. Kennedy Company. 
The 6-ft radial extensions were designed and con- 
structed at NBS. The antennas, separated by 475 
m, are equatorially mounted and equipped with 
sidereal drives for tracking radio stars. The dipole 
feed, with four reflecting elements, is designed for use 
simultaneously at both frequencies. A multicoupler 
permits the use of a single coaxial transmission line 
from each antenna to the central receiving shelter. 
No preamplifiers are used at the antennas. 

Receivers. At each frequency, the dual receiving 
channels are connected together as a Ryle inter- 
ferometer [Ryle, 1952]. A rotary phase shifter has 
been added to one channel so that the interferometer 
pattern may be swept across the sky. This phase 
shifter, usually driven at a uniform rate of 10 rpm, 
converts the system to a “‘phase-sweeping”’ inter- 
ferometer. A ‘“‘phase-tracking’ feature was in- 
cluded in the original design, whereby the inter- 
ferometer pattern was servo-controlled to follow the 
apparent motion of the radio source. This feature 
permits the direct recording of angular scintillations. 
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Ficure 1. Time of observation of Cygnus-A, and state of illumination of the ionosphere 


334 








com 
were 
ings 
anal 
each 
the 1 
The 

of 1. 
mani 


Th 
proce 
papel 
plana 

Th 
exper 
phase 
Boulc 
sound 

Ion 
10NOs} 





vt 












@ NORTH PLATTE 








© Fort Couums 
--— 100 km —4 


= eo LO: 








FicureE 2. The variation, with local sidereal time at Boulder, 
of the position at which the line-of-sight from Boulder to 
Cygnus-A penetrates the ionosphere. 


The three curves apply to ionospheric heights of 90, 100, and 300 km. 


However, it keeps a zero-crossing of the inter- 
ferometer pattern always on the source and so 
prevents recording of amplitude information. As 
a result, it is difficult to distinguish authentic angular 
scintillations from the random behavior of the servo 
loop which accompanies deep amplitude fades. 
Accordingly, all of the observations reporied here 
were made in the phase-sweeping mode. 

The receivers have bandwidths of 500 ke/s, and 
were operated with output time constants of one 
second. The double-conversion superheterodyne 
system uses intermediate frequencies of 31.5 and 
7.5 Me/s. 

Digital recorder. In order to permit the use of a 
computer for processing the data, chart recordings 
were supplemented by simultaneous digital record- 
ings on punched paper tape. A simple, meter-type 
analog-to-digital converter sampled the output from 
each interferometer once per second, and punched 
the two outputs alternately on a single paper tape. 
The six-bit binary code provides a full-scale accuracy 
of 1.5 percent, approximately the same as careful 
manual scaling of a pen recording. 


4. Data Analysis 


The following is a summary of data-analysis 
procedures. The appendices and the following 
paper by the same authors, contain a detailed ex- 
planation of the procedures outlined in this section. 

The observations made in the progress of this 
experiment consist of scintillation records from the 
phase-switch, phase-sweep interferometer located at 
Boulder, and of ionograms from an _ ionospheric 
sounder located at Ellsworth, Nebr. 

Ionograms. Vertical-incidence, sweep-frequency 
jonospheric records were taken at 5 min intervals 





for 1 hr, 1230 through 1330 LST each day during 
the time that the ray path from Cygnus-A passed 
through the FE region (100 km) near Ellsworth. 
After this period, ionograms were made each half- 
a for the next 3 hr. These ionograms have been 

caled for critical frequencies, minimum virtual 
heiahiee, and for the occurrence of special features, 
particularly for the absence or presence of spread F 
and of the various types of sporadic #. Each type 
of information from the ionogram was converted to 
an appropriate single-digit scale and punched in a 
separate column of a uniquely identified IBM card. 
Simultaneous scintillation information, described 
below, was later added to each card. 

Scintillation records. The scintillation records, 
both analog and digital, consist of a sine wave (‘‘car- 
rier’) having a frequency determined by the rate 
of lobe sweeping and an amplitude and phase which 
are modulated by the ionosphere. For the present 
experiment several methods were used to extract, 
from the digital records, the amplitude and phase 
modulation. However, all of these methods depend 
upon the fact that the envelope of the carrier is 
affected only by ionospheric amplitude effects and 
that the times of zero-crossings are affected only 
by ionospheric phase effects. Intermediate points 


on the carrier are affected by both. The power 
index 
(4F)=| _|P—P| 
| P 


and the amplitude index 


(e2y -(454) 


have been computed to measure amplitude effects. 
P is the power from the discrete source and may be 
read directly from the envelope of the carrier since a 
square-law detector was used in the scintillation 
receivers. A=~+7P, and a bar denotes an average 
over the record. The reason for computing both 
of these quantities will be shown in a later section. 

The variance of the angular scintillations, 6°, was 
computed as a measure of phase effects. For any 
given elevation angle, 6? can be deduced from the 
observed value of oz, the standard deviation of the 
phase difference, a, between the two ray paths to 
the interferometer. Variations in a may be deter- 
mined directly from the variations in the interval 
of time between successive zero-crossings. (See 
following paper by the same authors.) 

The first attempt to measure amplitude scintilla- 
tions was to assign visually one of five possible cate- 
gories to each 5-min section of the scintillation 
records (see appendix 1). Those sections for which 
there were no visible variations in the envelope of 
the carrier were assigned a ‘“‘depth index” of 1, those 
for which there were shallow variations in the enve- 
lope were assigned a depth index of 2, etc., to the 
depth index 5 for records with amplitude effects so 
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deep that the lobes disappeared. Several examples 
of each of these categories of records were scaled 
by hand so that an estimate could be made of the 
range in values of (AA/A)? and of (AP/P) associated 
with each category (see the table in appendix 1). 

A second approach was to recover P and a from 
the carrier and then to plot these two quantities as 
functions of time. The resulting graphs are referred 
to as power plots and phase plots respectively. 
Two different methods were used (‘‘zero-crossing”’ 
and “‘filter’’?) to make this recovery (see following 
paper). The power and phase plots may be used to 
compute the quantities (AP/P), (AA/A)’, and 6, 
and it is also possible to use them to determine the 
power spectra of amplitude and angular scintillations. 

Some of the records, because of large amplitude 
and/or phase effects, are not amenable to analysis 
by either the “zero-crossing” or “filter”? methods. 
These records were analyzed using a statistical ap- 
proach which does not give detailed power and phase 
plots but which gives the quantities P, and oz, 
quantities which may be related to (AA/A)? and 
(see following paper by the same authors). 


5. Results and Discussion 


The relation between amplitude scintillations and | 


ionospheric parameters. 
bution with elevation angle of the Ellsworth iono- 
grams and the coincident scintillation records which 
were used to investigate the correlation between the 
two types of observations. Figure 4 shows a definite 
correlation at both radiofrequencies between the 
depth of amplitude scintillation and the occurrence 
of spread F. Here the Fisher Z, function [Walker 
and Lev, 1953] is plotted at the weighted mean 
elevation angles of the 10-deg class intervals shown 
in figure 3. The upper scale gives, for each eleva- 
tion angle, the ionospheric height at which the line 
of sight to Cygnus-A comes closest to the zenith at 
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Figure 3. The grouping for statistical analysis, with elevation 
angle and sidereal time, of the scintillation data. 


The “effective ionospheric height” is the height at which the line-of-sight 
from Boulder to Cygnus-A comes closest to the zenith at Ellsworth, Nebr. 


Figure 3 displays the distri- | 


CORRELATION COEFFICIENT, Zr 





Ellsworth. N is the number of samples used in the 
computation of each value of Z,, and the vertical 
bar through each point has length 20,, twice the 


| expected standard deviation of the correlation co- 


efficient. Each graph is derived using the total data 
sample. Division of the data into sunlit and dark 
ionospheric groups results in correlations (not shown) 
which have, in general, a value of Z, which is 0.1 
higher (dark ionosphere) and 0.1 lower (sunlit iono- 
sphere) than the correlations shown, with corres- 
pondingly larger values of o,.. 
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Ficure 4. The correlation between amplitude scintillations 
and spread F, as a function of source elevation for each 
radio frequency. 
N is the number of cases. 


336 





Fe 


CORRELATION COFECICIONT 


Fr 











In the same manner, figure 5 displays the correla- 
tion, using the total data ‘sample, between the depth 
of amplitude scintillations and the occurrence of 
sporadic £. Division of the data according to sunlit 
and dark ionosphere, as well as by season, resulted 
in no significant change from the results shown. 

Other ionospheric parameters showed little or no 
correlation with the depth of amplitude scintilla- 


tions. In particular, there was no significant corre- 
lation, at either radiofrequency, with foF2, foEs, 


h’Es, or fain. Correlation with h’F, the minimum 
virtual height of the F region, resulted in a Z, of 
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Ficure 5. The correlation between amplitude scintillations 
and sporadic E, as a function of source elevation for each | 
radio frequency. 
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about 0.2 to 0.3 at 53 Me/s and a Z,of about 0.1 to 
0.2 at 108 Me/s, at F-region effective heights. 

Many, but not all, previous observers have found 
scintillations to be correlated with the phenomenon 
of spread-F ionospheric echoes. In the early days 
of scintillation study, Ryle and Hewish [1950] and 
Little and Maxwell [1951] both found sail relation- 
ships, while Mills and Thomas [1951] found positive 
correlation with an ‘‘F-region index’”’. Later Bolton, 
Slee, and Stanley [1953] reported positive correla- 
tion during the winter months; as did Wild and 
Roberts [1956a] for nighttime scintillations, although 
in both cases the ionosonde was located many 
hundreds of kilometers from the point of F-region 
penetration. Both Duefio [1955] and Hartz [1958] 
failed to find correlation with ionosonde data taken 
at points appreciably far from the optimum for F- 
region penetration, although Hartz and, more re- 
cently, Dagg [1957] noted that the diurnal variations 
are similar. Briggs [1958], on the other hand, has 
recently studied the correlation between scintilla- 
tions and spread F to deduce the spatial extent and 
the height of the region responsible for scintillations. 
Wright, Koster, and Skinner [1956]; and Koster, 
[1958] have noticed a quite high positive correlation 
with equatorial spread F, but this may be a differ- 
ent phenomenon from the spread F' at temperate 
latitudes. 

The diurnal variations of spread F and scintilla- 
tions are generally reported to be similar, and 
comparison between figures 6 and 7 indicates that 
this is clearly the case for the present experiment. 
This fact alone might explain some of the correla- 
tions which have been reported. In other cases, a 
disadvantageous location of the ionosonde may have 
reduced or destroyed the correlation. We have at- 
tempted to avoid both pitfalls, first, by removing 
the average diurnal and seasonal variation from both 
the scintillation data and the spread F’ observations 
(see appendix 2), and second, by locating the iono- 
sonde at Ellsworth, Nebr. Although Ellsworth was 
chosen specifically for H-region observations, the 
location should provide quite reliable F-region data 
for elevation angles of Cygnus-A in the range from 
25 to 65 deg above the northeastern horizon at 
Boulder. 

The positive correlation between spread F and 
the depth of amplitude scintillations, as shown by 
figure 4, agrees with many of the results described 
previously. The location of the ionosonde, although 
not ideal for F-region effects, is quite satisfactory if, 
as has been reported [Briggs, 1958; Koster, 1958; 
Yeh and Swenson, 1959; and Reid, 1957], the irregu- 
larities responsible for spread F occur in patches 300 
to 500 km in lateral extent. The data revision proc- 
ess, Which removed the diurnal and seasonal varia- 
tion from the ionospheric and scintillation param- 
eters, insures that any correlation which remains 
must be the result of detailed hour-to-hour or day- 
to-day correlation between the two variables (see 
appendix 2). 

The correlation, for a given elevation angle of 
Cygnus-A was generally higher at 53 Mc/s than at 
108 Me/s. This result might have been anticipated 
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since the refractivity of the ionosphere varies in- 
versely with the square of the radio-wave frequency. 
The nearly constant value of correlation, at either 
radiofrequency, over the whole range of elevation 
angles, indicates that the patches of spread F cover 
large areas. This result has the unfortunate effect 
of making height determination, by the method used 
in this experiment, impractical. The correlation 
was considerably higher at nighttime than in the 
daytime, and this result tends to reinforce the idea 
that scintillations are related to spread F, since 
spread F’ reaches a maximum shortly after local 
midnight. It is possible that this variation of cor- 
relation with local time is partly a result of our 
inability to remove completely the effect of diurnal 
and seasonal trends. It is clear that correlations of 
this type could be found simply by comparing trends, 
even though there is little, or no, detailed correlation 
present. 

The correlation between scintillations and spread 
F is significant, but still far from perfect; however, 
two features of our experimental technique have 
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Ficure 7. Hourly average values of scintillation depth index 
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and to the ‘‘noise-to-signal ratio”, k. 
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tended to reduce the apparent correlation. First, 
Ellsworth is not ideally located for F-region com- 
parisons; second, and perhaps more important, we 
have used a spread-F index which depends only 
upon the absence or presence of spread F’, not upon 
its intensity. By themselves, these two effects 
would be insufficient to reduce an otherwise perfect 
correlation to the observed value of about 0.3. 
Therefore, it seems probable that scintillations 
sometimes arise in ionospheric regions, such as F# or 
upper F, which do not produce observable spread-F’ 
echoes. 

A number of scintillation observers have searched 
for correlation with E-region phenomena. Mills 
and Thomas [1951] found no correlation with 
sporadic #, whereas Duefio [1955] and Bolton, Slee, 
and Stanley [1953] reported a general, not detailed, 
positive correlation between the occurrence of 
sporadic / and the occurrence of scintillations. In 
each case, the ionosonde was located more than 200 
km from the /-region penetration point. Wild and 
Roberts [1956a] found positive correlation with day- 
time fEs, but no correlation at night. Hartz [1958] 
working in Canada, found a general negative corre- 
lation with sporadic F. 

Because sporadic EF seems to occur in smaller 
patches than does spread F, it is significant that the 
Ellsworth station was located at precisely the correct 
place to study correlations at £-region heights. 
The ray path from Cygnus-A to Boulder passes very 
near the zenith at Ellsworth at the 100 km level and 
is within 100 km of the zenith for elevation angles in 
the range from 10 to 22 deg above the northeastern 
horizon at Boulder. As mentioned previously, 
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the data revision process retains any correlation 
between detailed events but markedly reduces the 
effect of diurnal and seasonal trends. The results 
of other workers [Hartz, 1958; Dyce, 1955] as well 
as those shown in figures 6 and 7 indicate that the 
diurnal variation of the occurrence of sporadic 
and the diurnal variation of the occurrence of scintil- 
lations are inversely correlated. 

Interim results [Lawrence and Jespersen, 1959] of 
the present experiment indicated a negative correla- 
tion (at H-region heights) between scintillations and 
sporadic FE. Those early results were obtained from 
only the first six months (Feb. 58—July 58) data, 
and the diurnal and seasonal trends had not been 
removed. It appears in retrospect that these corre- 
lations were produced by the inverse diurnal trends 
in the two variables. They, as well as the present 
results, are derived from data taken at 5-min 
intervals in the relevant range of elevation angles. 
Accordingly, the sensitivity of the result to any real 
relationship should be high. 

The present results, shown in figure 5, indicate no 
significant correlation, at either radiofrequency, with 
sporadic #. Although this contradicts the weak 
negative correlation mentioned in the interim report, 
we infer that there is no observable, detailed con- 
nection between sporadic # and amplitude scintilla- 
tions at this latitude. 

Excepting spread F, the only ionospheric param- 
eter to exhibit significant correlation with scintilla- 
tions was h’F, the minimum virtual height of the 
F region. This positive correlation was higher at 
53 Me/s than at 108 Me/s and tended, even more 
than did the correlation with spread F, to be re- 
stricted to a range of elevation angles which indicated 
F-region heights. This suggests that the properties 
of the electron-density profile which determine h’F 
are related to scintillations, but do not have good 
serial correlation over such large distances as does 
spread F. In addition to this, the correlation may 
be related to the diurnal and seasonal variation of 
radio star scintillations, a variation which has been 
the object of much study [Bolton, Slee, and Stanley, 
1953; Hartz, 1958; Dagg, 1957; Briggs, 1958; 
Koster, 1958; Yeh and Swenson, 1959]. Figure 7 
shows the observed variation in (AA/A)?. This 
variation contains a mixture of both seasonal and 
diurnal effects since only one source, Cygnus-A, was 
used. Theory indicates [Booker, 1958] that, in the 
Fresnel (near zone) region, (AA/A)? is proportional 
to the square of the distance, z, between the observer 
and the diffracting screen (see eq (2)). If h’F is 
related to z, then one would expect a correlation 
between amplitude scintillations and h’F; also, any 
seasonal or diurnal variation in h’F’ would produce 
a correponding variation in (AA/A)*. However, 
hmaxf’ is probably a better measure of z than is 
h’F. In view of this possibility, the observed diurnal 
and seasonal variation in hmax/’, at this latitude for 
the period corresponding to the times that scintilla- 
tion data were taken, was used to obtain 2 for eq (2), 
and the resulting variation in (AA/A)? was com- 
pared with the observed variation (fig. 7). At 53 
Me/s (using the P,, curve) the agreement was quite 








good. No comparison was made at 108 Me/s since 
the available portion of the P, curve shows little 
variation. The agreement at 53 Mc/s shows that 
a hypothetical diffracting screen which always re- 
mains at the height of maximum electron density, 
and whose statistical properties never change, will 
produce the same diurnal variation in amplitude 
scintillations as is observed at Boulder. However, 
such a model seems quite unreasonable since it 
implies that all the irregularities are forced to move 
with Amax/’, and that their intensity, AN, is inde- 
pendent of electron density, N. 

Amplitude scintillations. The power plots ob- 
tained from the ‘“‘zero-crossing”’ and “‘filter’’ methods 
provide a convenient basis for detailed analysis 
of amplitude scintillations. Using the method 
described by Tukey [1949], power spectra extending 
over the frequency range 0.005 to 0.17 c/s have 
been calculated from the power plots. Four such 
spectra, together with the corresponding values of 
(AA/A)?, are displayed in figure 8. The ordinate 
value is the common logarithm of the variance per 
unit frequency interval, expressed in units of the 
square of the average power. Since the spectra 
come from power plots, and these cannot be ob- 
tained for the highly disturbed records, the two 
upper curves in figure 8, although they are the 
highest spectra available, should be taken as typical 
of only moderate scintillation activity. The two 
lower spectra were derived from data with little 
scintillation activity, and are close to the noise 
level. wd 
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Fraure 8. The figure shows the highest and lowest observed 
amplitude power spectra for each radio frequency. 
The lower level at each frequency indicates the background ‘“‘noise’’ level. 
The units of W are the variance per unit frequency interval expressed in units of 
the square of the average power. 
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Each pair of spectra in figure 8 was derived from 
simultaneous power plots at the two radiofrequen- 
cies. Comparison of these spectra, and of many 
others like them, indicates that there is no detailed 
correlation in the frequency domain. (The quasi- 
periodic variation superposed on the spectra is 
caused by the finite length of the data samples 
and should not be misconstrued as an indication 
of correlation between the spectra [Blackman and 
Tukey, 1959].) No correlation is to be expected 
for scintillations which result from a diffraction 
process because the distance, measured in wave- 
lengths, between the observer and the diffracting 
screen is very different for the two radiofrequencies. 
Refraction, on the other hand, should cause detailed 
correlation. Thus, we conclude that, at 108 Mc/s, 
refraction plays no observable role in the produc- 
tion of the amplitude scintillations we observe, 
in contrast to the observations made at lower fre- 
quencies by Wild and Roberts [1956b.]._ The lack of 
correlation also indicates that there are no correlated 
instrumental errors, such as might be introduced 
by variations in the power-line voltage or by move- 
ments of the antennas. 

Excepting the short-period fluctuations which 
are the effect of finite sample length, the power 
spectra follow a definite pattern. Typically, log 
W decreases linearly with frequency, up to a fre- 
quency of about 0.1 c/s, and becomes nearly con- 
stant thereafter. The negative slope of the decreas- 
ing portion of the spectrum is directly related to the 
“fading period,” ¢t, which was used in a previous 
study of the amplitude scintillations [Lawrence and 
Jespersen, 1959]. (In that study, ¢ was defined as 
twice the average interval between crossings of the 
average power level.) Although the upper example 
in figure 8 does not seem to concur, the negative 
slopes of pairs of spectra derived from simultaneous 
scintillation records made at the two radiofrequen- 
cies, are generally very similar. Such independence 
between fading rate and frequency is to be expected 
for scintillations which arise from a diffraction 
process. 

The power spectra are nearly flat for frequencies 
greater than 0.1 ¢/s, but the shape of this portion of 
the spectrum is severely influenced by the discrete 
sampling of the original data. When a function, 
z(t), is sampled at equal intervals of time, At, then 
the frequencies 2fy+f, 4fv+f, 6fv+f, ete. are not 
distinguishable from the frequency f [Blackman and 
Tukey, 1959]. Here, fy, the Nyquist frequency, is 
defined as fy='4At. Thus, if frequencies higher than 
fw are present in the data they will contribute their 
energies to lower frequencies in the spectrum. 
This phenomenon is_ frequently described as 
“aliasing” or “folding” of the higher frequencies into 
the range of interest of the power spectrum. 

As mentioned previously, the power spectra could 
be computed only for records which displayed mod- 
erate, slow scintillations. On such records, most of 
the variance due to scintillations is composed of 
frequencies less than fy, and those frequency com- 
ponents greater than fy will be largely due te white 
noise. The effect of white noise, when folded back 





into the spectrum, has been examined and is found 
to increase the general level of the spectrum by a 


nearly constant amount. In order to avoid aliasing 
the power spectra, one would ordinarily use appro- 
priate low-pass filters to remove the higher frequen- 
cies before sampling. However, in the case of the 
“‘zero-crossing’”’ method (see following paper by the 
same authors), the necessary filter would have an 
upper frequency cutoff equal to the phase-sweeping 
frequency and would thus remove the signal. 

The power plots have also been analyzed in terms 
of (AA/A)? and (AP/P), the amplitude and power 
indices [Lawrence and Jespersen, 1959]. The rela- 
tion between (AA/A)? and (AP/P) depends upon the 
distribution law obeyed by the amplitude, A, of the 
radio wave. Thirty-one records selected at random 
were used to produce the dots shown in figure 9, 
where it may be noted that the observed values of 
amplitude index (AA/A)? and power index (AP/P) lie 
close to a smooth curve. 

Since the relation between these two quantities 
depends upon the distribution function of the 
amplitude of the radio wave, this result gives us 
some information about the distribution functions 
themselves. Any particular probability density 
function of amplitude should produce a single value of 
(AA/A)? and a corresponding single value of (AP/P). 
The fact that the points on figure 9 lie near a single 
smooth curve implies that the family of probability 
density functions from which our scintillation samples 
have been taken can be described by a single param- 
eter, the position along the line depending upon the 
value of that parameter. 

The solid line in the figure shows the theoretical 
relation which would occur if the observed signal 
were composed of a steady (undiffracted) component 
and a random component consisting of a number 
(three or more) of independent, randomly phased, 
diffracted waves. The probability density function 
from which this curve has been derived was first 
described by Rice [1944 and 1945] and has been 
discussed in terms of its application to ionospheric 
problems by Ratcliffe [1956]. This probability 
density function, which we shall call the “Rice 
distribution”, is really a family of distribution 
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Figure 9. The relationship between amplitude and power 
indices, compared to the theoretical curve from the Rice 


distribution. 
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functions parametric in the quantity 


amplitude of the steady component 





@e= - : . 
rms instantaneous amplitude of the random component 


__v2(rms instantaneous amplitude of the steady component) 





rms instantaneous amplitude of the random component 


Norton et al., [1955] have defined the quantity 


rms instantaneous amplitude of the random component 








which has the advantage of being a ratio of mean 
squares. k may be thought of as the rms noise-to- 
signal ratio of the radio wave. In the present 
application “signal” means the steady component 
unaffected by the irregularities of the ionosphere, 
and “noise”? means the combination of the various 
diffracted waves. 

The Rice distribution degenerates into the familiar 
Rayleigh distribution in the limit as k approaches 
infinity, and approaches a displaced Gaussian dis- 
tribution as k approaches zero. The details of the 
derivation of this curve have been given elsewhere 
[Lawrence and Jespersen, 1959]. 

The apparent agreement between the amplitude 
scintillation records and the Rice distribution means 
that the arbitrary amplitude and power indices 
may be replaced by the physically more meaningful 
parameter k. (This agreement is also justification 
for the method chosen to produce artificial test data, 
as described in the following paper.) Each depth 
index naturally covers a considerable range of values 
of k. Extreme and average values are approximately 
as shown below. 














| k 
Depth index | 

| min max avg 
1 | 0. 04 0. 12 0. 05 
2 . 05 .13 . 07 
3 13 . 30 . 20 
4 . 28 . 35 ol 
5 | AT >3 | 1. 41 








The family of Rice distributions as defined by 
Rice [1944 and 1945] applies to the amplitude of the 
radio wave and is thus quite convenient for the 
study of communications signals. For radio-star 
observations it is usual to employ a square-law de- 
tector and thus obtain records which are linear in 
power. The Rice distributions, G (P), plotted in 
terms of power, are shown in figure 10. These curves 
accurately match the observed distributions of 
chart deflections. On these curves the long tail 
which exists for large values of k is particularly 
noticeable. The corresponding large deflections of 
the record are indeed prominent, and may easily 
give the erroneous impression that they are caused 





?=2/a*= : : . : 
rms instantaneous amplitude of the steady component 


by some special mechanism unrelated to random 
diffraction. 

Many simultaneous pairs of records at 53 Me/s 
and 108 Me/s have been used to compare the values of 
(AA/A)? at the two frequencies. Booker [1958] has 
pointed out that, in the Fraunhofer region (i.e., in 
the far-zone region, where (AA/A)? is independent of 
the distance between the observer and the diffracting 
screen), (AA/A)? should vary as the square of the radio- 
frequency, while in the Fresnel region (i.e., the near- 
zone region) the variation should be as the fourth 
power of the frequency. Figure 11 seems to indicate 
that, for values of (AA/A)? less than about 0.2 at 53 
Me/s, we are on the border between these two regions. 
For large values (not shown in the fig.) (AA/A)? ap- 
proaches the Rayleigh limit (0.274) independent of 
frequency [Lawrence and Jespersen, 1959]. This is 
expected, since the Rayleigh limit represents the 
vase where all of the signal (undiffracted component) 
has been removed. 

As has been indicated, the noise level may be esti- 
mated from the power spectra. This estimated level 
has been used to compute the contribution of noise to 
the value of (AA/A)? asa function of (AA/A)? at the two 
frequencies. The observed values of (AA/A)? were 
then corrected for noise effects and these are the 
values plotted in figure 11. 

A rough indication of the diurnal variation, mixed 
with an unknown amount of seasonal variation, of 
amplitude scintillations comes directly from the sub- 
jective depth indices. The lower curve in figure 7 
was drawn by smoothing the hourly average values 
of depth index, and converting the result to (AA/A)? 
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FiaureE 10. The Rice distribution function plotted in terms ¢ 
power. 
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FicureE 11. Amplitude scintillations at the two radiofrequencies. 


The expected wavelength dependence for Fraunhofer (A?) and Fresnel (A‘) 
diffraction is shown by the solid lines. 


by means of the table (see appendix 1). These 
curves are certainly too low because of a limitation 
at the high end of the scale of depth index. Index 5 
has been used to describe records in which the scin- 
tillations were “‘so deep that the lobes disappeared’’. 
Since many different degrees of scintillation activity 
are included in this category, the effects of the most 
violent scintillations are not properly represented 
in the lower curves of figure 7. In order to assess 
the importance of this saturation effect, a portion of 
the data has been analyzed by the more elaborate 
procedure of evaluating P, (see following paper). 
The results, when transformed to equivalent val- 
ues of (AA/A)? and then averaged by months, 
are shown by the upper curves in figure 7. These 
curves give a more accurate estimate of the true 
variation of (AA/A)?. The data used in this result 
were a composite of two groups centered around the 
elevation angles 15° and 35°. 

Figure 12 shows the elevation-angle dependence of 
(AA/A)? as computed from P,. Only data from 
August 1958 to February 1959 are involved here. 
For comparison, a curve proportional to sec*x, the 
cube of the secant of the zenith angle, is included on 
the same figure. The observed variation does not 
agree very well with the sec*x curve which is pre- 
dicted by a theory based upon plane-earth geometry 
and spherical irregularities [Booker, 1958]. How- 
ever, most observers agree that the irregularities are 
alined along the earth’s magnetic field and are 
therefore not spherical. Also at the elevations con- 
sidered, the curvature of the earth is not negligible. 
Perhaps the most important reason for the dis- 
crepancy is that (AA/A)? approaches some constant 
value independent of frequency, and hence, for 
large scintillations, any theoretical elevation-angle 
dependence should be modified to include this effect. 

Using the noise-to-signal ratio, k, which indicates 
the fraction of incident energy that suffered diffrac- 
tion in the ionosphere, we can make a rough compari- 
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Ficure 12. The variation, with source elevation above the 
northeastern horizon, of the amplitude index for each radio- 
frequency. 

The broken lines show a theoretical elevation angle dependence proportional 
to the cube of the secant of the zenith angle. 


son of our results with those of other observers. 
Interpolating between our two frequencies, we find 
that, for a given elevation angle, the value of k ob- 
served at Boulder is 2 or 3 times greater than the 
value observed at 80 Me/s at Cambridge, Englana 
[Hewish, 1952]. Since measurements at College, 
Alaska [Little et al., 1958] have shown little or no 
dependence upon zenith angle, direct comparison 
with the Boulder results is less meaningful. The 
value of k at College was consistently 2 or 3 times 
greater than that at Boulder, even near the horizon, 
despite the fact that the College measurements were 
made at 223 Me/s. (The Boulder and College ob- 
servations were made when the relative sunspot 
number was about 190; for the Cambridge observa- 
tions it was 85.) 

Angular scintilations. When normalized to ac- 
count for the spacing of the interferometer lobes, the 
phase plots can be used to calculate power spectra 
of angular scintillations over the frequency range 
0.005 to 0.17 c/s. Figure 13 shows examples of such 
spectra. The ordinate value is the common loga- 
rithm of power density, or variance per unit fre- 
quency interval, in units of tenths of a minute of 
arc, squared. Here, as with the corresponding 
spectra of amplitude scintillations, the reader is re- 
minded that spectra can be obtained only for moder- 
ately disturbed records. Records containing intense 
scintillation activity cannot be used to derive spectra. 

Hewish [1951] has shown that for meter wave- 
lengths, the scale of the ionospheric irregularities is 
the same as the scale of the phase irregularities on the 
ground whenever the rms phase variations are less 
than 1 radian. That this is generally true for our 
data is shown by the fact that the shape of the power 
spectra of the angular scintillations are similar for 
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FictrRE 13. The highest and lowest observed angular power 


spectra for each radiofrequency. 


The lower level at each frequency indicates the background “noise” level. Wp 
is the variance per unit frequency interval expressed in tenths of a minute of are, 
squared. @ is indicated, in units of minutes of are squared. 


records taken simultaneously at the two radiofre- 
quencies. However, as in the case of amplitude 
scintillations, the high-frequency end of the power 
spectra is at times masked by the presence of noise 
which is folded back into the spectrum. Figure 13 
shows the highest and lowest observed levels of the 
power spectra for angular scintillations at the two 
frequencies. As for amplitude scintillations, the 
lower two curves are indicative of the noise level. 
The corresponding values of the variance of the 
angular scintillations, 6¢, are given for each of the 
spectra. 

The shape of the spectrum varies sufficiently from 
day to day to obscure any variation in spectral shape 
which may be related to the elevation angle of the 
source. Inspection of simultaneous power spectra of 
angular scintillations at 53 and 108 Me/s indicates 
that the details of the spectra do not coincide. The 
general shape may vary from hour to hour, usually 
in the same way on both radiofrequencies. 

Inspection of the phase plots of records taken 
simultaneously at the two radiofrequencies indicates 
(except for a special case to be described) that there 
is no detailed correlation. Here again, where dif- 
fraction is involved, no correlation is expected. 

Forty-one pairs of simultaneous phase plots (repre- 
senting moderate scintillation activity) were used to 
compare 6", the mean square angular scintillation at 
the two frequencies. Figure 14 displays this fre- 
quency dependence*of angular scintillations. 
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The solid line is the theoretical wavelength dependence. 
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FicurReE 15. The variation, with source elevation above the 


northeastern horizon, of the rms angular scintillation for each 
radiofrequency. 


The broken lines show a theoretical elevation-angle dependence proportional 
to the square root of the secant of the zenith angle. 


Theory [Booker, 1958] indicates that there should 
be a fourth-power wavelength dependence for 6, and 
our data fits this quite well, Here there is no tend- 
ency to approach a limiting value which is independ- 
ent of frequency, since 6? depends upon the number 
of ionospheric irregularities which the radio wave 
passes through, each irregularity introducing a cer- 
tain phase change dependent upon the frequency of 
the incident wave. 

Figure 15 shows the elevation-angle dependence of 
6 at the two frequencies and figure 16 gives the 
seasonal (and/or diurnal) dependence of @ over a 
limited portion of the year. Both of these curves 
are based upon og (see following paper by the same 
authors) which was then transformed to the appro- 
priate value of 6 Simple theory indicates that & 
should vary as sec’/?x and this seems to be true in 
figure 15. 

Relationships between amplitude and angular scintil- 
lations. About two hundred scatter diagrams of the 
instantaneous power, P, versus the instantaneous 
phase, a, have been examined. In no instance is 
there an indication of a detailed correlation between 
these two quantities. This is the result to be ex- 
pected if a diffraction process is involved. However, 
a definite correlation does appear for the magnitude 
of power and phase variations averaged over 10 min 
intervals. 
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Ficure 16. The seasonal and/or diurnal variation of rms 
angular scintillation. 


Theory [Booker, 1958] indicates that the statistical 
relation between amplitude and angular scintillations 
depends upon the observer’s distance from the iono- 
spheric disturbing region. In the Fraunhofer region 


2 _— 
(2Ay 1 (EY a) 


whereas in the Fresnel region 





(= . = (2) 


w 3 Sal) , 


L is the ionospheric irregularity size and z is the 
observer’s distance from the disturbing region. These 
results are valid only for rms phase scintillations less 
than 1 radian. From (1) it is evident that the rela- 
tion between (AA/A)? and @ is linear, with slope pro- 
portional to (Z/A)? in the Fraunhofer region. In the 
Fresnel region the relation depends upon z, the ob- 
server’s distance from the disturbing region, and, 
furthermore, a plot of (AA/A)? versus @ should coin- 
cide for two different radiofrequencies. 

Figure 17 shows the relation which was obtained 
by fitting a line through the observed scatter of 
points. Because (AA/A)? increases with zenith 
angle, the horizontal scale corresponds to a non- 
linear scale of increasing z. The curves diverge, 
indicating a dependence on A, but they are also 
nonlinear, indicating a dependence on z. We must 
conclude that our observations were not completely 
within either the Fraunhofer or the Fresnel regions, 
but in the transition zone. 
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Fiaure 17. The observed relationship between mean square 
angular scintillation and amplitude index. 


Examination of the spectra of the power and phase 
scintillations from a given record indicates that there 
is no detailed correlation between the two, but that 
they have the same general shape, viz, steep at the 
low-frequency end, and leveling off toward higher 
frequencies. This is in agreement with other 
observations that the rate of fading of the amplitude 
and phase are similar for a given record [Ryle and 
Hewish, 1950]. 

Slow, irregular angular fluctuations. In addition 
to the angular scintillations described above, slow 
irregular phase variations having autocorrelation 
times of some minutes have been observed at both 
frequencies. Figure 18 shows two typical examples 
of these irregularities as observed simultaneously at 
the two radiofrequencies. If one were to remove 
the slow variation from them they would look like 
ordinary scintillation records in which no irregu- 
larities are present. Figure 19 shows one of the slow 
irregular structures measured for a period of 1 hr. 
This figure is a composite of consecutive records 
of the type shown in figure 18. It should be noted 
that the steady angular variation of the radio star, 
due to its diurnal motion, has been removed, as it 
is from all phase plots. Since no attempt has been 
made to determine the absolute phase of the inter- 
ferometer system, the zero levels are arbitrary on 
these diagrams. 

The slow angular variations show high correlation 
between records taken simultaneously at the two 
frequencies. Forty-seven pairs of records, each ten 
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Figure 18. Typical ecamples of slow irregular angular varia- 
tion observed_simultaneously at the two radiofrequencies. 
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Ficure 19. A sixty minute record exhibiting slow irregular 
angular variation. 


108 Mc/s; 1100 a.m., October 17, 1958. 





minutes long, showing slow irregular variations at 
the two frequencies, have been examined to deter- 
mine. the average total angular excursion for each 
frequency. The ratio of the average excursion at 
53 Me/s to that at 108 Me/s ranges from 2:1 to 6:1 
with the average 4.1:1. A four-to-one ratio would 
be expected if ionospheric refraction were responsible 
for the angular excursions. The close agreement 
between the observed ratio and the theoretical 
ratio expected for refraction, together with the high 
cross-correlation, indicates that these slow angular 
variations are due to ionospheric refraction. Tests 
with a target transmitter were made to eliminate 
the possibility that equipment malfunction might be 
responsible for these slow angular variations. The 
negative results of these tests, together with the 
frequency dependence of the variations, indicate 
that the effects are of ionospheric origin. 

Examination of the amplitude variations of pairs 
of records for which there were slow angular scintilla- 
tions, indicates that there is no correlation between 
the amplitude fluctuations on the two frequencies, 
nor is there any long-period variation in amplitude. 
In fact, the amplitude scintillations seem to be un- 
affected by the presence of the slow angular 
scintillations. 

This lack of correlation suggests that the iono- 
spheric irregularities responsible for the refraction 
are acting as prisms, or as weak lenses whose focal 
lengths are long compared with the observer’s 
distance from them. For the case of a prism we can 
use the observed total angular excursion to estimate 
the amount by which the electron-content of the 
ionosphere deviates from its mean value. Smith 
[1952] has shown that a prism will introduce an 
excursion in apparent direction given by 


we. 
2rmf? dx 





Here dn/dz is the gradient, transverse to the line of 
sight, in electron content per unit column along the 
line of sight; f is the radiofrequency; e is the charge 
of an electron; and m is the mass of an electron. 

A typical angular excursion at 53 Me/s is about 
0.5 deg. If we assume that such a gradient is 
maintained over a portion of the irregularity equal 
to, say, 10 km, then from the above formula, the 
total electron-content variation is a few percent of 
the mean value. 

Figure 20 shows the seasonal and/or diurnal varia- 
tion in the frequency of occurrence of these irregu- 
larities. The ordinate value is the percent of obser- 
vations which exhibit slow irregular variations, 
relative to the total number of cases available for 

xach month. The shape of the curve suggests that 
the variation is more diurnal than se: asonal, since 
the peak is centered about two seasons (fall and 
winter). Furthermore, in the months from May 
to December the data were taken during daylight 
hours. The shape of the curve seems to indicate 
that the irregularities build up rapidly as the iono- 
sphere becomes illuminated, and then gradually 
disappear. 
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Ficure 20. The seasonal and/or diurnal variation in the 
frequency of occurrence of slow irregular angular varia- 
tions. 


The total angular excursion is found to decrease 
as the source rises, while the probability of observing 
these slow irregularities seems to increase. Since 
ordinary amplitude and phase scintillations decrease 
as the source rises, the slow irregularities are most 
likely to be seen on records for which the ordinary 
scintillations are small, (AA/A)? less than 0.06. 
However, they have occurred on several occasions 
when (AA/A)? was as large as 0.20. 

It is reasonable to assume that the elevation- 
angle dependence in the probability of occurrence 
of the slow angular variations results from the fact 
that at low elevation angles amplitude scintillations 
are most severe and are likely to limit the ability of 
the ‘‘zero-crossing’’ method to recover the slow 
phase variations from the data. 

The variation in the total angular excursion with 
elevation angle probably arises from the changing 
path length through the ionosphere, although there 
is some possibility that the structure of the irregu- 
larities varies rapidly with latitude. No attempt 
has yet been made to check these ideas quantita- 
tively. 

At the low elevation angles of 10 to 20 deg, the 
distance which the line of sight moves through the 
ionosphere in a period of 20 min (typical of the length 
of time during which any one irregularity is observed) 
is of the order of 200 km. This distance may be 
assumed to be indicative of the horizontal size of the 
ionospheric irregularities in the appropriate direc- 
tion if ionospheric drifts are neglected. 

Examination of the distributions of ionospheric 
parameters scaled from the Ellsworth ionograms 
reveals no significant difference between times of 
slow angular scintillations and other times. A 
possible exception is a marginally significant increase 
in foF2. 





Vitkevich and Kokurin [1957 and 1958}, using a radio- 
astronomy technique, observed large-scale angular 
variations with no corresponding amplitude effects. 
Other authors have apparently observed a similar 
phenomenon [Bramley and Ross, 1951; Whitehead, 


1956; Jones et al., 1957]. The large-scale travel- 
ing disturbances which Munro [1950] has studied are 
possibly related to the disturbances observed here. 
Satellite observations made at the Boulder Labora- 
tories of the National Bureau of Standards [Little 
and Lawrence, 1960] have also disclosed the presence 
of these large-scale, ionospheric irregularities. 


6. Conclusions 


After removal of the diurnal (and/or seasonal) 
trends from the data, amplitude scintillations still 
correlate with spread F. The variation of correla- 
tion with elevation angle of the source is small, 
indicating that spread F is in patches whose sizes 
are at least 400 km. No significant correlation is 
found between amplitude scintillations and sporadic 
E. This result is particularly interesting since the 
experiment was designed especially to measure and 
detect possible #-region effects on scintillations. 

The distribution of flux received from a scintil- 
lating radio source closely matches the theoretical 
distribution to be expected from the superposition of 
several randomly phased, diffracted waves upon a 
constant, undiffracted wave. The zenith-angle de- 
pendence of these amplitude scintillations does not 
agree with a theory based upon isotropic ionospheric 
inhomogeneities. At Boulder, the fraction of the 
radio-star energy diffracted by the ionosphere is 
about twice as great as at Cambridge, England, but 
less than half as great as at College, Alaska. How- 
ever, the Boulder and College observations were 
made when the relative sunspot number was about 
190; for the Cambridge observations it was 85. 

The root-mean-square value of angular scintil- 
lations is found to be proportional to the square of the 
wavelength, in accord with a theory of diffraction by 
ionospheric irregularities. 

Comparison of amplitude scintillations with angu- 
lar scintillations indicates that, for elevation angles 
of 15° to 50°, the region of the ionosphere responsible 
for scintillations lies near the border between the 
Fresnel and Fraunhofer diffraction regions for both 
frequencies. 

Slow irregular angular variations having autocor- 
relation times of 5 to 20 min are commonly observed 
at both frequencies. These variations occur chiefly 
by day and are attributed to large (of the order of 
200 km) lens-like ionospheric irregularities whose 
electron content per unit vertical column varies by a 
few percent from the mean value. The resulting 
angular displacement of a radio source is frequently 
0.5 deg at 53 Me/s, and varies as the square of the 


wavelength. 
7. Appendix 1 


Amplitude scintillation index. A visual estimate of 
the depth of scintillations can be made from a pen- 
recorder trace of the interferometer output. Figure 
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21 shows 12-min sections of two typical scintillation 
records. The modulated sine wave is formed as the 
interferometer lobes sweep past the discrete source 
(Cygnus-A) at the rate of 10/min. The amplitude 
of the sinusoidal trace is proportional to the flux 
received from the discrete source; the fluctuations in 
the envelope are ionospheric scintillations. The 
depth of these amplitude scintillations was classified 
by eye into 5 categories as follows: 








Depth 
index Appearance of record 
1 No scintillations. 
2 Shallow scintillations. 
3 Medium scintillations. 
4 Deep scintillations. 
5 Scintillations so deep that 


the lobes disappear. 





For example, the records in figures 21a and 21b have 
been assigned depth indices 3 and 5 respectively. 








The estimated depth index for each 5-min period was 
inserted upon the appropriate IBM card. 

Many of the scintillation records were supple- 
mented by simultaneous digital recordings on 
punched paper tape. A random sample of 108 Me/s 
digital records was selected for each depth index and 
subjected to detailed analysis on an electronic com- 
puter. This analysis yielded both the power index 


AP _|P—P 


a AF 





and the amplitude index 


CY-Ce) 


where P is the power received from the discrete 
source, A=¥yVP is the amplitude of this radio signal, 
and a bar denotes an average over the record. The 
results are summarized in the following table: 














AP/P (AA/A)? 
Depth | = a 
index 
min max avg min max avg 
| | | | 
1 0.04 | 0. 13 0. 09 0.0008 | 0. 0071 0. 0025 
ae 06 | .14 | . 10 0016 | . 0083 | 0045 
3 ae | eo | . Ba . 0092 . 0448 | 0238 
4 .o0 | . 39 a . 0377 . 0575 0512 
5 . 44 48 iy . 0838 « 3237 . 1689 








FicurE 21. Twelve-minute samples of typical 108 Mc/s phase-sweeping records. 


The modulated sine wave is formed as the interferometer lobes sweep past the source. 


source; the fluctuations are ionospheric scintillations. 


The.amplitude of the record is proportions! t> che flux received from the 
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8. Appendix 2 


Correlation of amplitude indices with vonograms. 
The IBM cards containing the simultaneous iono- 
spheric (Ellsworth) and scintillation data are referred 
to as “statistical index cards”. Tests for correlation 
between any two parameters, z and y, on the cards 
may be performed easily on IBM statistical machines. 
The statistical-sorter output consists of a bivariate 
distribution table, each variable, z and y, being 
allowed a scale of one through nine, zero being re- 
served for cases where there are no data. These 
ten-by-ten distribution arrays are then punched onto 
IBM cards for processing through a computer pro- 
gram which determines the correlation coefficient 
between the two variables. The normalized corre- 
lation coefficient is defined in the usual manner as 


7-77 
FY 
O70y 


(1) 


where the variance o2=2?—7Z*. For independent 
data which are normally distributed, the standard 
deviation of the correlation coefficient is 


ae a (2) 


op" air? 
Tr VN 


where JN is the number of samples. If the data are 
not independent, but possess a significant serial 
correlation over an interval n samples in length, then 
N in (2) should be replaced by N/n. A visual inspec- 
tion of the data indicates n to be approximately 6 
for data taken at 5 min intervals, but nearly 1 for 
data taken at half-hour intervals. For the correla- 
tion coefficients reported here, the standard devia- 
tions have been calculated using these values of n. 

The distribution of the data is approximately 
normal for most of the ionospheric parameters but 
for scintillations it is strongly skewed towards the 
high depth indices. Figure 22 gives the distribution 
of scintillation index for each radiofrequency, both 
before and after revision of the data. The methods 
of, and reasons for, data revision will be discussed 
in the following paragraphs. In order to assess the 
error introduced by using (2) with skewed data, a 
more sophisticated analysis [Cramer, 1954] was 
undertaken. This calculation, involving higher order 
moments of the actual distribution, was made on a 
typical correlation coefficient taken from our data. 
The difference in o, amounted to less than 0.001; 
hence (2) has been used in all subsequent calculations 
of o,. 

In cases where |r| >0.5, or where N is small, there 
is some advantage to use of the Fisher Z, function 
[Walker and Lev, 1953] and its corresponding stand- 
ard deviation ¢,,. Z, is defined as 


= | l-+r ; 
Le =< In i) (3) 
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Figure 22. Frequency distribution of observed scintillation 
index. 


and o,, is obtained from the formula 


| 
va, V¥N—3 (4) 
In the present case where correlation is small and 
N is large, Z, and o,, are practically equal to r and 
o, respectively, and they can be used interchangeably. 

The geometry of the experiment indicates that 
any correlation between ionospheric and scintillation 
parameters should vary with local sidereal time 
(LST) as the ray path from Cygnus-A moves through 
the different effective ionospheric levels over Ells- 
worth. Also, scintillation activity is found to 
decrease with increasing elevation angle (increasing 
LST) of the source. Consequently, the statistical 
index cards were divided into groups according to 
the elevation angle of Cygnus-A. Figure 3 gives 
these groupings as well as a complete distribution 
table of the cards as a function of the elevation angle 
of Cygnus-A. Figure 3 also gives the relationship 
between the elevation angle of Cygnus-A, local side- 
real time (LST), and effective ionospheric heights. 
The average elevation angle of each major (10°) 
group has been computed and is used for plotting 
the correlation coefficients in section 5. 

Early correlation studies, reported previously 
[Lawrence and Jespersen, 1959], consisted primarily 
of bivariate correlations as a function of LST, and 
used only six months’ data. With the cessation of 
recording, the thirteen months (February 1958 
through February 1959) of available data were 
placed on statistical index cards. These data were 
then tested for correlations, but with few consistent 
results. By contrast, the previously reported 
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results [Lawrence and Jespersen, 1959] of correlation 
of a portion of the data had appeared internally 
consistent. Prominent diurnal and/or seasonal 
effects in the full year’s data are presumably 
responsible for this difference. 

Division of the data into two classes, sunlit and 
dark ionosphere, proved helpful and indicated that 
complete removal of diurnal and seasonal trends 
would be worthwhile. To determine the trends, the 
hourly average value of each variable was plotted, 
except that in the case of spread F' and sporadic EL 
the frequency of occurrence was used, as shown in 
figure 6. The hourly average values of scintillation 
depth index, converted to (AA/A)?, appear in figure 7. 
These graphs are referred to as diurnal curves, 
although each contains an unknown amount of 
seasonal variation. It is impossible to distinguish 
with certainty between diurnal and seasonal effects, 
since only one radio source was used. 

The IBM 650 computer was used to subtract the 
appropriate hourly average from each numerical 
value on the statistical index cards, and to normalize 
the remainder to a scale of one through nine. The 
resulting cards are referred to as “revised statistical 
index cards’, and have been used in all the correla- 
tion results reported here. Figure 22 shows the 
effect of this revision upon the distribution of the 
scintillation indices. This redistribution cannot, of 
course, remove the effect of skewness in the original 
depth indices, as this was a result of limited equip- 
ment sensitivity. 

The data in final form consist of several thousand 
revised statistical index cards. The extreme excur- 
sions of the inospheric and scintillation indices from 
their average value for any particular local time, 
now appear as extreme values on a scale of one 
through nine. Most of the data are clustered about 
the median of the scale and so they contribute little 
to the correlation coefficient. Thus the revised 
values are in a form which may be tested for detailed 
correlation between individual events, regardless of 
trends in the averages of the original variables. 


This experiment could not have been performed 
without the cooperation of several people. Among 
them are E. R. Schiffmacher and H. A. Erickson 
who were responsible for the maintenance and 
operation of the scintillation receivers, R. F. Carle 
who established and operated the ionosonde at 
Ellsworth, and Mrs. C. J. Enfield who assisted with 
the analysis of the data. This project was supported 
by the Ballistic Missiles Division of the United 
States Air Force under D. O. AF 04(647)-134. 
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Digital Methods for the 


Extraction of Phase and 


Amplitude Information From a Modulated Signal ' 


R. S. Lawrence, J. L. Jespersen, and R. C. Lamb 


(November 3, 1961) 


A description is given of three digital methods which have been used to recover ampli- 
tude and phase information from a modulated sinusoidal signal sampled at equal intervals 


of not more than one-sixth of a period. 


The first method, the ‘‘zero-crossing’”’ method, is 


economical of computer time and, for modulation which is not too deep and does not contain 
frequencies near the carrier frequency, accurately recovers the phase and amplitude modu- 
lation. The second method, the “‘filter’’ method, is more laborious but it gives better accuracy 


and will operate with deeper and more rapid modulation. 


The third method, a statistical 


approach, will work with severly overmodulated signals, but it yields only a statistical sum- 


mary of the modulation. 


The methods were designed specifically for analysis of radio-star 


scintillation records but they may be applied to many other modulated signals. 


]. Introduction 


Three digital methods are described which have 
been used to recover amplitude and phase informa- 
tion from a modulated sinusoidal signal sampled at 
equal intervals of not more than one-sixth of a 
period. Although these methods were designed 
specifically for the purpose of analyzing the records 
described in the preceding paper, ‘Amplitude and 
angular scintillations of the radio source Cygnus-A 
observed at Boulder, Colorado”, they may be ap- 
plied to many other modulated signals. In general, 
the desired modulation is imposed on a sinusoidal 
carrier of constant frequency and amplitude. The 
sample interval and the carrier frequency must be 
chosen carefully with regard to the expected range 
of modulating frequencies, in order to avoid biasing 
the results either by omission of or folding at fre- 
quencies of interest. 

All three of the methods described herein were 
designed for use with the IBM-—650 digital computer 
and its peripheral equipment which included an 
automat ¢ graph plotter. The input data were 
digital samples of the modulated carrier, sampled 
to an accuracy of 114 percent of full scale (six bits) 
and punched on paper tape. The first method of 
analysis, the ‘‘zero-crossing method”’, uses the prop- 
erties of the signal only at the times when it passes 
through its average value. Phase modulation is 
derived from variations in the intervals between 
these times; amplitude modulation is derived from 
the slope of the signal at these times. The second 
method, the ‘‘filter method’’, passes the signal 
through appropriate digital filters which recover the 
phase and amplitude modulation. The last method 
is a statistical one, used to analyze records which, 
because of overmodulation in amplitude and/or 
phase, are not amenable to treatment by either the 





' Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 
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zero-crossing or filter methods. This statistical 
method cannot give amplitude and phase as a func- 
tion of time, but produces results which describe 
some statistical properties of these quantities. 


2. Zero-Crossing Method 


In this procedure the computer must deal with 
the modulated carrier at the times when it crosses 
its average value. The first step is to replace the 
original samples by deviations from their average 
value. Since drifts in the digital recording system 
might cause the carrier to undergo slow changes in 
average (d-c) level, changes unrelated to the modu- 
lation, it is usually better to use a moving average 
several periods in length than to use an average over 
the entire record. If a moving average is used, the 
process of replacing the original samples by devia- 
tions from the average is, in fact, a digital high-pass 
filter which removes the d-c level and its slow varia- 
tions, but does not appreciably affect the carrier or 
its modulation. For the radio-star observations, a 
moving average 30 samples in length (approximately 
5 periods) proved satisfactory. 

After the filtering operation, the computer can 
find the points of interest by searching for adjacent 
samples with opposite signs. Each such ‘crossing 
pair’ of samples yields, by linear interpolation, a 
good estimate of the exact time of zero crossing. 
In order to extract phase information, first differ- 
ences are taken of the list of times of zero crossings. 
Each of these intervals is then replaced by its devia- 
tion from the mean interval. The cumulative sums 
of these deviations provide a list of numbers called 
the phase table. The phase table is simply a tabu- 
lation of instantaneous phase variations at half- 
period intervals. 

From the crossing pairs the computer derives 
another set of numbers called the amplitude table. 
This tabulation contains values proportional to the 
amplitude of the modulating signal. It is the sum 











of the absolute values of the two samples involved 
in each crossing pair, and measures the instantaneous 
amplitude at half-period intervals. It will be 
demonstrated later that this simple method, which 
might be expected to be inaccurate, is, in fact, an 
effective means for extracting the envelope of the 
sinusoidal trace. Figure 1 gives examples of actual 
amplitude and phase plots, made by an automatic 
plotter from the amplitude and phase tables. Each 
plot represents 600 points (10 min) of data and has 
been expanded to occupy the full scale. The full- 
scale range of each phase plot is shown in electrical 
degrees. The reliability and sensitivity, together 
with a brief description of the drawbacks, of the 
zero-crossing method will be discussed in the follow- 
ing sections. 


3. Filter Method 


Amplitude and phase tables may also be obtained 
by passing the data through a system of digital 
filters. This analysis procedure is based upon a 
method suggested by N. R. Goodman [1960], of 
Space Technology Laboratories. 

Consider a modulated sine wave which is sampled 
once every t, seconds to produce a time sequence of 
n data points, which we designate z;. We are inter- 
ested in the amplitude (i.e., envelope), and phase 
of this record. A possible representation of the 
amplitude A(¢) and phase a(t) is 


A(t) exp [ia@)]=X(®)+iY@ 


where X(t) and Y(¢) are the real and imaginary parts 
of the complex representation of the power and 
phase. We may let the x, represent a set of sampled 
points from the X(t) process. If we can somehow 





construct from the z; a set of y; which will be a good 


representation of the Y(t) process, then we may 
compute A(t) and a(t) from the relations 


A®=L?' 0+ PO)" 
a(t)=sin-" [Y()/A®]=cos“' [X()/AQ]. 


One possible representation of the y; in terms of 
the zx; Is 


Y= DO, Ver, (1) 
3 Om 


where the a, have suitable values. There are several 
ways in which one might evaluate the a,. Goodman 
shows that for the present purpose the a, should be 
chosen in such a way that, over the frequency range 
of interest, the spectral density matrix of the result- 
ing (x;, y;) process is a good approximation to the 
spectral density matrix of the actual [X(¢), Y(6] 
process. Following this line of thought, a set of a, 
were chosen to satisfy that criterion. For the case 
where the sampling rate was about 6 samples per 
period, the following a, proved satisfactory. 


ap = 0.16 
a = 15 
a2 =—.13 
a3 = — .20 
As = — .06 
a; = .03 
ag = .00 
az = — .02 
ag = .03 
Ag = 05 
Ajo = .02 
a= — 01 
a2= — 03 
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Each record is 600 samples (10 min) in length. 
full scale. 


SEPTEMBER 21,1958,CYGNUS—A,ELEVATION ANGLE =38°, DEPTH INDEX=2 108 Mc/s 


Examples of amplitude and phase plots obtained from the computer. 


The bottom edge of the amplitude plot represents zero signal, and the plot has been normalized to occupy the 
The scale of the phase plot, a, given on the figure, is in electrical degrees. 
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These a, were then used in conjunction with (1) to 
evaluate the y;. 

Since the y; are a moving linear combination of the 
z;, they are in fact the result of passing the x; through 
a linear digital filter. This filter (y filter) has a 
narrow pass band, approximately unit gain, and 
7/2 phase shift in its pass band. The phase shift, 
uniform over the pass band, is accomplished by mak- 
ing the filter anti-symmetric; i.e., a,=—a_,, and 
a=0. Noise outside the pass band of the filter 
may be removed from the x; by means of a similar 
filter (x filter) with zero phase shift. This will be 
a symmetric filter, a,=a_,. 

The amplitude and phase tables produced by either 
the “‘zero-crossing”’ or “‘filter’’ methods may be used 
for such subsequent purposes as spectral analysis 
and cross correlation between amplitude and phase. 

In order to test the ability of the ‘‘zero-crossing”’ 
and ‘filter’? methods to reproduce accurately the 
amplitude and phase fluctuations inherent in the 
data, an artificial record with known amplitude and 
phase variations was analyzed using each method. 
The known amplitude and phase variations are 
shown by the heavy curves in figures 2 and 3. The 
light curves labeled ‘‘two-filter method” show the 
results of using both the z and y filters. To assess 
the possibility of noise contamination, the same 
artificial record was analyzed using only the y filter. 
The result of this analysis is shown in the figure 
labeled ‘‘one-filter method.” 

Although both methods reproduce the known 
amplitude and phase variations quite well, the one- 
filter method includes a high-frequency component 
which is not evident in the two-filter method, indi- 
cating the existence of noise outside the pass band 
of the filter. 

The heavy curves (input) are tilted slightly with 
respect to the light curves (output). This is because 
the input data is passed through a high-pass filter 
which removes trends. 

The figures labeled “‘zero-crossing method’? show 
the corresponding results produced by the zero- 
crossing method when it is applied to the artificial 
data. 

The principal difference between the filter method 
and the zero-crossing method is that the former 
utilizes consecutive data points whereas the latter 
uses only pairs of consecutive points which have 
opposite signs (crossing pairs). For this reason the 
filter method has a twofold advantage: (1) mod 
2r ambiguities are less likely to occur between 
consecutive points than between crossing pairs, and 
(2) higher frequency components are preserved by 
the filter method since its uses data points which 
are spaced more closely in time. The two-filter 
method involves about four times as much computing 
time as does the zero-crossing method, but the 
fidelity of the more lengthy method is better. 


4. Statistical Method 


Both of the foregoing methods have been applied 
to artificial data containing various amounts of 
amplitude and phase modulation. Although the 





filter method is more capable of producing accurate 
amplitude and phase plots than is the zero-crossing 
method, both methods fail for some records. These 
highly modulated records contain significant phase 
and amplitude modulation with frequencies com- 
parable to the carrier frequency. 

In order to analyze such records, a statistical ap- 
proach was taken. This procedure does not give 
detailed amplitude and phase plots but instead gives 
two quantities, P, and o,, related respectively to 
the magnitudes of the amplitude and phase fluctua- 
tions. 

P,, is evaluated by first finding the points on the 
actual modulated sine wave, P(t), which lie above 
the peaks of a sine wave having the same average 
amplitude. The average value of these points, 
divided by the average amplitude, is defined as P,. 
The sampled points which enter into the evaluation 
of P, all result from amplitude modulation, since 
phase modulation can never increase the instantane- 
ous amplitude of the record above the peaks of the 
equivalent sine wave. Thus, P, is a statistical 
measure of the amplitude modulation on the record 
and is unaffected by the phase modulation. The 
standard deviation, og, of the intervals between zero 
crossings, is a measure of the phase modulation alone 
since amplitude modulation cannot affect the time 
of zero crossings. 

It is desirable to relate the quantities P,, and o, to 
quantities which are more commonly associated with 
amplitude and phase effects. For example, the 
quantities 


and oq were of interest in the preceding paper. 
Here A is the amplitude of the signal, A the mean 
amplitude, and o, is the standard deviation of the 
phase fluctuations, a. 

One possible way to relate P, and og to (AA/A)? 
and oq respectively, is to compute (AA/A)? and og 
from the variations in amplitude and phase for a 
particular record as given by the ‘“‘zero-crossing”’ or 
“filter”? method, and then to find P, and oe, for this 
same record. If enough records are treated in this 
manner, it is possible to obtain regression curves of 
(AA/A)? versus P, and of o. versus og. This pro- 
cedure not only gives the desired relations but yields 
confidence limits for P,, and o,as measures of (AA/A)? 
and oq respectively. However, it cannot be applied 
to the highly modulated records for which the sta- 
tistical quantities P, and oq were originally intro- 
duced. A better method is to compute P, and og 
from artificial records which have a predetermined 
value of (AA/A)? and og, for the particular example 
under consideration. 

The artificial records were constructed in the fol- 
lowing way. A pair of random normal deviates (m,, 
n;) are used to compute the quantity 


P,=(a+ m)?+n3 


P, may be thought of as the 





where a@ is a constant. 
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square of the envelope of a signal at time ¢; which 
consists of the sum of a sine wave with amplitude a 
and noise with mean square amplitude of unity. 
The P; are then passed through a low-pass filter 
which auto-correlates them over a length of 30 points. 
The artificial data points, d;, are then computed 
from the following expression: 


d,=P;, sin (wot ;+a:), 





where a;,, the phase fluctuation, is normally dis- 
tributed with mean zero, standard deviation o, and 
auto-correlation length of 30; and w, is some constant 
angular frequency, approximately 7/3 per sample 


interval. The d,; are then normalized so that the 
artificial records will have the same range in value 
of the d; for any a. 

The probability distribution function of the (P;)'/? 
constructed in this manner has been derived by S. O. 
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Figure 2. An illustration, using artificial data, of the effectiveness of the zero-crossing and filter methods in reproducing 


amplitude fluctuations inherent in the original record. 
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Rice [1944 and 1945]. (The reason for constructing 
the artificial data in this particular way was, of 
course, dictated by knowledge of the sort of signal to 
be analyzed. Signals with different modulation 
characteristics would require different techniques for 
producing the artificial data.) 








The method involving artificial data has the ad- 
vantage that it cannot be biased by mod 27 am- 
biguities since o, is a predetermined quantity. Thus 
the artificial records may be used to extend empiri- 
cally the relation between o, and o, to regions where 
o is much greater than 1 radian. P,, measured on 
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Figure 4. The relationship, derived from artificial data, 
between the statistical quantity P, and the actual value of 
(AA/A)?. 





the artificial records, can be related to (AA/A)? 
through the parameter a which is a known function 
of (AA/A)? [Lawrence and Jespersen, 1959]. 

For values of o as great as 3 radians, the linear 
relationship 


oa=0.55 oa 


appears to hold. Figure 4 shows the observed 
relationship between (AA/A)? and P,,. 





The development of these digital methods was 
sponsored by the Ballistic Missile Division of the 
United States Air Force, under D. O. AF 04(647)-134. 
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Comparison Between Mode Theory and Ray Theory 
of VLF Propagation’ 


H. Volland 


(November 10, 1960) 


It is shown that the field strength according to mode theory and ray theory in the 
VLF band are derivable from the same expression of the original vector potential, and the 
result of one theory is the analytic continuation of the other one in another range of con- 


vergence. 


In fact, both ranges of convergence overlap. Estimates of these 


ranges are 


made and an example shows that within this overlapping region (between distances of 300 


and 2,000 km) both theories give the same result. 


Using this fact calculations of frequency 


spectra are possible which in the case of a white noise dipole show some similar features 
to measured frequency spectra of lightning discharges. 


1. Introduction 


Electromagnetic waves at frequencies below 100 
ke/s are reflected during daytime at the bottom of 
the ionosphere at about 70 to 80 km height. At 
this height the collision number is still greater than 
the angular frequency w, and the depth of penetra- 
tion of these waves into the ionosphere is limited. 
For distances greater than about 500 km and fre- 
quencies below about 50 ke/s the actual ionosphere 
during the day can be replaced to a good order 
of approximation by a homogeneous isotropic plasma 


; aie N.e? 
with conductivity i and a sharp lower 
v 
boundary. (N, and » are mean values of electron 


density and collision number respectively, and «; is 
of the order 107§ Q-! em™.) 

For the ‘‘medium” distances between about 500 
km and the limiting geometric skip distance of the 
first hop wave (1,700 km at 70 km reflection height) 
higher order terms both in mode theory and ray 
theory may be important. It is the purpose of this 
paper to compare mode theory and ray theory in 
this distance range. 


2. Setting up the Equations 


We wish to calculate the vector potential of the 
commonly assumed vertical dipole within the system: 
flat earth and isotropic homogeneous ionosphere 
with sharp lower boundary. Following the method 
of H. Weyl [1919] it can be shown that the vector 
potential of this vertical dipole is 


T/2+i@ 
I(p,$,2)=— ikp * ett $005 (6-4") SQ’ dp! 





! Contribution from Heinrich-Hertz Institut, Berlin-Charlottenburg, Germany 
(Western Zone), 
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+R, @ ~ FRC (2anh+29—- 24 Pe 
+Rie-acoment \') osesh (1) 


with 
S=sin 6’ 


C=cos 6’ 


n? C—./n?—S? 
e 





Ve 
ae. ba ee 
e neC+../n2?—S? 
e We 
JG: 
e Epw 
2) ionosphere 
a,conductivity of < — 
f earth 


(p,¢,2) cylindrical coordinates 
Z : receiver 
height of the antenna of ; 
Zo transmitter 
h virtual height of the ionosphere 
> 9 


p=py e/*' moment of the dipole 


Qr 
k=~ wave number 


The first two terms of the integral (1) represent the 
vector potential of a vertical ‘dipole above a con- 
ductive earth without taking into account the 


ionosphere, a well-known problem that has been 
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solved by A. Sommerfeld [1948] and which gives 


5 ek VeP+(2—20)? 


2 
P p?+(2—2)? 


W (a) (2) 


W(a)=|W\e-* is the Sommerfeld absorption factor 


a=—jkp[{1—n,(1+n?)-'”] is the numerical distance. 


Again, following Weyl [1919], the terms of the sum 
within the integral are approximately 
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yn r 

yn 

This (saddle-point) approximation is valid when 
ryn>d and C,, is not too small compared with unity. 
|n;| should also be somewhat greater than unity. 

Setting z=z,=0 and neglecting higher order 
terms of 7,, the vertical component of the electric 
field strength becomes 


E,=2E, [ Wie) +2 >» sin’ 6, 


n=1 


R2(0,)R2-*(6,) CEE were | (4) 


— jkp 
with E,= = *? ~ aed (free space value of field strength 


of the pil ‘al dipole at the distance p) 


sin 6,—=P 
rn 


r2=p*+(2nh)? 


Equation (4) was derived explicitly by Wait and 





with 


RAA+R je7 286% *0?) 








j= 1— R,R,e~ 2° 
a R,(Ree~280C + ¢ — 28 C0-20)) 
ag 1—R,R,e~ 2" 


Following Budden [1951] and Alpert [1955] (who 
set R,=1) equation (5) may be approximated by a 
sum of modes. For z=z,=0, it has the form 


a? RS 
n=0 E 
[1+abe Sc dC 


Here the values C,, are determined from the eigen- 
value equation 


27 De 


TI(p,9,0)=—F 2) (6) 








RAC eee : a=0, 1, 2, . . (7) 


The vertical component of the electric field is then 
(using the approximation of the Hankel function 
H§?(x) valid for large values of x) 


E,=2Ey Ve" e7sr/4 > 5, S3/2¢I-S,) (A<<p) (8) 
n=0 


ee 1 _ Sf i(n=1) 
ane 7 dR, ~1 0(n=0) 


(+a ac , C=C 


with 











for f>5 ke/s and = is of the order of unity. 
€ 


Equation (8) in the language of wave optics con- 
sists of the possible propagation modes. Actually it is 
the approximate continuation of the ray equation into 
another range of convergence. More general forms 
of eqs (6), (7), and (8) are given by Wait [1957a] 
and Wait and 'Spies [1960] for an imperfectly con- 
ducting and curved earth. 


3. Estimate of the Ranges of Convergence 


We wish to confine the sums in eqs (4) and (8) to 
the first five terms. This is reasonable provided 
that the sixth term is smaller than 1 percent of the 


Murphy [1957], using a geometrical-optical approach. = 
It can be interpreted in the language of ray-optics | first one. Thus, for =1, we get from eq (4) 
as consisting of the ground wave, a wave reflected ™ 
from the ionosphere, and further waves reflected - 3 9 |P5(a.\| -2 
between the ground and the ionosphere any number | __ Pain Ry < 10 
of times. “ 
Carrying out the integration over ¢’ and forming _ PP van 
the sum, (1) becomes - — (9) 
jkp . From eq (8) we find 

M(0,6,2)—" | HP (kSp) { e~ #C'2—20l + fe BCU+20) 

- | S \ 3/2 a | x 

t+ fret #O—*0) } Sda’ O<c <h (5) (st) eFkel'S) */<10 2 (10) 
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S2+C2=1. (11) 
Equation (11) is the first approximation of the 
eigenvalue-equation (7) according to J. R. Wait 
{1957a]. 

Because we have 


Se 
S; 


~1 





eq (10) becomes 


bo] 


e~%< 107? with y= = 


| 


(72 
6 


or 


o> 0A (12) 


The range of small convergence of eq (8) simul- 
taneously is the range of good convergence of eq (4) 
and vice versa. In fact, both ranges overlap. 


4. Comparison Between Ray Theory and 
Mode Theory 


Equations (4) and (8) give only approximate 
results of equation (1). To get an idea of their 
accuracy examples have been calculated with the 
following parameters: 


C7 ~ 
R,.=1; —=1; f=15 ke/s. 
€gW 


For the determination of the eigenvalues C,, (eq (7)) 
tables have been used which are reported elsewhere 





[Volland, 1959]. Table 1 contains values of the 


argument 
E, 
[2E 
and the phase of 
E, 
2K, 


depending on the distance p with an ionosphere of 
constant height h=70 km, calculated from eq (4) 
and eq (8), respectively. Table 2 shows the same 
values depending on the height A and with constant 
distance p=1,000 km. Figure 1 is a graphical form 
of tables 1 and 2. 

From tables 1 and 2 it can be seen that between 
300 and 2,000 km distance and between 35 and 100 
km height, respectively, both equations give the 
same result to within about 5 percent in argument. 
Figure 2 illustrates the influence of the higher order 
terms for the case h=70 km. In figure 2a, corre- 
sponding to eq (4), the distance O—A represents the 
constant ground wave. The vector drawn from A 
to one of the circles represents the first hop wave, 
the thin lines being equivalent to the vectors of 
second, third, and fourth hop waves, respectively. 
The vector drawn from the point O to the solid line 
which connects the end points of the thin lines gives 
the values of argument and phase for the different 
distances tabulated in table 1, second and fourth 
row. 

Figure 2b shows the first five terms of the sum 


3/2 
n=1 S; 


in the same sequence as in figure 2a. The distance 
O-A is equal to the constant first term 1, ete. The 
vector drawn from the point O to the solid line must 
be multiplied by the first term of eq (8) 


a e7 48/4 §3/2 gike(1—S)) 


to get the values in the first and third row of table 1. 


TABLE 1. Argument a and phase p of E,/2E, as a function of the distance p with constant h=70 km according to equation (4) (symbol 
“r”) and equation (8) (symbol “‘w’’). 
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TABLE 2. Argument a and phase p of E./2Eo as a function of the height h with constant p=1,000 km according to equation (4) 
(symbol ‘‘r’’) and equation (8) (symbol “‘w’’). 
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Fiaurefl. Argument a (solid line) and phase p (dashed line) of 
E,/2E, depending on (a) distance p and (b) height h. 


From figure 2 can be seen the way in which the 
higher order terms of eq (4) increase and the higher 
order terms of eq (8) decrease with increasing 
distance. The minima of argument in figure la, 
beginning with the chief minimum at about 500 km, 
result from the interference between (n—1)™ and n™ 
hop waves and from the interference between the 
first and second order mode, respectively, as can be 
seen from figure 2. This is in accordance with an 
idea from H. Poeverlein [1959]. 


5. Discussion 


The strong equivalence between ray- and mode- 
theory as shown above permits the calculation of 
field-strength values over a wide range of frequency, 
height and distance using eq (4) and eq (8) alter- 
natively. In the case of medium distances it is also 
possible to replace the mode picture by the ray 
picture when taking into account higher order 
corrections of the earth’s curvature and an inhomo- 
geneous and anisotropic ionosphere. In such cases 
mode theory becomes very difficult, and even the 
calculation of the first order mode is extremely 
laborious [Wait and Spies, 1960]. In ray theory it 
is sufficient to replace the Sommerfeld ground wave 
by the ground wave above a curved earth according 
to Van der Pol and Bremmer [1949], to exchange the 
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b 700 


Figure 2. Sum of (a) ray-vectors and (b) mode-vectors depend- 
ing on distance p at constant height h=70 km. 


The ratio of the scales of figure 2a and figure 2bis1:2. (For details see text.) 


reflection coefficient and to correct for the earth’s 
curvature, for the diffraction [Wait and Murphy, 
1957]. Furthermore, the factor 1+F,(@,) must be 
replaced by a more complicated function which is 
valid when @ is near 90°. 

Examples of calculations covering a wide range of 
frequencies involving the use of both eqs (4) and (8) 
are shown in figure 3 and figure 4. Figure 3 gives the 
argument of #,/2H, depending on frequency for the 
constant heights A=70 km (solid line) and A=60 km 
(dashed line). These curves are equivalent to the 
complex transfer ratic of a white noise dipole. 














Bm © 


-— OO =m OY O ™ => 











oO 10 20 30 40 50 
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Ficure’3. Complex transfer ratio of the vertical components of 
the electric field of a white noise dipole for a distance of 1,000 


km. 
Virtual height of the ionosphere h=70 km (solid line) and h=60 km (dashed 
line). 
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Figure 4. Mean frequency spectrum in the proximity of a light- 
ning discharge multiplied by the complex transfer ratio of the 
vertical electric field of a white noise dipole for different 
distances. 
(solid line: h=70 km; dashed line: h=60 km) 
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These values, multiplied by the mean frequency 
spectrum of a lightning discharge in its proximity 
according to Taylor and Jean [1959], give figure 4a. 
They should be equal to the measured spectrum of a 
mean discharge at 1,000 km distance. Figure 4b 
and figure 4c give similar calculations for distances 
of 2,000 and 3,000 km, respectively. The calculated 
values, especially figures 4a and 4c, show some 
similar features to the mean frequency spectra of 
lightning discharges measured by Obayashi et al. 
[1959], particularly the double maximum between 
10 and 20 ke/s and the relative maximum at about 
40 ke/s. Very similar curves have been presented 
by Wait [1957b] who employed the mode sum 
exclusively. 

During a solar flare the reflection level descends. 
This is represented by a 10 km lower level of h=60 
km in figure 4 (dashed line). Measurements of 
field strength during such an event generally show 
a decrease below about 15 kc/s and an increase above 
15 ke/s [Obayashi et al., 1959]. In figure 4 this is 
only true up to about 30 ke/s. The disagreement 
between calculated and measured values above 30 
ke/s may be caused by neglecting second order terms 
of the earth’s curvature, which according to Wait and 
Spies [1960] become more important with increasing 
frequency. 
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Antenna Coupling Error in Direction Finders ' 
Charles W. Harrison, Jr. 
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Antenna coupling or scattering error in an interferometer-angle-measuring system 
consisting of two identical base-loaded quarter-wave vertical antennas over a perfectly 


conducting earth is investigated. 


A curve is supplied for the error in angle as a function of 


antenna spacing for an incident vertically polarized electric field arriving in the plane of 
the antennas, when the ratio of antenna length to radius is 75 and the base loads are resistors 


of 50 ohms, 


1. Introduction 


An 
missile trajectory determination—in rudimentary 
form—consists of four identical base-loaded vertical 
antennas, one-quarter wave or less in length, erected 


interferometer-angle-measuring system for 


at the vertices of a square. The phase angles be- 
tween the voltages developed across identical load 
impedances in diagonal antenna pairs are measured. 
These angles establish the direction cosines of the 
missile from which the signal emanates. In addition 
to electronic missile tracking equipment, many fixed 
antenna direction finders utilize this principle. 

There are several sources of error in interferometer- 
angle-measuring systems. It is convenient to class 
these errors under two headings: geometrical and 
electrical. Simmons [1959] has studied geometrical 
errors in a radio interferometer. Under electrical 
errors, Carswell and Flammer [1957] investigated the 
error due to multipath. In the present paper antenna 
coupling error in a two-antenna interferometer is 
investigated theoretically and numerically. The 
model employed consists of two isolated parallel 
nonstaggered antennas of the same length and 
radius, center-loaded by identical impedances. The 
incident electric field is assumed to be polarized 
parallel to the axes of the antennas, but this restric- 
tion in the theory may be removed. The azimuth 
angle of wave arrival that is indicated by the equip- 
ment is computed and compared to the assumed 
angle of wave arrival. The computed error angle 
then applies to a base-loaded two-antenna array 
over a perfectly conducting earth for which the 
loading impedances are halved. Second order 
accuracy is achieved in the theoretical development 
and in the numerical results of this report. In 
another paper Harrison [1960] studied the more 
general case of the four-antenna direction finder and 
obtained results accurate to zeroth order. 

Some insight into the nature of antenna coupling 
error may be gained from the following remarks: 
Each antenna in the array continuously samples the 
scattered field from the other antennas comprising 
the array. The resultant field acting on each antenna 
in the interferometer is the vector sum of the field 
maintained along its surface by the distant telemetry 
transmitter on the missile and the several scattered 


! Contribution from the Sandia Laboratory, Albuqerque, N. Mex. 
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fields maintained by the currents in the other 
antennas. Thus the coordinates of the missile, 
determined by the phase-comparison method, must 
bein error. The error is a function of the polarization 
of the electric field, the coordinates of the signal 
source, the antenna dimensions (length and radius), 
the spacing between elements, the load impedance, 
and the properties of the earth. In this paper the 
scattering error in an interferometer that is far from 
the earth, or that is erected on a very highly conduct- 
ing plane, is determined. The effect of an imper- 
fectly conducting earth is not treated. 


2. Integral Equations for the Currents in an 
Unloaded Two-Element Receiving Array 
and Their Solutions 


Consider two identical non-staggered parallel 
receiving antennas, labeled 1 and 2, oriented as shown 
in figure 1 with respect to a system of rectangular 


z 


ANT ANT | 


Tag 


ye D=2a 








ANT TERMINALS 
(SHORTED) 














FicureE 1. Two-element receiving antenna array (load terminals 
short-circuited) . 





and spherical coordinates with common origin. The 
half-length of each antenna is hf, the radius a, and 
the spacing between them is b. It is assumed that 
h>>a, and Ba<<.1. Here £8 is the radian wave 


number. The terminals of the receiving elements 
b — 
are located at z=0, y= +5) and z=0. The incident 


electric field EF arrives from the distant source at 
the azimuth angle #, measured from the positive 
z-axis, and is directed parallel to the axes of the 
wires. The currents in antennas 1 and 2, when the 
terminals are short-circuited, are designated I 2(2) 
and Jog(z), respectively, and are functions of E, ®, 
and the electrical dimensions of the array Ba, Bh, 
and 8b. From symmetry considerations it is clear 
that Tiz(2)Aloz(z) except when @=nz, n=0, 1, 

~ When the electric field is directed parallel to the 
axes of the wires, the simultaneous integral equations 
for determining the currents J\z(z) and J:,(z) are 
[King, 1956, p. 517, eqs (2a) and (2b)] 


[i tusle)Kaleede'+ [ese VK2,2") de! 


=-iF{o cos B2+Ue’ =). (1) 


% Tox 2’)Kel2,2/)d2’ +] Iie 2’)Ko(2,2")d2! 
y —J > sine 
=~ 5 F{ Cocos: +Ue reel, (2) 


8 V (2-2) +02 


ieee, (3) 


Here 
K,(2,2’ )= 


K,(z,2’) is obtained from (3) by writing 6 for a, 
¢ is the characteristic resistance of space, C, and 
C, are constants, and 


U=—Eje. (4) 

The currents J,z(z) and Jz2(z) may be written 
Tix(z)=Te(z)+Te(2) (5) 
Tze(2)=Ie(z)—Ie(2) (6) 


where Ji(z) and Jg(z) are the symmetrical and 
antisymmetrical components of current in either 
antenna. The assumed directions of these currents 
are indicated in figure 1. 

Substituting (5) and (6) into (1) and (2), the 
following independent integral equations for the 
currents Ts (z) and J£(z) are obtained: 


= {c cos6z+U*} (7) 


{ Py(2")K*(2,2')de’=—j ; 








"\de!=—5 F {C* eos Be +U"}. (8) 


*h 
i I8(2’)K"%(z,z 


In (7) and (8), 





K*(2z,2’)=K,(2,2’)+K,(z,2’), (9) 
K*%(z,2’)=K,(z,2’)—K,(2,2’), (10) 
— ‘Bb. 
U*=U cos (F sin ®) (11) 
U*=j7U sin (F sin ®), (12) 
1s__AtCl, , 

C wane, ead (13) 
(r=, (14) 


The integral equation for the unloaded current in 
an isolated receiving antenna consisting of a single 
conductor is identical in form to (7) and (8). It 
follows that solutions for the currents J%(0) and 
I#(0) may be obtained from investigations of the 
receiving dipole by appropriate substitutions for 
the kernel and the U-function. These currents 
may also be obtained from the relations 


2Bh3U cos (5 sin v) 














a a ™ 
I;,(0)= Z ies Latin , (15) 
and 
‘ 2Bh2U sin (5 sin ®) 
2ph2U* 3 
ag ee 


where 2h$ is the symmetrical effective length of 
either antenna; 2h% is the antisymmetrical effective 
length of either antenna; Z,; is the self-impedance of 
either antenna in the presence of the other, and 
Zi. is the mutual impedance. Equations (15) and 
(16) follow directly from the equivalent circuit of a 
receiving antenna, for the unloaded condition. 
The no-load currents J;z(0) and J,2(0) may be 
obtained from (5) and (6) using (15) and (16). 


These are 
cos (5 sin ®) 


» ZatLie 


. sin (# sin ®) 
he | ® 2 
0351 2-8, ‘ 


I,z(0)=—4rU 








(17) 
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and 


hi cos (F sin o ) 
nN ZatZLiz 


he sin (5 sin ») 
~J nN Zai—ZLiz 


Ix (0)=—4eU 





(18) 


Second order mutual- and self-impedances are 
available in the literature [King, 1956, ch 3]. If 
the symmetrical and antisymmetrical effective 
lengths can be found to the same accuracy, second 
order no-load currents can be obtained, using (17) 
and (18). It turns out to be much simpler, from the 
computational point of view, to determine the sym- 
metrical and antisymmetrical effective lengths to 
second order accuracy than to evaluate numerically 
the solutions of the integral eqs (7) and (8) to 
obtain second order short-circuit currents. A second 
order solution for the scattering error in direction 
finders is sought in this paper. 

The effective lengths 2h{ and 2h% are like 2h, for 
an isolated dipole except that the expansion param- 
eters are Y and y instead of ¥x;. These differences 
should be expected from an examination of the 
integral eqs (7) and (8). The various expansion 
parameters are defined as follows when ph=7/2: 


—_ Y r é r | 
V=|0.($0)+0($0)) (19) 
[King, 1956, p. 268, eqs (12) and (14a)], 
| nN r | 
weal £5 ie a 2 the | 9 
v=|C.(}0)-C ($0)| (20) 


[King, 1956, p. 273, eq (4a)], 


Yn=|C.(}0)| (21) 


[King, 1956, p. 102, eqs (27) and (29a)]. 
Here 


Cy (> 0) ~2 sinh! (¢)-G r—] Siz, (22) 


C, (}0)=Ciap—Cisq—jSisp+jSis, (23) 


Ci X=0.5772+In X—CiX. (26) 
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An important relation is 


Vri=2—k (27) 
where 
2h 
o=2in( —~ 92 
Q 2In( - ), (28) 


and k=1.445 when 7<2<20 for bh=> (King, 1956, 


p. 105, fig. 20.7]. It is reasonable to suppose that 
y'=2'—k, (29) 

and 
y=0"—k. (30) 


[King, 1956, p. 270, fig. 3.3]. 
From (27) and (29), 


Y=0+yY'—Yr, (31) 
and from (27) and (30), 
0° =0+-y*— x1. (32) 


Q* and Q* may now be computed from (31) and 
(32) using (19), (20), (21), and (28) when Bh=7x/2. 
The values of |Aé|/A and |h2|/A are then read directly 
from figure 2, a curve relating |h,|/A to Q@ when 
Bh=7/2. It should be mentioned that in general 


he=hi' +jhi. (33) 


Reference to the literature [King, 1956, p. 488, fig. 
9.2a] shows that h;~0 for all values of 2 when 
Bh=7/2. It is a legitimate procedure, therefore, to 
replace hé and h? by |As| and |A?|, respectively, in (17) 
and (18) provided the theory is restricted to the case 


Bh=n/2. 
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FIGURE 2. Effective length Rs for a half-wave dipole 


as a function of 2(Bh=7/2) 











In table 1 are listed second order values of Z,. and 
Z,, as functions of 6/\ when Bh=z/2 and Q=10 [King, 
1956, p. 305, tables 8.10 and 8.11]. Table 2 pro- 
vides computed values of Y*=1/Z* and Y*=1/Z* for 
Bh=x/2 andQ=10. This table facilitates error com- 
putations when @¥72/2. Table 3 furnishes computed 
values of y’, ¥*, 2’, and Q* and the values of |Ai|/A 
and |h%|/A, as read from figure 2 when Q=10 and 
Bh=-x/2. These data permit evaluation of J;z(0) and 
Isz(0) from (17) and (18), as functions of ®, to second 
order accuracy. 
































TaBLe 1. Mutual and self-impedance for coupled antennas 
2h 
a=2n( =) =10, Bh== 
(All values are in ohms) 
b Zi2 2a 
A | 
0.1} 75.1-j7.8 82.0+j 37.7 
.2 | 55.2-$33.7 81.44) 38.7 
3 29.9—j 46.3 82. 0+] 39.9 
4 | 4.9-j43.2 | 84.6+741.7 
5 —15. 2—j 30.4 85. 8+) 42.7 
6 | -26.2-j125 | 86.2+742.9 
.7 | —27.0+74.9 86. 1+) 42.9 
-8 | —19.5+ 17.1 86. 2+) 42.8 
‘9 | —7 247 21.6 86. 3+) 42.8 
1.0 5. 14+j 18.2 86. 4+] 42.9 
TABLE 2. Symmetrical and antisymmetrical admittances for 
coupled antennas 
rT 
=—, 2=10 
o 2 
b = ae pi ae 
: | my 
0.1 0. 6143—j 0. 1169 0. 3258—j 2. 149 
- . 7311—j 0.0268 - 4420—j 1. 221 
3 - 8909+-j 0.0510 . 5136—j 0. 8498 
4 1.112 +] 0.0187 . 5878—j .6261 
5 1.375 —j 0. 2395 -6498—j .4703 
-6 1.326 —j 0. 6720 .7160—j .3529 
on 1.023 —j 0. 8273 -7946—j . 2670 
.8 0. 8299—j 0. 7453 .8933—j . 2172 
9 - 7603—j 0. 6190 1.017 —j .2307 
1.0 . 7559—j 0. 5048 1.126 —j .3421 














*Multiply all numbers by 10-2 


TaBiEe 3. Computed values of the expansion parameters and 
the symmetrical and anlisymmetrical effective lengths for 


ph=5 and 2=10 














: y ya Qe Qa [heli |nclin 
‘Fee Mee AiR ae AGRE 
| ae 
0.1 | 10.55 6.78 12.00 8.23 | 0.1884 | 0.2077 
21 306) 2 10. 63 9.57 1944 |. 2000 
34 ea). - 88 9. 86 10. 26 1984 . 1964 
21 ae 9.12 9.48 10.57 2005 . 1948 
5 7.94 9.17 9.39 10. 62 .2019 | = 1944 
6] 809 | 9.02 9.54 10. 47 2002 |. 1955 
; oe 8.77 9.81 10.22 | .1988 | .1965 
8 8. 62 8.52 10. 10 9.97 | .1975 | .1979 
‘So: 4 Ue 8. 32 10. 25 9.77 | .1964 | .1990 
1.0 8.87 8.25 10.32 9.70 | .1961 |. 1995 
' 














3. Circuital Relations in Coupled Receiving 
Antennas 


Given the short-circuit currents [,z(0) and J2z(0), 
the problem is that of finding the voltage V,; and 
V., that is developed across identical load impe- 
dances Z;, connected across the terminals in antennas 
1 and 2, respectively. Writing Viz,=|Viz! /diz and 
Vor=|Voz| /dor, the angle sought is ¢6=¢:z—¢2;. The 
angular error in the interferometer due to antenna 
coupling or scattering is then 





56=¢—8b sin ®. (34) 


The terminal currents in a two-element transmitting 
array immersed in an electric field may be obtained 


by superposition, as indicated in figure 3. Clearly 
11(0)=he0)+ hv(0) (35) 
and 
T21(0) = Io2(0) + hv (0) (36) 


Here J,2(0) and J2,(0) are the terminal currents 
when the antennas are excited only by the incident 
electric field, and the impedanceless generator volt- 
ages V, and V, are zero; J,y(0) and J,,(0) are the 
terminal currents when the incident electric fields F, 
and F, are zero and the antennas are center-driven 
by the voltages V,; and V2. 

The fields, currents and voltages shown in figure 3 
are separated into symmetrical and antisymmetrical 
parts as follows: 


E,=E+E (37) 
E,=E'—E* (38) 
11(0)=L,(0) + 70) (39) 
Tx4(0)=1i(0) — 170) (40) 
Vi=V4+Ve (41) 
,=V—Ve (42) 


This separation permits the symmetrical problem 
(superscript s) portrayed by figure 4 to be discussed 
independently of the antisymmetrical problem (su- 
perscript a) portrayed by figure 5. 


Symmetrical Problem 


Refer to figure 4. The following relations are true: 
T;(0)=I;(0)+17(0), (43) 
Ix(0)=5, (44) 
VAL 
Vi=—I7(0)Zz, (45) 
B=ZyntZLie. (46) 
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The principle of superposition applied to obtain the terminal currents in a 
two-element transmitting array immersed in an electric field. 
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The symmetrical problem. 
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The antisymmetrical problem. 
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Hence 
ZatZi . (47) 


1i(0)=15(0) tat, 


Equation (47) establishes the relation between the 
symmetrical part of the load current and the sym- 
metrical part of the short-circuit current. 


Antisymmetrical Problem 


The relations (43) to (45) apply 


Refer to figure 5. 
The antisymmet- 


with superscript @ written for s. 
rical impedance is required. It is 


L!=Z2—Zi. (48) 
It follows that 
Z1—Z 
HO= TO 777} a 


Equation (49) establishes the relation between the 
antisymmetrical part of the load current, and the 
antisymmetrical part of the short-circuit current. 

Figure 6 illustrates the two-element receiving an- 
tenna array with identical impedances Z;, connected 
to the load terminals of each antenna. The load 
currents in antennas 1 and 2 are given by (39) and 
(40), respectively. The voltage drop across Z, in 
antenna 1 is then 








Viu=Zz {Ii (0)+J7(0)}, (50) 
and the voltage drop across Z; in antenna 2 is 
Vir=Z, {Ti (0) —I7(0) }- (51) 
ANT | 
ANT 2 
I,,0) = T0)- IO) | X,(0) = TO) + TO) 
Vp. = I,(01Z, Vz = ,(00Z, 
h ws 
| Vie $2, 
ss 
pe 
$7, Vat 
<2 —S~—OS 
" ns 
—— 
Ficure 6. Two-element receiving antenna array. 


Impedances of value Zz are connected to the load terminals, 








7i(0) and J2(0) are expressed in terms of J%(0) 
and J%(0) by (47) and (49). Also, from (5) and 
(6), #(0)=(1ie(0) + L22(0))/2 and J%(0)=(hez(0) 
—Ip(0))/2. It is now a simple matter to obtain 
Viz and V2;, as given by (50) and (51), in terms of 
Tiz(0) and J2,(0). The final expressions are 


Viz Z, 
(Za tZitZz) (Za—Zi2tZz) 








{Tie (0)(Z32,+2Z1Z,—Z}.]+e2(0)Zi2Z,}1 (52) 
and 
Vez= Z1 
(Za tZitZz) (Za—Zi2tZz) 
{Ize (0) (Zi, +ZnZ,—Zj2]+ Ne(0)Z2Z,}- (53) 


These equations provide a way of computing the 
voltages across the identical load impedances in 
terms of the short-circuit currents J,,(0) and J2,(0). 

One check on the validity of either of the above 
expressions is to let b> so that Z,.—-0. Then (52), 
for example, becomes 


1e (0) Zs, Zz 


Dnt he ai 


, 
Vii= 





Here the primes indicate that the various parameters 
change value when one of the antennas (in this case 
antenna 2) is moved to infinity. Equation (54) is 
recognized immediately as the correct expression for 
the voltage across the load impedance Z, of an 
isolated center-loaded receiving antenna. 

Let the following definitions be made: 


Zr = b 
(Za t+Zw+Zz)(Zi—Zi2t+Zz) 
Z23:+ZnZp—Zi,= Ae s 
Z12.Z,= Be® 
I,2(0)=Cha(0)e” j 











With this notation (52) and (53) become 
Viz=T' L(0) { Ae**+ BCe?4*” } 
Voz=T'11(0) { ACe**t” + Be#}. 


Using these expressions and the formula 


—1 » —14,— -1 zy 
tan~* z—tan~* y=tan {; x 


it is readily deduced that the angle ¢ indicated by 
the interferometer is 


(B?—A?)C sin y+ AB(1—C?) sin (a—B) 
(B?+ A*)C cos y+ AB(1+C?) cos (a—8) 
(59) 


(58) 








¢=tan-! 
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4. Numerical Analysis 


Figure 7 exhibits a second order curve for 6 as 
a function of 86 for an isolated interferometer when 
Q=10, Bh=/2, 6=7/2 and Z,=100+) 0.0 ohms. 
To facilitate computations for other values of ® and 
Zz, when 2=10 and Bh=7/2, tables 2 and 3 are pro- 
vided. Table 4, giving values of A, a, Band Basa 
function of b/d is useful when constructing curves of 6 
against 6b with as parameter. This table applies 
only when 2=10, Bh=x/2 and Z,=100+7 0.0 ohms. 
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Figure 7. Antenna coupling error in a two-antenna inter- 


ferometer. 


2=10, Bh=x/2, P=x/2, and Z_.=100+j 0.0 ohms. 


. b 
TABLE 4. A, a, B and B versus x 


(These data apply only when Bh= 2/2, Q=2ln(2h/a)=10 and Z,=100+/ 0.0 ohms) 




















o Bb A* a | B* B 
x | | 
\_ | et ee ae 
| 
0.1 0. 6283 1. 366 0. 9520 0.7551 | —0.1037 
a 1, 257 1.794 . 8853 .6467 | —.5479 
3 1, 885 1.974 . 7396 . 5512 | —. 9975 
4 2.513 1. 956 . 6371 .4348 | —1.458 
5 3.142 1. 826 . 6245 3399 | —2.034 
6 3.77 1. 757 . 6783 . 2903 | —2. 696 
4. 398 1.801 . 7252 . 2744 2. 962 
a 5. 027 1. 875 . 7174 . 2593 | 2. 422 
9 5. 655 1. 893 . 6849 . 2277 1. 893 
1.0 6. 283 1, 857 .6687 | 13990 | 1.300 
| 





*Multiply A and B values by 104, 
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5. Concluding Remarks 


It is of interest to observe that when 6= 72/2, 5=0 
for Bb=nz, n=0,1,2,... Evidently, when 6¥7/2 a 
spacing 8b exists such that 6=0. 

The error curve shown in figure 7 applies to base- 
loaded vertical antennas over a perfectly conduc ting 
earth when the load resistors are halved.. If all im- 
pedances are divided by 2, a, 8, y, and C remain 
unchanged, but A and B are one-fourth their free 
space values. Accordingly ¢ (and therefore 6) are 
the same as for the isolated array. These statements 
are readily verified using (55) and (59). 

A precise theory has been presented for the scatter- 
ing error in a two-antenna direction finder. It 
applies with equal accuracy to radio interferometers 
employing antennas of length other than h=)/4 
(but with slight increase in complexity). Further- 
more, the theory is readily generalized to include a 
wave tilt with respect to the antenna axes. 


The writer is indebted to R. W. P. King of Har- 
vard University, and R. H. Duncan of New Mexico 
State University, for reading this paper. He wishes 
to thank Mrs. Dixie S. Baca for doing the typing. 
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The Electrically Short Antenna as a Probe for 
Measuring Free Electron Densities and Col- 
lision Frequencies in an Ionized Region ' 


Ronold King,? C. W. Harrison, Jr.,* and D. H. Denton, Jr.* 
(February 19, 1961) 


If the admittance of a missile, satellite, or drone-aircraft antenna is monitored as the 
vehicle traverses an ionized region, it is possible to determine the free electron density and 
the collision frequency of the region if theoretical relations between these quantities are 
available. In this paper formulas are developed that relate the admittance of an electrically 
short center-driven dipole or a base-driven monopole when immersed in a conducting dielec- 
tric to the effective dielectric constant and conductivity of the medium. From well-known 
formulas relating these quantities to the free electron density and the collision frequency 
of an ionized region, these !atter may be determined directly from measured admittances. 
The results obtained when the antenna is treated as a lumped capacitor are considered. 
It is shown that when the conductivity of the medium is increased to a value that is still 
quite small, the effect of radiation on the input admittance becomes negligible. The 
electrically short antenna immersed in sea water is discussed briefly. 


1. Introduction 


It has been suggested that electron densities and collision frequencies in the ionosphere 
may be determined by measuring the impedance of a base-driven electrically short missile 
antenna, using a lumped-element impedance bridge, as the vehicle traverses the ionized region 
[Jackson and Kane, 1959 a and b; Haycock and Baker, 1959]. A theoretical relation between 
the reactance of the RF probe and the electron density then permits the evaluation of the 
electron densities from the measured reactances. It is the purpose of this paper to derive an 
approximate expression for the admittance of a center-driven dipole or of a base-driven probe 
antenna in a infinite ionized homogeneous medium. The admittance so computed is approxi- 
mately the same as for an antenna in a medium of finite extent, provided its boundaries are 
sufficiently far away. It should be mentioned that the impedance of a stub type antenna 
protruding from a missile differs negligibly from the impedance of the same antenna when 
mounted on an infinite ground plane, provided the dimensions of the carrying vehicle are not 
too small in terms of the wavelength. 

The formula for the admittance is derived from the integral equation for the distribution 
of current along the probe by means of a technique that achieves sufficiently high accuracy in 
the leading terms of an iteration procedure so that higher order terms are not required. The 
leading terms in the current consist only of linear and parabolic functions with complex coeffi- 
cients. A method for finding these coefficients is described. No account is taken in the 
theoretical development of residual charges on the missile [Kraus and Watson, 1958; Jackson 
and Kane, 1959a; Krassovsky, 1959]. 


2. Integral Equation for the Current 


Figure 1 illustrates a cylindrical electric probe of length 2A and radius a, center-driven at 

=0 by an idealized delta function generator of emf V. The axis of the antenna coincides 

with the z-axis of a cylindrical system of coordinates 7,6,z._ The quasi-one-dimensional theory 
employed requires the following inequalities to be satisfied : 


2 
h>a, Ba mit (1) 


! Contribution from the Sandia Laboratory, Albuquerque, N. Mex. 
2 Consultant to the Sandia Corporation and Gordan McKay Professor of Applied Physics, Harvard University. 


3 Member of the Scientific Staff, Sandia Corporation. 
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Figure 1. Probe antenna in a conducting dielectric. 


where \ is the wavelength in the medium. The antenna is assumed to be perfectly conducting 
and in contact with the dissipative medium characterized by uy, ¢€, and o, the absolute perme- 
ability, dielectric constant and conductivity, respectively. 

The integral equation for the current along a perfectly conducting symmetrical center- 
driven antenna of finite radius immersed in an infinite dissipative medium of arbitrary attenu- 
ation is [King and Harrison, 1960} 





m . Ank Bins ae 
= A,(2)= * L.(2’)Kp(2,2’)dz’= — sat |g V Mott F,, | (2) 


The only restrictions that have been imposed are given by (1). 

In (2), /,(2’) is the unknown complex current in the element dz’ at a distance 2’ from the 
origin. The distance along the antenna from the origin to the point on the surface of the 
antenna, where the vector potential A,(z) is computed, is z. The kernel Kp(2,2’) is 

eR 9 —IKRp ; 
(2 a = — a ’ (3) 
‘ R R, 


where R= y(z—2’)?+a*?, Ria= y(h—2’)?+e’, and k is given by the expression 


k=B—ja=wyue(1—jp)=wy ue [ f(p)—Jjg(p)); (4) 








w is the radian frequency, f(p)= 6 sinh~! p) and g(p) sinh! p), where p=a/we. 
Tables of the function f(p) and ste are eueaiies in the literature [King, 1945]. wu—k¢, where 
¢ is the characteristic impedance of the dissipative medium given by ¢=[jwu/(o+jwe)]!” 


=f,/[1—j(a/8)|. M,.=sin k(h—|2|) and /,=cos kz—cos kh. The function U is defined by 








e IkRh 
ee an 1.(2") de! =—j2 A). (5) 


4nk 


Equation (2) is applicable to an antenna of any length. An RF probe is defined for present 
purposes as an antenna for which 6h <0.3 radians. For such an electrically short structure, 
the several functions may be expanded in powers of kh since it is postulated that 


Bh<l, ah<l. (6) 
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4 Atal 


The kernel Kp(z,2’), as given by (3), may be written 


- | , m 
Ko(2,2" ROH e- —R,) +35 (Re Rj). (7a) 
Also, 
— lol) — he (RI ol _Ra—lal)* , lel keh? 2|2| 2 \ , 
Moe=sin k(h—|2|) =k (h—|2|) =A 2 ata 1—F) (1 [:- : +5) - () 
and 


2},2 4 2 x ~2 
Fy,=cos kz—cos ha (1 5) (- Ra i) => (1-@) (1-45 | iti :]) ore 


At the driving point (z=0), MM, and Fy, have the amplitudes Skh and Ckh?/2, in the order 
named, where S=1—(k?h?/6) and C=1—(k*h?/12). It is assumed that [1—(|z2//A)] and [1— 
(2?/h?)] adequately represent the variation, respectively, of M4, and Fy, along the probe. Ac- 
cordingly, the simplified representations of 14, and F), are 





My. = Skh (1 —*), (Sa) 
Fu~3 keh (1-5) (Sb) 
The function U becomes 
° (7- “7, wi m4. 
and 
cos pied (10) 


In these expansions terms of the order of magnitude k*h* have been neglected. 
g g 
The final form of the integral equation for the current along a probe antenna in a con- 
ducting dielectric is obtained from (2) with (6), (8), and (10). It is 


Ps) 
‘. 12 aR al (RR n+5% (RR) | ae’ 
a (1) ae i-)+uckh (1-5) } (11) 


l~ 


3. Some Remarks Pertaining to the Solution of the Free-Space Antenna Equation 


King and Middleton [1946; see also King, 1956, Chapter II] begin the process of obtaining 
an approximate solution of the free-space antenna equation 


4 h ee | oe ae 
 4()={- T(z’) : PR dz’=—} «lc, cos Be+5 sin ala\| (12) 





by expressing J,(z) in terms of a reference current J,(2)) and a distribution function f(z) that is 
unknown. Let J,(2z)=J,(2)f(z), so that J,(2’) a zo)f(2’). Then I,(2’)=1,(2)[f(2’)/f(2]= 
I,(z)g(z,2’). A function ¥(z) is now defined in terms of g(z,z’) by the expression 


*h 2 —IBR 
Wa=J 9(e,2") Ge de Wa) +102), (13) 
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where 2; is a suitable reference point; namely, z;=0 when h<)/4, 2,=h—(A/4) when h> d/4. 
With (13), the integral in (12) may be written in the form 


# A={ bi) © dy 1d +O), (14a) 
M —h u 
where 
*h — jpR 
D(2)=1(2)v(2) +] [ 2) L2g(,2") | Cae PR dz’. (14b) 


It is evident that (14a) is an identity regardless of whether g(z,2’) is the correct distribution 
function or not. However, the better the choice of g(z,2’) the smaller is the difference function 


D(z). If g(z,2’) were precisely correct, the integral in (14b) would vanish, and y(z) as defined 
in (13) would become 





<a (15) 


Nevertheless, D(z) in (14a) could not vanish simply because the ratio of the actual vector 
potential to the actual current is only approximately equal to the constant ¥(z,). The difference 
term D(z) is evaluated by successive approximations and its contribution to the expression for 
I,(z) is contained in the higher-order terms of the solution by iteration. 

Owing to the analytical complexity of the higher-order terms in the King-Middleton 
solution and the fact that its leading trigonometric term includes only a very rough approxi- 
mation of the susceptive part and no contribution whatsoever to the conductive part of the 
current, it is desirable for practical purposes to derive an approximate solution for the current 
in which the leading trigonometric terms are a much better approximation. This may be 
accomplished if the vector potential and the current are separated into components and these 
are then properly associated to permit the introduction of several expansion parameters of 
the type ¥(z2,) in (13). Each of these is the approximately constant ratio of a particular com- 
ponent of vector potential to its associated current. Several integral equations are obtained in 
this way which may be iterated independently. A study by King [1959] has shown that the 
sum of only the leading trigonometric terms of these component currents with suitably defined 
complex coefficients yields the total current and driving point impedance with reasonable 
accuracy. A precise solution for the distributions of current and charge, and the impedance of 
an electrically short antenna in free space has been given by King [1952; see also King, 1956, 
ch. II, sec. 31], using a method that is not particularly convenient when the antenna is in a 
dissipative medium. 


4. Solution of the Integral Equation for the Current Along an Electric Probe 
in an Ionized Medium‘ 


In the King-Middleton antenna theory, the leading reactive term in the complicated 
series for the current distribution is the familiar sinusoid. Only one g(z,2’) and its associated 
¥(z) are used in solving the problem, but several iterations are required to achieve accurate 
results. On the other hand (11) is arranged in such a way that the leading term consists of 
two components; one varies as [1—(|z|/h)] and a second varies as [1—(z?/h?)]. The linear 
distribution may be considered to be maintained principally and directly by V (which is localized 
at z=0); the parabolic distribution is induced by the retarded field distributed along the 
entire antenna and expressed in terms of the equivalent voltage U (as in a receiving antenna 


4 The procedure employed in this section to solve the integral equation is discussed in somewhat greater detail in Sandia Corp. Tech. Memo. 
No. 53-60(14), ‘“‘ Approximate solution of the free-space antenna equation based on a determination of the complex coefficients of the leading com- 
ponents and current,’ February 1, 1960. 
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situated in a uniform field). Since the leading terms consist of two components, which are 
different functions of z, these may be expressed separately in terms of the appropriate distri- 
bution functions gy(z,2’) and gy(z,2’). These relat*ve distribution functions prove useful 
in the definition of two expansion functions yy(z) and Yy(z). In order to introduce these func- 
tions, it is necessary to separate the total current along the antenna and the vector potential 
on its surface into component parts, and associate these parts according to their analytical 
dependence on z. The ratio of a particular component of vector potential to its associated 
current is the constant ~y(z,) or yu(2,). Here z, and 2, are the values of z where the given 
part of the total current has its maximum value. 

The form of the right-hand side of (11) suggests that it should be possible to approximate 
the current along an electrically short antenna in a conducting dielectric by 


1,(2’)= AM,» + BF,, (16) 


and evaluate the complex coefficients A and B in such a way that when (16) is substituted into 
(11), the equation is approximately satisfied. With this substitution, 


h h 
A il Mur Ko(2,2')d2’ +B { Foe Ko( 2,2") 


_; 2h fy [2 iF }. ‘ 
Pan s(1- = )+UCKR (1 s) fF (7) 


The integrals that occur in (17) are 


:. M,.: (R-7 rakes (18a) 
r F (p—a ee’ (18b) 
ae ar 
f. M,,(R"— Rp) dz’, (18¢) 
*h 

| Foe(R"—R2dz’, (18d) 
-h 


where n=1 or 2. When these have been evaluated it is found that the leading term in (18a) 
varies approximately as 14), and the leading terms in (18b, ¢c, d) vary approximately as F.. 

Equation (17) may now be separated into two parts: One part has an MM), variation; the 
other an Fy, variation. Thus, 


nV SKA z 
A Mo. {i a} d'= =U ay i 
and 


h k? | ro , Y a ie 
Al” Mw —5 (R—Ry)+5-— (RF) d+ BY Fed RoR 


a Fs) 7/113 h2 ~2 
2 (RR) +i 5 (RD } del=j ee { 1-73 }- (20) 
. wut 1 


—*) 
2 
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For convenience let the following notation be used: 


ns , a By 
va=| 9o(2,2’) { Fe paz 


va= | ul2")4 —5 e- Ry) +3 (R— Re) baz 


and 
tare { otz,0?)4 41 —E e—2y gE e—aep bar’, 
du “a “ R R, 2 h 6 h eal 
where . 
to’. 5 
- ie, Se 
J\2,2 uy ( |2| i 
1_—— 
= 28( are) 
Jur 2,2 = > 





Fos -omn(1—2) 
92,4?) Ete OS i) 
u\Z, Fiz (1- 0-2) 


Equations (19) and (20) now become 








, 2rVk 
Ayn= 2a heh? 
mr 
and 
: 4rkU 
Aeot Blau=) = _) 
"on (1—2®) 


respectively.5 
The substitution of (16) in (9) gives 


U=—j pe (Ad (h) + Bhulh)}, 


h Be 
¥a(h)= Shh | = 4 Fe iy GEE st bade 
ae tp 
HRs, ER 
wi =O8F |" (1—Fr ) { sts EF hae 


The use of U from (27) in (26) gives 


where 





and 





B=AT(h), 
where 
3 2h? 
Y, (h) —Wa2 1—— 7) 


eae (aa re 


5 Equations (19) and (20) are valid for any value of z; in particular, they are valid for z=0. 
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(30) 


(31) 





; 








| 
i 
' 


Rte oe aes 


earcest 





Also, from (25), 
QnVk 


desis (1—") 


(32) 





A=j 


A tentative expression for the current is obtained when (30) and (32) are substituted in (16) 


It is ; 

. 2wnVSkh ke ; 

[,(2)=Jj ma (2 ERY ae 7; i )+ 70h) S 39 ni > (= ~i) } (33) 
dl 





where use has been made of the relation wu=k¢. Readers are reminded that in deriving (33), 
all difference integrals have been neglected. The parameter J(h) is evaluated from the 
integrals (21), (22), (23), (28), and (29). It is found that with terms of the order k°h5 





a 
Ya ~2 In (=)—2, 

pa 2hh_; RY_R 

— Ss * Ss Se 


; 2: 
ae 2In (= *) 3423 hj = BR, (34) 


es ig 


Yo(h) ~2kh In 2— jreh—-(5 +?) eh, 





du(h) Ej : eh, 


4 


An approximate formula for 7(h) may be obtained if (34) is substituted into (31) and 
terms of the order (kh)* are neglected. It is 


(21n 2) th 5 koh ee (9 +2 
2 In 3 kent (In 4 is) +35 4 paps 





T(h) = (35) 


In the derivation of (33) terms of the order k*h* and higher have been neglected in the 
individual functions. Therefore, when the current is expressed as a series in powers of kh, 
all terms of order higher than k*h? without the common factor kh, or higher than k*h* with this 
factor, must be dropped. As a first step in this expansion and simplification 7(h) as given in 
(35) may be expressed as follows when terms of the order k*h? and higher are neglected (since 
T(h) occurs only where multiplied by an additional factor kh): 


T(h) = 3) 13 In 2—1)kh—jk*h?), (36) 


2 

3(Q—3) 
2h ‘ 

where Q=2 In 7 It follows that including terms of the order k*h’, 


Ckh 


Th) $5 TO) (1455) See 


gy [(3 In 2-1) kh? jh] (37) 


If (37) is inserted in (33) and only terms up to and including the order k*h? are retained exclusive 
of the common factor, the current (33) is given by 


j2nkhV keh? oe ae | _2 
1,(2)= Wa L0H) 3(Q—3) a) } (38) 











In (37) and (38) use has been made of the relations S=[1— (k*h?/6)], C=[1—(k*h?/12)], C/S= 
[1+ (k7h?/12)], and [1— (kh?/6)]/[1— (k7h?/2)] ~ [1+ (kh?/3)]. The admittance obtained from (38) 


with z=0 is 





. 2xkh kh? (3In2—1)] . FA 
epi Gee 9 E  @—3) J 736-3 — 
These are the final general formulas for the current in and the admittance of an electrically 
short center-driven dipole of half-length h. The corresponding formulas for a base-driven 
monopole of length A erected vertically on an infinite ground plane are obtained from (38) and 
(39) by multiplying the right sides by 2. 

It is now clear why it is necessary to include cubed terms in the small quantity kA. For 
if a=0so that k=6 and ¢=¢, are real, the only contribution to the real part of Y(k) comes from 
the term with kh’ as coefficient. Thus, with Y(k)=@(k)+ 7B(k), the admittance for a dipole 
in air with k=£o,f=f) ~ 1207 ohms, is 


3 ack 2rBih'* . 2rBoh er) 
} (Bo) ~ 3(Q2—3) Coa Jj ton 1+ 3 , (40) 





where 
F=14[(3 In 2—1)/(@—3)]=1+[1.08/(Q—3)]. (41) 
The corresponding expression for the impedance is 


wae CoWar Boh? ee War , 
3 — 3 








For many purposes when £2’ is sufficiently small the correction terms B7h?F/3 may be neglected 
; 0 ; 0 ; 


compared with unity. 
5. Accuracy and Numerical Example When a =O 


As a specific numerical example consider an antenna for which h/a=75, 2~10 so that 
¥a=6.635. With ¢)~1207 ohms it follows that 








Y (B) = {0.120 Bih'+j 2.51 Boh [1+-0.385 Bzh2]} X 10-3 mho, (43) 
19.06? 398.1 
4(6o)=17 0.385 Be! Bh [1+ 0.385 Behe] OMS: (44) 


This last formula may be compared with the following more accurate result obtained by another 
method [King, 1952; see also King, 1956, ch. II, sec. 31]: 


o=18.3 Bih?(1+0.086 Bih?)—7 ae (1—0.383 Bh?) ohms. (45) 
0! 


The leading terms in (44) and (45) are clearly in good agreement. Specifically, when B=0.3, 
the numerical results from (43) to (45) are: 


Y (Bo) =0.972 X 10-*+-7 0.779 X10-* mhos, (46) 
Z(B) =1.60—j 1283 ohms, (47) 
7 )=1.66—j 1274 ohms. (48) 
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The leading terms alone in (44) and (45) give 


7 (Bo) =1.71—7-1327 ohms, (49) 
Zy=1.65—7 1320 ohms. (50) 


It is clear that the corresponding numerical values are in very good agreement. Moreover, the 
simpler formulas with leading terms only yield results that differ by only about 5 percent from 
t the more accurate values. The accuracy of the formulas (40) and (42) has thus been verified 
when a=0 and BA $0.3. 





ee 


6. General Case When a0 


When k=8—ja and rel(1—-J <) are complex because a0, the ratio of conductance to 


susceptance in (39) may differ greatly from its value when a=0. This is a consequence of the 
fact that contributions to the conductance may now come from the leading term with kh as 
coefficient and not only from the term with k*h* as coefficient as when k=8). If the complex 
values of & and ¢ are inserted in (39) and the real and imaginary parts separated, the result is 


Y (k) =G(k) +jB(h), (51) 


orf 2a i+ oh F(1—% i) tangy (1-10 +55, =) } - 
fa LBL 3 3(2—3) “4 


3 9% oe) Bh af, 27, a! ” 
BW= = { sh(1-$)+4 Bn F ( 1 65 +5)— Pear G 105 +5.) } (53) 


and F=1+[(3 In 2—1)/(Q—3)] as in (41). 


In many applications a is small compared with 6 so that with (4) 


a —- Dp — ‘ 
k=p(1-J§) =o we( 1-75 )=wr we(1—J sx) 


Let it be assumed that the condition 
9 
Fil -|= f<: (54) 
WE 


is well satisfied and (51) with (52) and (53) reduces to 
_ oat 37,3 B*h* il 1 4375 eel} Eas 
Y(k)= ot { 2s | a+? Bh Fits, (03) tI | Abts Oh F—350—3 (55) 


Note that (55) reduces to (40) when a=0. 

As a numerical example consider the short dipole, previously considered when in air, 
immersed in a region of the ionosphere where the number density of the electrons is V=1.5 10" 
electrons per cubic meter and the collision frequency is y=1.110° sec~!. At a frequency of 
6 Me/s the real effective dielectric constant ¢ and conductivity o are 


where 


G(k)= 














f Ne? " . 
€= €3€- — €p (1a )=0.068 €o farads/m, (56) 
POR: 2610-7 mho/m (57) 
m(y?-+w*) : 
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where ¢ is the charge and m, the mass of an electron. It follows that B=) ye, =0.1025 radians/ 
= > « 1 — we ; +* 2a is 
m, ¢=o/Ve,=462.3 ohms, a=50f-=7.534X 10-° nepers/m. Since ra 70 X10-*, tie 


simpler formula (55) is applicable. With 8h=0.3, Bh=0.2447. When substituted in (55), 


these values lead to 


Y(k)=1.12 10-*+-7 0.513 X 10-3 mho (58) 
and 
Z(k) =4.26—7 1949 ohms. (59) 


If these results are compared with the corresponding values (46) and (47) for the same antenna 
in air, it is clear that a significant change in admittance and impedance occurs when the antenna 
is moved from air into the specified region in the ionosphere. 


7. Determination of the Dielectric Constant and Conductivity With a Short 
Antenna 


The measured change in the driving-point admittance of an electrically short dipole 
when it is moved from air into the ionosphere or another conducting dielectric may be used to 
determine the effective dielectric constant and conductivity of the region. If the susceptance 
of the antenna when in the ionosphere is B(k), and when in air B(@)), it follows from the 


imaginary parts of (55) and (40) that 


t 2 F793},3 /« "92h 2/< 
Blk) _«Boh+2F BiH /3_ eee (60) 


B(B) Boh FRN/3 U1 FBR 





Use has been made of the relation ¢,={o/ye, and B=Byye,.. Since the terms with 3h? as a coef- 
ficient are small, the susceptance ratio differs only slightly from the value e,._ If e,~B(k)/B(B») 
is taken as the approximate value for substitution in the small correction term, the corrected 


formula for e, becomes 
Bik) 1 nsf, BO ) 
«Big LI+3 Ha (1 Big) f (61) 


where, from (42), F=1+[(3 In 2—1)/(Q—3)]. When £,A is sufficiently small, the leading term 


in (61) may be adequate. 
If the admittance of the antenna when it is in the ionosphere is G(k)+7B(k) and the con- 


ductance in air is G(B)), it follows from (55) and (40) 





2a _ o i G(k)—!? G(Bo) g (62) 
B ile B ! 1 7 - 
(k) ( 1+5 Fe 


; 1 we ; : 
As in (61) the small term 3 F*h? may be negligible when Bh is sufficiently small. 


8. Determination of Electron Density and Collision Frequency 


The formulas (56) and (57) that relate the dielectric constant e=«e, and the conductivity 
o of the medium to the electron density N and the collision frequency v may be solved for v 


and N. Thus, 
i lees (63) 
N= = oe Al id (64) 
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If eee, is obtained from (61) and o from (62), v and N may be evaluated. Thus, both the 
number density and the collision frequency may be determined from measurements of the 
impedance of a short dipole successively when immersed in air and in the ionosphere, provided 
the relations (56) and (57) adequately represent the properties of the latter. 


9. Static Capacitance Formula 


Since it has been assumed that the antenna is electrically short, its admittance may be 
determined approximately from its lumped capacitance. When the antenna is immersed in 
air, the approximate susceptance is given by the leading term in (40). That is, 





2rBoh Qreoh sna 
B =wl' = = ( (6: 
(Bo) w Coal @ Vai ’ dD) 
so that with (34), the capacitance is 
J 

C=. (66) 

in=—1 

a 





Qreoh 
<o. (67) 


The susceptances and reactances obtained at a frequency of 6 Me/s with 8h=0.3, and h/a=75 
are, from (67), 


wC>= 0.666 < 10-3 mho, —, =—1502 ohms; (68) 
0 
irom (66), 
wCO=0.754 107° mho, — =~ 1326 ohms. (69) 


More accurate values according to (46) are 
B(Bo)=0.779X10~* mho, — Xo(8))=—1283 ohms. (70) 


It is seen that the use of the conventional formula (67) leads to an error of nearly 25 percent 


in the susceptance of the short antenna in free space. 
If the antenna is immersed in a medium with relative dielectric constant ¢,, the appropriate 


expressions for the capacitances C, and C. are simply (66) ard (67) multiplied by e,. It 
follows that 
Ce Cn 
fe ee (71) 
( 0 
This is the same value given by the leading term in (61). Thus, the relatively larger error in 
susceptance that results from the use of the formula (67) for the static capacitance is canceled 
out when the ratio (71) is evaluated. It is this ratio that was used, for example, in the meas- 
urements reported by Jackson and Kane [1959a]. Note that the above discussion applies 
when (54) is satisfied. 
10. More Highly Conducting Medium—Plasma and Sea Water 
The generality of the derivation of (39) and (51) has been limited by the conditions (6), 


Bh<1, ah<1, but the relative magnitudes of a and 8 have not been restricted. In the 
important special case corresponding to sea water, a=8, the leading real and imaginary terms 


9 
Y(k) one {i —J 3 FB } (72) 
dl 5 


in (51) are 


CM, 
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where F is defined in (41). The appropriate formulas for ¢, and 8 are obtained from the 
defining relations, with a=, ¢>>we. Thus, 


-_f,_;%][ sow 7? ff 2en]? a 
n=[1-i5 |] LST (73) 


1/2 
k=olu(e—jo/w)}!/2~[—jopo}”*=(1—J) | | > (74) 
2 
so that 
1/2 
a=p=|"] . (75) 


With Y(k)=G(k)+7B(k) as given in (72), it follows that 


9 
pee, (76) 
Var 
So fo 
B(k)= 3a WU (77) 


The susceptance of the same antenna in air is 


2r 2Qrh 





3 —— —=—— WE. 78 
E (Bo) ton Bh Van WEp (78) 
It follows that 
o  G(k) gi 
wen B(B) (79) 


This formula may be used to determine o from measured values of the conductance G(k) of 
the short antenna in the conducting medium with a=8 and the susceptance B(@)) of the same 
antenna in air. 

Note that the radiation term with the coefficient 6*h*/3(Q—3), which constitutes the 
principal contribution to the conductance when the antenna is in air, does not appear in (72) 
since it is negligible to a high order. Thus, the entire admittance is independent of radiation 
and the dipole behaves like a pair of electrodes at low frequencies. Moreover, the conductance 
is very great compared with the susceptance, which is a small quantity of higher order. 

The admittance of a dipole immersed in media for which the attenuation constant is in 
the broad range 0<a<8 may be considered in three ranges as follows: 

1. 2a/B< <6°h?/3(Q—3)<<1. The dipole is essentially a radiating structure with a conduct- 
ance that is very small compared with the susceptance. 

2. Bh 3/3(Q—3)<<2a/B<<1. The dipole is like a pair of electrodes with negligible radia- 
tion. The conductance and susceptance may be comparable in magnutide. 

3. a/B~1; PA<<1. The dipole is like a pair of electrodes with negligible radiation. The 
conductance is large compared with the susceptance. 

Graphs of the input admittance and impedance of an antenna with Q=10, Bh=0.3 are 
shown in figures 2 and 3 as functions of a/@. 


11. End-Effects and Admittance 


Throughout the analysis it has been assumed that the antenna is center-driven by a so- 
called delta-function generator which provides a discontinuity in the scalar potential at the 
center z=0 of the antenna. The antenna consists of a perfectly conducting thin tube with 
adjacent knife edges at z=0 along which an electric field E,(z)==V6(z) is maintained where 
5(z) is the Dirac delta function. The entire tube with its idealized generator is immersed in 
the dissipative medium under study. Strictly speaking, the generator is short-circuited by an 
infinite capacitance and an infinite conductance at the knife edges, so that the current must be 
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infinite. It has been shown by Wu and King [1959] that this large eurrent that is associated 
with knife edges exists only in an extremely short distance on each side of these edges. More- 
over, if the current is approximated by a series of continuous functions, this part of the current 
can contribute nothing if the solution is not of too high order. Thus, the finite current obtained 
in this paper is the total current minus the infinite current at the driving point. 

An actual antenna used for probing a medium such as the ionosphere may be a dipole 
center-driven from a two-wire line, or a half dipole base-driven from a coaxial line that pierces 
a conducting plane. The current entering such an antenna is well approximated by the current 
supplied by a delta-function generator (minus the infinite knife-edge current) and a small 
reactive current characteristic of end-effects and coupling effects near the particular junction 
region between antenna and feeding line [King, 1956, pp. 33-69]. If the line spacing is very 
small compared with the length of the electrically short antenna, it is necessarily an extremely 
small fraction of a wavelength. Hence end and coupling effects are negligible. 


12. Conclusion 


Formulas have been derived for the admittance of a perfectly conducting electrically short 
dipole immersed in an infinite medium that is characterized by a conductivity o and a dielectric 
constant ¢ for a range of ratios o/we that extends from zero to values that are large compared 
with unity. The application of the formulas to the determination of both o and ¢ from measured 
values of the admittance of the antenna when immersed in the medium and when immersed in 
air has been discussed together with the evaluation of electron densities and collision frequen- 
cies of media in which these quantitites may be related directly to o and «. It is to be noted 
that the more complicated problem of a dipole near the interface between two media has not 
been treated. 
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Of the various mechanisms put forward in recent years to explain long range tropospheric 
radio propagation the most important seem to be scattering from atmospheric turbulence 
and partial reflexion from high level subsidence inversion layers. In this paper the writer 
extends earlier theory appropriate to the case of a single elevated inversion layer to cover the 
case of multiple layers. In some cases it is found that the theory and formulas can be applied, 
with very little modification, to the latter. The special cases in which this can be done are 
worked out in detail and it is found that, if certain conditions are satisfied, several weak high 
level inversions can produce a similar effect on the propagation to a single strong inversion. 


1. Introduction 


Of the various mechanisms put forward in recent years to explain long range tropospheric 
propagation the most important seem to be scattering from atmospheric turbulence and partial 
reflection from high level subsidence inversions. In 1952 the writer, after first reviewing earlier 
experimental work, showed theoretically that a well developed inversion could cause the phe- 
nomenon. In a three part paper in 1955, after first reviewing the major theories which had 
been put forward in the interim, he considered in greater detail the special characteristics of 
the field under a high level subsidence inversion (part II). It is now fairly generally agreed 
that the major mechanisms responsible for the long distance fields in tropospheric propagation 
are scattering from atmospheric turbulence and partial reflection from high level subsidence 
inversions. 

Refractometer measurements have shown that several minor high level inversions can 
coexist at different levels and that this situation is probably present for most of the time. It 
would therefore seem desirable to make the previous investigations complete by extending the 
theory to take account of the case when several such inversions are present. This is what will 
be attempted in the present paper. 


2. Propagation Under Several High Level Inversions 


It is possible to show [Northover, 1952] that, provided the distance of the field point from 
the exciting source is small compared with the radius of the earth, the Hertzian function II 
representing electromagnetic disturbances excited by a horizontal electric dipole (or system of 
dipoles) in a stratified atmosphere outside a spherical earth is expressible in the form 


21 2n(kr) P,(cos 8) (1) 


1 
r (n) 
where 2, (kr) is a solution of the equation 


n(n+1) 
r 





a+ { K(é)— 2,(:)=0,  £=2nr)d (2) 


1 Contribution from Department of Mathematics, Carleton University, Ottawa, Canada, 
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and that this representation is valid whether the dielectric constant « or the refractive index p 
is a single analytic function of height (as for a surface inversion) or has different forms for 
different regions (as when high level inversions are present). Further, it appears that, subject 
to the same restriction on range, the condition at the surface of a perfectly conducting earth 


can be taken as II=0. 
Consider the case when several high level inversions are present. Let there be p of them 


a 


and let the boundaries of the jth be given by £=y;, §=2,;(2;>y,;); let also 
2o(=2) Yar <ya<ee 2s SYp<p<.- (3) 


Let the dielectric constant below the lowest? be x, and above the highest be unity; let also 
K=Ky IN 2,<CE<Yo, K=Kg IN 22<CECY3 . . . K=K; IM 2j_1CE<Yy; and K=ky iM 2p) <CE<CYp. 

These values of « refer to the space between the layers; throughout a}given’ layer « is 
supposed to change smoothly, e.g., through the second layer « changes smoothly from x at 
the bottom (=y2) to xz at the top (=z). See figure. 








CWQe 








Wr 
DWF 


Remembering (1) we assume that the secondary field below all inversions (i.e., in the 
space r<é<y,) is of the form 








K 
K 
-' 2(2n+1) { PagnlEr Ku) +G ann (Ev xi) }P,,(cos 6) (4) 


where 


"Oe | = 
fa=(Frr) HL (0); me=(Sx2) HP nla) (4) 
andjwe are assuming the time factor e“‘, and the primary contribution from the dipole is 
g ’ p a 


— "BS (2n+-1)f4(br)va( EV )Pa(cos 6), when Ey m<b, 


and (5) 


— FD 2n+ 1)va(bi) fu (EV) Paco 6), when Evi >by 


where 
ae at Ne ; 5A 
by =29(R+HA)y i/d; 2n(x)=Fn(x) + n(2). (5A) 
2 The normal slow lapse rate of u will be approximately allowed for by taking an “effective’’ value for the earth’s radius. This is areasonable 
procedure if the inversions are not too high [cf Norton, 1959; Bremmer, 1960; Wait, 1959]. 
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The constant 8 is proportional to the square root of the power radiated: (see reference 1, 
page 111), R and H are, respectively, the earth’s radius and antenna height. 
Let the total field in the space 2;_;<&<y,, i.e., the region where x=x; be 


1B Sy(ant1) AY fale) + BY ng (EG) JPa(008 0) «+ (6) 
and within the jth inversion layer (y;<&<z;) let it be 
2S (2n+1) (Cu? Q+DLe" @}Pa(cos 6) (7) 
where u” (£), v (€) are any pair of fundamental solutions of (2) within the jth layer (yj;< 


E<z;). 


The field above the complete set of inversion layers is of the form 
iB , 
Tp (on +1) E,&,(&)P»(cos 0): (8) 
1 


It will be useful to note here that, in view of the continuity of the rising and descending field 
contributions at the lower boundary of the first layer, we have the formulas 


Ay =F,+¥n(b1) 


BY=4, (9) 
for b,<é<y;, and also, that 
Att =i: | eal (10) 


in view of the radiation condition at infinity. 


3. Boundary Conditions 


These conditions, expressing the continuity of the field and of its normal derivative at the 
interfaces of the successive inversion layers, read as follows: 


At §=2,.1 
AP fu 24-aVq)+ BP nal 24-1VG) =CYP Ul (25-1) + DEPOL (2) 
Vj { AL Fn (2s-1V¥ Ej) + BY ng (Z5-1Ve,) } =COP-? (uy? (2y-1) }’ +D, {og (25-1) }’. «os rs 
where j=2,3,..., p. 
At é=y;-1 
CEP ue —? Yy-1) ADE POR? (Yj—1) = ANA? Sa (Ys-rVkj-1) + Be? tn (Ys-rVej-1), 
Ce-? fun? (yg) }? ADE? {o-? yy-1) } Ve {Ae-? Fa(Ys—1V 69-1) (12) 
+ Be-? nn(ys-1vikj-1) }. oo 
where j7=3,4,..., p+. 
If we let 
BY /AY=—H,... (13) 
we find, after somewhat involved algebra 
H.-G et (as j<pti) (14) 


Py HQ;-1 
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where 


Qj 1 = NOI 5, (yy reg) Sn Sprig ELI” $4 (yy ava) £4 2-1-¥e,) 
—ME-” 56 (yy_aVigg—a) Sa Z5—1Veg) —ME ~? 6 (yy—rVieg a) $4 (2 5-1-V 0g) 
= NIV (yy y—1) Fn (2) -1V ey) FLO nn (ya ky 1) 64 (25-1-V 3) 

— MEY 9p (Yyrv 5-1) Sn 25-1V 6) —ME~? ny (Yy—aV 1) Fa 2y-aV 4) (14a) 


the bar means that we interchange the ¢ and 7 functionality symbols in these formulas, and 


LO = Kj 411U? (25)0P (ys) — Us? (YO? (24)] 


MP = {uP (2) }’0P (y) —uF? (yy {vo (2,)}’ (14b) 
MY =a)" {ul (z)(0P (yd) {uP (ys)}0 (29)} 
NO= Eula}! {On (ys) }’— {Un (yp) }? {On (2,)} 7]. (14B) 


In cases where the earth behaves effectively as a perfect conductor the surface condition is, 
as already mentioned, II=0, by (4) and (5) this amounts to 


{ Fut tu(bs) péatey w+{ Gut Fa(b:) paulo) =0. (15) 


Hence, by some rather involved algebra, it is possible to derive an expression for the total 
Hertzian function at points below the lowest inversion similar in form to that developed in 


the writer’s 1952 paper, namely, 


in(b | 
m= — 18S (ont 1) COD (¢(c)na(&— malts) En(B)}Pa(eos@) (16) 


in which 
(16a) 


ss 
b=tVei3 m=z; c= 2eR/r 


and H, is the value of H; when j=1: it is expressible in terms of F’, and G, by (9). 
We have now arrived at the mathematical solution of the problem, for, since H,; can now 
be determined (see the argument immediately following), then (9) and (15) now amount to 


two simultaneous linear equations for F,, and G,. Determination of H,: 


We have, from (10) 


Ber 
Hy1=Aor 
=0 by (10) (17) 
we have 
H,=Q,/P, (17a) 


and then all the H’s from H,_; down to H, are given in succession by (14). 
4. Approximate Evaluation of the Field 


The field series (16) cannot, as it stands, be used to discuss the behavior of the field 
below the inversions because of its slowness of convergence. To make further progress, it 
is necessary to transform it into a more rapidly convergent series (series of residues) by complex 
-contour integration (Watson’s transformation), as was done in the writer’s 1952 paper on the 
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single high level inversion. As in that paper, the contribution of the imaginary axis to the 
complex contour integral involved is zero, and it appears that the series of residues giving 
the field is of the same general form as that obtained therein. It is as follows: 


_7if © §y, -1/2(0) my ,—1/2(@1) — $,-1/2(21) my, -1/2(0) 
sais r a" 21— {ny ,-1/2(21) }?[0Ai/08](5=» 


x { $y ,—1/2(21) my ,-1/2(€1) — mv ,-1/2(21) $y ,-1/2(E1) } sec v mP, -12(—Cos 6) (18) 





where the summation is over the roots v; of the equation 


$s—1/2(41)/ns—1/2(01) = Mi (19) 


considered as an equation in s. 


5. Thin Layers—Meter Wave Propagation 


It can be shown that, as in the case of a single inversion layer, the high level inversions 
have their maximum effect upon the propagation when they are sufficiently thin to act effec- 
tively as dielectric discontinuities. The condition that this may be so is that J 2p must be 
small compared with unity for each inversion, where J=2r7/k, p=h/R being the layer thick- 
ness and / the inversion height. We shall, therefore, confine ourselves to a discussion of this 
case only. 

When the layer thickness may be neglected, the preceding results simplify; we have 


L9-») ==(); NG-Y=—0 
Mi; MY? =—Vejaives (20) 


and since only the mutual ratios of Z9-», AZ,9-», M,9-” and N&» are significant, we may 
take 


Qa= VKj-1 t-172 (yj—1-Vj—1) Sy-1/2 (2y-1-V«y) —-vVk; $v-1/2 (y)-1-VK;-1) $12 (25-1 vj) 


Pji= Vkj-1 Ny-1/2 (y,—1 v K;-1) Sv-1/2 (25-1 V Kj) a VK; M172 (yy—1VKj—1) f-172 (25-1 V«;) 
= Vkj-1 Ny—1/2 (yj—1 Vkj-1) Nv—1/2 (2-1-V Ky) — Vk; M12 (yj-14 ky-1) My—1/2 (2-1 VK;) 


P3i= Vkj—1 (y,-1 v kj—1) Nv-1/2 (21 v kj) — V ky Sy-1/2 (yy-11 ky-1) Ny —1/2 (25-1 V k;) (21) 


It is possible to make further progress only when the heights of the high level inversions 
are large compared with R'*)?*. The inversions are then to be supposed as situated, above a 
well-known critical height which is of great importance in all diffraction theories for propa- 
gation along a sphere; asymptotic approximations may be applied there to all height gain 
functions. Fortunately, this condition is usually satisfied in practice, even for meter length 
VHF waves. It is then possible to show that the only roots v which could possibly give im- 
portant contributions to the field series (18) by reason of smallness of —Im (v/z’”*) fall into the 
following types: 

(A) those near * =z 


(B) those near é=y; (j=1, 2, ... , p) 


(A) corresponds to waves nearly grazing the horizontal level at the receiver and (B) to waves 
nearly grazing the lower boundaries of the successive inversion layers. 

It can further be shown that the contribution of roots of the group (B) to the field series 
is unimportant compared with that of group (A). Although limitations of space forbid the 


3 By “& near 2”? we mean |&—z| =O(zr'’). 
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inclusion here of the exact analysis, it may be remarked that the reason for this is to be found, 
not in any wide discrepancies in the magnitudes of their imaginary parts, but in the fact that 


- 3/2 
for v in group (B), {m,~1/2(21)}? (0H1/0s),s;-», becomes large of order Exp 5 2x (> } 
i.e., of order Exp (5 v2apy :) 
where 


pj=(y;—2)/z, by virtue of the factor {7,~1/2(21)}? {we have written e=(v—z)/x in the above}. 
This makes the leading terms of the series (18) which correspond to values of v in this 


group exponentially small. 
We shall therefore confine ourselves to consideration of values of v for which |y—z,| is 


O(x*).4 For such values, we may apply well known asymptotic properties of Bessel functions 
(cf Watson’s treatise or the writer’s previous paper, 1) and we obtain, from (21)* 


3 1/2 
FQ, ret | (2Dy1-+2)}™—( 4)" (2(Dy-a-0) 
J 
+0 ye]. : Exp —4 8x(p yt 
oi {2(pj1—e) } 3 a 
1 


re \L2 
Kp? P5151 {2(p4-1-) +2,}"2+( =) {2(py-1—€) +0j-1} 7 | + —e) yi 
Kj, {2(pj-1 €) | 








eo 1/2 
Kp? ast [ (20ry1-9) +2, } va_( 1) {2(pj-1—e) 
J 


—1 Se ' 
$24.2)" | - tap ayy Ex { 3 V20(Dy-1—€)?” i} 
1 


N12 
Kp? B,_y=—i [ 21-9 +2,) va4( 1) {2(py_1—€) + 05-1} nt ; (pase } (22) 
| _ J 








where 


kj;=1+0,, ete. (22A) 


Using these approximations in (14), we obtain, since v/r may be taken as effectively unity 





FL (2 mp1) +2j)"— (200-12) +2y-1}"4] Exp { —3 sBri(p, 9, } 
H _ L+H, [{2(pj-1—€) +05} "4+ {2(pj-1—€) + 95-1} 7] = (23) 
rs TL {2(py-1—€) +0; } "2+ {2(pj_1—€) + 0,1} /?]—1H, [{2(pj_-1—€) +05}? q 
— {2py-1— 9) +2y-1}"] Bxp { § v2zi(p,—9"" } j 








It transpires that Im (e) {which is always negative, so that p,-;—e is situated in the first quad- 
rant} is never very small compared with 2-2", but p>>2~?" and hence 


Exp 5 2a (pj-1—e)*/*|<1. 





Exp {-$ V2x1(p;-1—e)*”” D>; 


Again, it appears that H; is never large compared with unity. The denominator may then be 
effectively replaced by 
{2(pjy-1—€) +0} 72+ {2(py_1—e) +0,-1} 


4 See also note 2. 
5 Since zp3?>>1. See note 3. 
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and then (22) and (23) give 
Hj-1:=k;-1+H, (24) 


which yields at once 
A,=hith+ dai ie +k, (25) 


where 





, {2(p)— e042} F— {2(p;—e)+0,} 1”? . .. woes __ (3/2 
aa {2(py—€) O54. } 2+ {2(p; —e) +0; pie Exp v2 xi(D; e( 


=Q,/P;, by (22). (25a) 


But k; is the value which H, would have in the field series (18) if the jth inversion only were 
present. We obtain, therefore, the following important result: 

Addition Theorem. ‘When the inversion layers are at heights large, or moderately large 
compared with R'*)?%, and when also they are sufficiently thin to act effectively as dielectric 
discontinuities, then the value of H, to be used in the field series (18) and in the mode equation 


(19), is Dh; where k; is the value which H, would have in (18) if the jth inversion only were 
present.” 


6. Discussion of Results 


Even when the conditions of the above theorem are fulfilled, the resulting equation (19) 
for v is, unfortunately, still too complicated to allow of much further progress being made. 

There i is, however, a special case of the theorem in which the characteristic properties of 
the propagation can be clearly perceived. It occurs when the layers are relatively close 
together. As far as is known, this case is likely to occur in practice. It might, therefore, be 
useful to conclude by considering it. 


7. Case When Inversions Are Close Together 


In practice the dielectric constant differences ’ 0; are all small compared with p;. Since 
by hypothesis rp,?”>1 and «¢ is to be O(2~?*), (25A) gives the following approximate formula 
for k;: 


=— = 0441 — ¢)3/2 26 
k; ve es {- £ [axi(p, €) } (26) 

When the inversions are close together, we may take p, in (26) to be effectively independent 
of j; that is, we may take p; equal to an “effective” value pg. It would be appropriate to call 
Rpg the “cluster height” or “group height” of the system and this height would not differ much 
from the heights of any of the inversions in the case we now have in mind. Thus we have, 


in this case 
f - 3/2 97 
H,=— oer je p{-3 \2ri(pr—e) } (27) 


where A is the total drop in dielectric constant across all the inversion layers. 

Hence, in this case, the inversions behave like a single inversion layer across which the 
refractive index drop is equal to the algebraic sum of all the differences of refractive index 
across the separate inversions of the group. 


61.e., their widest separation is small compared with the height of the lowest. 
7 These are twice the corresponding refractive index differences. 
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8. Conclusion 


In the writer’s 1952 paper, it was found that a sharp high level inversion whose dielectric 
constant difference was about 2.58>107> could cause considerable anomalous propagation of 
meter radio waves. Inversions of this size do not, perhaps, often occur, but under settled weather 
conditions it is more likely that sharp high level inversions of much smaller dielectric constant 
difference would occur much more frequently. If several such small inversions occurred simul- 
taneously {e.g., four such that 0;—0,;,,—0.6 10-5 would do}, the above theory shows that 
the propagation would be effectively the same as that characteristic of a single large inver- 
sion. A mechanism has therefore been found which is capable of explaining the long distance 
fields without having to invoke the existence of a single strong elevated high level inversion. 


Note 1. The influence of scattering from atmospheric turbulence. The present paper has shown that 
when several high level inversions are simultaneously present their effect on propagation can be equivalent 
to that of a single well-marked high level inversion aloft, and Saxton [1951] in his experimental work on 9 and 
45 Me/s waves has found that the effect of scattering is, in general, likely to be unimportant compared with reflec- 
tion effects from (the equivalent of) a single inversion aloft, at any rate over the path considered by him. On 
the other hand, experiments elsewhere seem to indicate that scattering from turbulence can be at least as 
important a factor in the propagation [Northover, 1955}. 

The basic meteorological question is whether departures from the standard atmosphere structure, which 
are significant from the radio propagation point of view, occur in a stratified form (which may be viewed as 
anisotropic turbulence), or as isotropic turbulence. The present paper has attempted to make the theory 
sufficiently general to cover the most important possibilities of the first case. The best type of radio meas- 
urement for deciding between these alternatives appears to be space diversity measurement in mutually per- 
pendicular dimensions, The difference in the angular scattering characteristics of the two models could be 
ascertained, in theory, from such measurements. Experimental work along these lines has been carried out 
and results which have been obtained so far seem to favor the stratified model. 

Note 2. On the location of the roots v. When several high level inversions are present simultaneously 
we have found that there is one case only which is tractable mathematically, namely, that in which the inver- 
sions behave like a single high level inversion. This, therefore, is the only case which has been worked out 
in detail in the text and the disposition of v for it is known from the writer’s 1952 paper dealing with a single 
inversion layer. It was shown therein that the significant values of v were those “near’’ 2; 


fi.e., |»—2| =O(z"/8)}. 
Note 3. Asymptotic approximations to Hankel-Nielson functions. These are summarized on page 116 
of the writer’s 1952 paper. Using them we obtain, since 
xp¥!?>> 1; 3j-1<<pi-nj ve*/?=0 (1) 
—1/4 


; ir _ 1 2 -. 1 3/2 
- pie a bo—1/2(ys—1¥ “)~} Pi-1—€ ts Oj-1 exp —3 Vaix (Pie + 52-1 


1/4 3/2) 


2 <. 1 
V2ix sel f 


ae a 1 rh 
au! ron sila V2 exp{ —% 
1 —1/4 _2 1 3/2 


ix, — = gt se 2 j-. = Re 
2M Ee agave} i- “ ] i} v2 + nee 





in which 7; is the thickness of the jth inversion. The corresponding expressions in 7, 7’ are obtained from the 
above by changing the sign of the coefficients of 7 in the above. 
This enables us to obtain (22) from (21) neglecting the term 7;-;/R and the terms 0;_,, 0; in comparison 
with p;-; in the exponentials. 
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A Few Observations of the Perturbations in the Phase of 


the Low-Frequency Ground Wave’ 
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76 


J. M. Ross and J. E. Kirch 


(December 5, 1960; revised March 6, 1961) 


The effects of anomaly-producing terrain features on the phase of the low frequency 
ground wave were measured in a nearly idealized environment to confirm theoretical predic- 
tions: Results for an isolated butte gave good correlation with the theoretical model. In 
most cases the mean value of the phase perturbation approached zero. 


1. Introduction 


Low frequencies have been used for radio naviga- 
tion for some time with excellent results obtained 
over water and in the air. However, in use on the 
ground, disturbances caused by cultural and terrain 
features, and changes in conductivity cause errors 
in position. 

Many tests [Bateman, 1950, 1951; Schneider, 
1952; Pressey, 1956] of various systems have been 
conducted and considerable information is available 
in unpublished reports for these systems on the type 
of cultural and terrain features that cause propaga- 
tion disturbances. However, little information is 
found of sufficiently general nature to be applied or 
extended to other systems. In addition, propa- 
gational theory is usually developed only for idealized 
cases, so that measurements are necessary to deter- 
mine the limitations of idealized cases or mathe- 
matical models. 


2. Measurement Program 


A comprehensive program was conducted in the 
desert near Florence, Arizona to obtain a quanti- 
tative look at the order of magnitude and exact 
nature of these disturbances. The particular area 
was selected because it contained several anomaly 
producers in isolation. It contains mountains, 
buildings, power lines and telephones, all isolated 
from each other by relatively large areas of flat 
desert (see fig. 1). 

For example, the area of measurement was re- 
stricted so that the range to all transmitters was 
less than 80 km to minimize sky wave reflection and 
effects of earth curvature. Measurements were 
made with a double-heterodyne, CW, hyperbolic, 
phase-measuring equipment operating at 480 and 
450 ke/s. 

Further limiting of variables was achieved by 
taking readings along a single line of position in 

! Contribution from Motorola, Inc., Military Electronics Division, Western 
Center, Scottsdale, Arizona. 

?This paper was presented originally at the Seventh Annual East Coast 
Conference on Aeronautical and Navigation Electronics, October 24-26, 1960, 


Baltimore, Maryland, Portions of this work were performed for the Bureau 
of Ships under Contract N Obsr 72739. 
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lieu of using position fixes. Even with the single 
line of position, particular care was taken to position 
the receivers so as to measure the effect of phase 
changes from only one transmitter. 

The principle of operation is common to several 
different systems operating on different frequencies. 
It has been described in the literature and will not 
be covered here [International Hydrographic Office, 
1956 and Kirch, 1959]. It was found that equip- 
ment developed for tactical applications could be 
modified simply to serve as scientific measuring 
instruments. Three measuring stations were used; 
a fixed monitor to assure equipment stability, a 
ground van, and a helicopter. Measuring points 
were tied by first-order triangulation and trilatera- 
tion to the U.S. Coast and Geodetic survey of the 
area. All data were recorded by photographic 
techniques and reduced by an IBM 650 computer. 
Reported here is error in line of position data as 
determined by the phase of a 480-ke/s carrier. The 
errors are given in meters and degrees of phase with 
respect to distance expressed either in meters or 
wavelengths. 
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Figure 1. Map of test area. 





2.1. Obstacle-Free Area 


Measurements were made in a relatively obstacle- 


free area. A triangulation station near Randolph, 
Arizona was used as an obstacle-free area because 
the nearest mountain is 24 km and the entire area is 
surrounded by relatively flat desert thought to be 
homogeneous. 

Bremmer [1949] derived a rigorous solution suitable 
for numerical computation for the propagation 
of the low-frequency ground wave over a homoge- 
neous spherical earth. At short ranges, the series 
approximations involved converge very slowly. 
However, use can be made of the plane earth solu- 
tion due to Norton [1941]. Both solutions show 
the field to consist of the primary factor and the 
secondary factor. The primary factor, or undis- 
turbed field, propagates with a velocity determined 
by the index of refraction of the atmosphere at the 
earth’s surface. The secondary factor includes all 
disturbances to the field due to the presence of the 
earth and the antenna (near field effect). 

At the ranges used in this research, the only 
secondary effect is that due to retardation of the 
wave by the earth acting as a lossy boundary, since 
measurements were made at ranges large enough 
that the near field had no effect. Neglecting near 
field effect then, the phase correction due to ground 
absorption increases monotonically with distance. 

Since calculation of the secondary factor is quite 
complicated and because the conductivity of the 
earth is difficult to measure, only the primary factor 
was taken into account in the data reduction 
computation. 

Measurements were taken along a line extending 
914 m south, 305 m east, and 305 m north of the 
monument. Results appear in figure 2. No sig- 
nificant effect of terrain clearance was noted. The 
mean error of 2.1 deg is less than that expected had 
the secondary factor been taken into account. A 
small cyclic variation occurred with a period of a 
half wavelength. This variation is thought to be 
caused by a mountain at 24 km. A telephone line 
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Figure 2. 


Phase perturbations when anomalies are not 
apparent. 
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in the immediate vicinity caused no noticeable 
effect, since the error was identical for travel parallel 
to and perpendicular to the line. The small cyclic 
error can be eliminated by integrating over several 
wavelengths. 


2.2 Effect of a Mountain 


Of particular interest is the case of an isolated 
mountain in the transmission path. At high fre- 
quencies the effect of mountains obstructing the 
transmission path can be treated by methods of 
physical optics; but another approach must be used 
when the obstacles are no longer large compared to 
a wavelength. Such is the case considered by Wait 
and Murphy [1956 and 1958] for the effect of the 
amplitude and phase of a plane wave scattered by 
a semi-cylindrical boss. The boss of finite con- 
ductivity is placed on an infinite plane that is 
perfectly conducting. A source is assumed to be 
located many wavelengths from the boss; hence, 
the incident wave can be considered to be a plane 
wave. Attenuation of the incident wave is neglected. 

On the surface of the plane the incident wave is 


H‘“—=Hye™ (1) 


where H, is the amplitude of the magnetic field and 
k the wave number. The scattered field due to a 
boss of radius r and conductivity ¢ is 


Hee——H, >) K,,(kr,o)H® (kz). (2) 
m=0 


K,,.(kr, o) is a series of constants dependent on kr 
and o for the index m. H,,® (kx) is a Hankel func- 
tion of the second kind of order m. The ‘phase 
lag” of the resultant field is 


ae 


The amplitude and phase of the total field (ie., 
incident field plus the scattered field) at various 
points around the boss, are plotted in the report by 
Wait and Murphy [1958]. 

On the near side of the boss, the variation in the 
phase is essentially a damped sinusoid decreasing in 
a direction away from the boss. On the far side, 
the error decreases monotonically in an essentially 
exponential manner away from the boss. 

Unfortunately, mountains of the type and shape 
of the model do not often exist in nature, especially 
in isolation. The mountain chosen (Walker Butte), 
however, did answer the required description ad- 
equately enough to enable test of the theory. This 
butte rises 500 ft above the relatively flat surround- 
ing desert and is approximately hemispherical in 
shape. Its relation to the transmitters is shown in 
figure 3. Data was taken for a distance of 3,000 ft 
starting at a distance 2,665 ft from the peak. Sur- 
veying was done using accurately known points. 
Aerail readings were taken at terrain clearances 
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of 30, 150, and 500 ft with no significant difference in 
results obtained between ground and aerial readings. 

Figure 4 shows plots of phase error versus distance 
from the butte as actually measured, and for the 
mathematical model taken from Wait and Murphy 
[1958]. As the figure indicates, experimental results 
agree remarkably well with the theory. The 
calculated correlation index is 0.5. 

At very great distances from a long ridge, the 
asymptotic expression for H® (kx) can be used. 


ka @ HE (ka) =a) F- (o(). (4) 


This equation would indicate that the attenuation 
of the scattered field at great distances is propor- 
tional to the square root of distance. At these 
frequencies, the cyclic error noted at Randolph, 
Ariz., could be attributed to the mountain 24 km away. 


2.3. Effect of a Critically Placed Power Line 


The perturbation of phase due to vertical parasitic 
reradiators has been covered theoretically and 
experimentally [Bateman, 1950, 1951; Hufford, 
1950]. <A literature search failed to uncover any 
theoretical treatment of the effects of power lines 
on phase and amplitude of the low-frequency ground 
wave [Reynolds, 1953]. Considerable data was 
taken in the vicinity of power lines and telephone 
lines. In the general case the effects were local, 
with most telephone lines having negligible effect. 

It has been suggested that power lines might be 
considered a Beverage antenna.’ The theoretical 
effects of the line are, to a great extent, determined 
by terminations, capacity and other electrical 
properties which could not be determined or con- 
trolled for a line in service. 
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Figure 3. Sketch of Walker Butte. 





§ Private communication: R. F. Linfield, R. H. Doherty, and G. Hefley, 
Evaluation of propagation aspects of cytac system. 
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Relation of power line to test transmitters. 


FicureE 5. 


To measure this effect, a unique transmitter loca- 
tion was chosen. The green transmitter was de- 
liberately positioned so that the signal would travel 
for 614 km along a 115-kv line (see fig. 5). The line 
turns at an angle of 30°, making it possible to take 
measurements along the extension of the line. The 
common transmitter signal crosses the line at right 
angles. 
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Figure 6 shows how the relative error varies with 
distance away from the power line. The test grid 
system at this location was not tied into the 
U.S.C.G.S. survey, so that absolute errors cannot 
be reported precisely. From Land Office survey 
data and other readings taken in the vicinity, 
there is reason to believe that the magnitude of the 
absolute error is two or three times that of the rela- 
tive error. At a triangulation station, 414 km along 
the extension, the absolute error reduced to zero 
when the secondary factor calculated for earth 
conductivity of 6 millimhos/m was taken into 
account. 
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3. Conclusions 


Phase perturbations of the low-frequency ground 
wave occur in areas where the cause is not apparent, 
but the mean value of the perturbations approaches 
zero. Measurements tend to confirm the model 
used for the perturbations due to an isolated moun- 
tain. Tests made near a power line would indicate 
that more investigation of the Beverage antenna as 
a model are in order. 
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Smooth Earth Diffraction Calculations 
for Horizontal Polarization 
L. E. Vogler! 


(December 1, 1960) 


This paper presents a simplified method of determining the attenuation relative to 
free space in the so-called far diffraction region for horizontally polarized radio waves 


diffracted over a smooth spherical earth. 


A criterion is given which permits use of the 


method not only for far beyond line-of-sight paths but, in many practical situations, at 


line-of-sight or even slightly within. 


Although numerical procedures have been de- 
veloped by various authors for the calculation of 
field intensity arising from the diffraction of radio 
waves around a smooth spherical earth [Norton, 
1941; Bremmer, 1949; Norton, Rice, and Vogler, 
1955], the actual computation is often laborious 
and time consuming. Because of the number of 
independent parameters associated with the problem 
and the mathematical complexity of the solution, 
many intermediate calculations and graph readings 
are necessary for a given propagation path. How- 
ever, for certain particular cases of importance a 
simplified method, explained in this paper, may be 
used to good advantage; furthermore, an improved 
criterion is given for determining the ranges of the 
distance, heights, frequency, and effective earth’s 
radius factor, k, for which this method is applicable. 

In the so-called “far diffraction” region field in- 
tensity is determined by the first term of the Van 
der Pol-Bremmer residue series [Bremmer, 1949]. 
This term consists of four factors containing, 
essentially, the dependences on the distance, two 
antenna heights, and the electromagnetic ground 
constants and polarization of the transmitted wave. 
The residue series is strictly valid only for the case of 
a smooth “airless” earth of radius ka where a is 
the actual earth’s radius; however, recent investiga- 
tions [Millington, 1958; Norton, 1959; Wait, 1959; 
Bremmer, 1960] have shown that for certain models 
of the atmosphere, atmospheric effects may be 
satisfactorily included by introducing corrections to 
the antenna heights and defining the factor & as: 


kf +(25)_ E : (1) 


where dn/dh is the gradient of the refractive index, 
n, with respect to height, h, above the earth’s 
surface. Determinations of the proper height cor- 
rections for various model atmospheres are contained 
in the above references. 


' Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 








By defining the quantities: 
Lo,1,2=PoCG OJ filedo, 1,2) (2) 


where { is a function of the polarization and ground 
constants, fne denotes the radio frequency in Me/s, 
Co= (4/3k)'“, and do, d;, d, are indicated in the geom- 
etry of figure 1, the attenuation relative to free 
space A, (in decibels below free space) may be 
expressed as: 
A=G(%)—F(m)—F(a2)—Ci (3) 

The distance dependence is thus contained in the 
function G(r). The height dependences are con- 
tained in F(z,) and F(a.) through the relationship 
dy, 2 (mi)& y(3k/2) hy, 2 (ft) 

In the particular case of horizontal polarization, 
and the further restriction that: 


K=0.017774 Cofat!{ {e—1}? 
+ {1.8 10%(mhos/m)/fme}2]-/4<0.01 (4) 


where ¢ is the dielectric constant of the ground re- 
ferred to air as unity and ¢ is the ground conductivity 
in mhos per meter, the term C; in (3) is very nearly 
constant as is the parameter 8) appearing in the 
definition of « in (2): C,;=20.67 and 6)=1.607 
Thus where the above restrictions hold, the attenu- 
ation relative to free space is given by: 


A=G (a) —F (a1) — F(a) —20.67 (5) 
where: 
29,1,2=1.607 C3 fY/3dp,1,9(mi), (5a) 
Qo= (4/3k)"%, (5b) 
dy, (mi) ~[(3k/2) hi, 2(ft)}”?, (5c) 


and with the conditions: (a) horizontal polarization, 
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and (b) K<0.01; see (4). Curves of the functions 
G(x) and F(z,,.) are plotted versus z in figure 1 
Note that F(z) +1.356 is asymptotically equal to G(x) 


=0.0925542—10 log x for large values of z. It | 


might be mentioned that a brief study of the param- 
eter K will show that only for unusual combinations 
of k, fac, € and o will K ever exceed 0.01. 

Since (5) is based upon only the first term of the 
residue series, a fairly good indication of its range 
of validity may be gained from the ratio of the sec- 





Curves for the calculation of A, the attenuation relative to free space. 


ond term, 7:, to the first term, 7,, of the series 
Thus, if we require the error in A to be less than some 
given value, say 6, we have: 
20 log |1+ 72/7; |< 20 log{1+|72/7T;\]<6 

This is the procedure used to obtain the A(x, 2) curve 
shown in the lower right hand corner of figure 1. 
For A to be accurate to within 1 db (approximately) 
it may be shown that: 


Yo— 2X, A(x,) —22A (22) >200 mi. (6) 
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Notice that in certain cases (5) is valid (at least 
to within 1 db) even for just within line-of-sight 
paths as long as the condition (6) holds. For 
instance, with k=4/3, fme=100, d)>=65.38 mi, and 
h.=32.81 ft, then for an h,=1640.42 ft, which would 
correspond to a radio line-of-sight path, (6) becomes: 


487.67 —427.25 (0.65) —60.42 (0) =210>200 mi 
and from eq (5) 


A= G (487.67) — F (427.25) — F (60.42) —20.67 
=18.3—7.8+36.2—20.67=26.0 


The value of A as read from the C.C.I.R. Atlas 
[Atlas, 1955] (assuming horizontal polarization and 
for either “sea” or ‘“land’’) is approximately 27.3. 

It may be seen from (6) that the first term of the 
Van der Pol-Bremmer residue series is applicable not 
only far beyond the horizon but even in some cases 
at the horizon; (in this regard it is of interest to 
note that the lower limit of (6) appears to be ap- 
proximately the point at which the C.C.I.R. curves 
depart from linearity). Furthermore, the ease with 
which the attenuation may be calculated using (5) 
is apparent from the above example. Curves for 
the calculation of A by this method for vertical 
polarization are now being prepared for publication 
in the near future. Since A 1s defined as attenuation 
relative to free space, its relationship to E, the field 
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strength expressed in db above 1 uv/m for 1 kw 
E.R.P. from a half-wave dipole, is simply: 


E=102.79—20 log dy (mi) —A 


oD 
In terms of basic transmission loss, DZ: 


Ly, = 36.58 +20 log fmedo(mi) +A. 


(8) 
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On the Theory of Mixed-Path Ground-Wave Propa- 
gation on a Spherical Earth’ 


James R. Wait 
(March 9, 1961) 


The problem formulated concerns the mutual impedance between two vertical dipole 
antennas A and B located near the surface of a spherical smooth earth. The path between 
A and B is made-up of a number of homogeneous segments where the surface impedance is 
constant. Various formulas are developed, for two- and three-section paths, which are 
suitable for computation. Certain limiting cases are discussed and where possible a physical 
interpretation of the results is given. Comparisons with previous work are made. 


1. Introduction 


In the computation of the ground-wave field of a radio transmitter it is often assumed 
that the earth is a homogeneous and smooth conducting sphere. While this is a satisfactory 
approach for many applications, such as estimating the coverage of broadcast transmissions, 
there is often a need for more precise determinations. For example, the phase characteristics 
of the ground wave are significantly influenced by inhomogeneities in conductivity that occur 
at coast lines. This is an important consideration in the prediction of errors in radio naviga- 
tional systems. 

An excellent semiempirical approach to the problem of calculating fields over an inhomo- 
geneous earth has been given by Millington [1949]. Sometime earlier Feinberg [1946] had 
formulated a general method for treating mixed paths over a flat earth. Similar results were 
obtained independently by Clemmow [1953] and Bremmer [1954]. The author [Wait, 1956a] 
showed that the integral formula for a two-section path could be easily derived by an appli- 
cation of the compensation theorem [Monteath, 1951]. It was also indicated in this paper 
[Wait, 1956a] that the results could be readily generalized to a spherical earth. In a further 
paper [Wait and Householder, 1957] extensive numerical results were given for propagation 
over a two-section path on a spherical earth. Furutsu [1956] has also considered propagation 
over a spherical earth under mixed-path conditions. His general method is based on solving 
the dual integral equations for the problem by an iterative procedure. 

It is the purpose of this report to discuss the theory for mixed-path propagation over 
both two- and three-section paths on a smooth spherical earth. Some attention is given the 
various representations which may be used in practical field computations. The situation 
when one terminal is at great heights is considered in some detail since this solution describes 
the radiation pattern of an antenna when the ground plane is inhomogeneous. Where possible, 
a physical interpretation of the results is given. 


2. Formulation for a Two-Section Path 


The mutual impedance Z,, between two vertical antennas located at A and B over ¢ 
spherical earth of radius a is considered. The situation is illustrated in figure 1 where a 


yas) 


8 


4 
> 





FicurE 1. Two-section path on a spherical earth. 


1Contribution from the Central Radio Propagation Laboratory, National Bureau of Standards, Boulder, Colo. 
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vertical cross section of the earth is shown. The great circle distance between A and B 
(measured along the surface of the earth) isd. The earth medium to the left of the boundary 
line has a conductivity o and dielectric constant e. The corresponding constants for the medium 
to the right of the boundary are o, and «. A variable P on the great circle path is a distance 
afrom B. Then, for a>d,, the surface impedance is 


. ner ete ee Ps egw |’? 
Z| tpow/(o+ ew) |!” | 1. +; -| (1) 
o+ lew 
while for a<d, 
e ° ; ~ « LE gw a/2 
Z & | ipqw/(o1 + 2€w) |! ‘f1— us — | (2) 
‘ oi +2€1W 


where up is the permeability of the whole space which is assumed to be constant. 

The mutual impedance between the dipoles A and B for the two-section mixed path illus- 
trated in figure 1 is denoted Z’,,. It was shown previously [Wait, 1956a] that it could be 
expressed in the form 


I(Zi,—Zqp) = (Zi—Z) / { (Hy Hi, )dS (3) 


where Z,, is the mutual impedance if the surface of the earth were homogeneous with surface 
impedance Z everywhere. Hot is the tangential magnetic field of dipole A over the homo- 
geneous earth while Hr, is the tangential magnetic field of dipole B over the inhomogeneous 
earth. The current in the dipoles are both taken equal to Z for convenience. The surface 
of integration S extends over the region of the earth which is characterized by a surface im- 
pedance Z,. Equation (3) also follows directly from the work of Monteath [1951]. 


3. Reduction to a One-Dimensional Problem 


Equation (3) is essentially a two-dimensional integral equation for the fields. On a 
straightforward application of the principle of stationary phase, the surface integral can be 
reduced to a line integral from a=0 to d; along the great circle path. Before stating this result 
it is convenient to introduce certain attenuation functions, W and W,, as follows 


Zoya" e- WUd,Z) (4) 
and . 
Ziyi eee 6-8 Wd, ZZ) - 


where /, and 7, are the effective lengths of the dipoles A and B. It is to be understood that 
W and W’ are also functions of h, and hy, the height of dipole A and B above the ground. The 
functions W and W’ are normalized such that they would approach 1 if the earth were flat 
and perfectly conducting—and provided h,=h,=0. Consistently in what follows it is assumed 
that h, and hy<d«a. 


The resulting one-dimensional integral equation is given by 


W’ (d,Z,Z:)~W(d,Z) ford<0 








Z —2)/* W(d—a,Z)W"(a,Z1,Z) da for d,>0 


" [a(d—a) |'” 
where No= (uo/€9) #2 1207. (6) 


af 1/ 
W'(d,Z,2)=W(a,Z) (5 ‘( 


This result was obtained in an earlier paper [Wait and Householder, 1957] where certain 
numerical results were given for a two-section path in the frequency range 20 to 200 ke/s. 


402 








In the present work, it is desirable to introduce certain dimensionless parameters. These 
are defined by 
ig= (ka/2)'* (Z/no), 
igi: =(ka/2)'* (Zi/no), 
r=(ka/2)'* (d/a), 
= (ka/2)'* (a/a), 
=(ka/2)'* (d,/a). 
Thus W’ (2,9¢,q) =W(2,q@ for <0, and 


, . 2\ie 1 W( (c—7,q)W’ (7,91,9) 
WwW’ rq, \~W( y = ay j, ae <% 7 an dz 7 
(29,0) a d+(= ) (q Y . [2 (x—2 yy? (7) 
for 2, >0. It is immediately observed that W’(7,q:,q) in the integrand may be replaced by 
W(2,q) since ?<a, over the range of integration. Therefore 


W'(2,4.9)=Weo+(4) (ao of nW( eae (8) 


An alternate form of this equation is easily obtained by regarding the left-hand portion of the 
path as a modification to the homogeneous earth of surface impedance Z,. Thus 


> , fa\ie2 r-1, wW( 
W’(2,0,Q0=N (x,0) +( = ) (I-H |, ae 7 D4 dz. (9) 


4. Alternate Representations 


Equations (8) and (9) are explicit integral formulas to permit the computation of the 
attenuation function W’ in terms of the attenuation functions W appropriate for a homogeneous 
earth. The latter are well-known from the theory of van der Pol and Bremmer [1937] and 
Fock [1945]. From their work it follows that 


_(sz\" co e~ 8 wi (ts—Ya) Wi(ts—Yo) 
W (7.0) =( 4 ) poe oo _t.—¢ wits) wits) (10) 





where 


r= (ka/2)'’? (d/a), Ya= (2/ka)*!? kha, Yo= (2/ka)"? khy. 
The coefficients ¢, are solutions of the equation 
w(t) — qui (t)=0 (11) 


where w(t) is an Airy integral and the prime indicates a derivative with respect to¢. In terms 
of Hankel functions of order one third, w,(t)=exp (—2zi/3)(—t/3)'?A3 (2/3)(—t)°7].. The 
notation used here is closely akin to that used by Fock [1945]. 

In a similar fashion, 





ra\'/2 > emit? w(t? —yq) wilt? —yp) 
9 
W(x,0)= (77 >> =1t—q w,(t”) wi(t;”) — 


where ¢! are solutions of 
wi (t)—qwi(t)=0. (13) 


Using (10) and (12), the integration in (8) may be readily carried out to yield 
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- ; <t\"2 en Le = —1)_ 7] Wi(t,—Yo) W(t —Yp) 
Waa) -Wad-(F) @-9 EEE HHP) wey we) 





(14) 


The double summation converges quite rapidly if x, is of the order of unity or greater and 
provided that either y, and-y, are not large compared with unity. 

When 2; is small, the double series expansion given in (14) becomes very poorly convergent. 
An alternate expansion for the case h,=0 is obtained by using the following representation 
[Bremmer, 1953 and 1958; Wait, 1956]. 


Weg)= >) Anem™tgin(z)"” ae 


m=0, 1, 2,..6 


where 


— [" 1 4 1 
7 ae Ss aa ey Ay=iy7 (1475) mid 


_iyn A.= 7 ate 
A y VF (1455 A= —(+et De 344)’ ete. 


Using this result along with the residue series representation for W(z—7,q) enables (8) to be 
expressed in the form 
w,(t,— Ya) 


W’ (29,0) =W(2,q) —iz*(q—q) >> = A,,e'*/4 Gf In—1 ee (16) 
s=1 m=0 w(t) 
where 
"2 A A m-1 A 
In-= | t on te-De(4) 2 dh. (17) 
2 0 


The integral J,,-; may be reduced by the following recurrence relation 


1 Le IES Oe m 
In=—- " aatin 's(y )™/?2— >: Im_, id (18) 
2 8 i 2 


Successive application of this equation enables J(»_1)2 to be expressed in terms of J) and J_j,p. 
These are given by 





. iA e~ . 
I=ents} ets dy —=— (e'1's—1), (19) 
0 3 
and 
"2 pitts , im\V2 —— 
Iip= ae | —- di=(*) e~*s erf (./—7t,21) (20) 
0 (x)! ts 


where erf (Z) is the error integral of argument Z. 

The double series expansion given by (16) is highly convergent if z is somewhat greater 
than 1 and z, is somewhat less than unity. In fact, if {gjr;|<1, only the m=0 terms need be 
retained. Furthermore, if in addition z;<1, 


T_syom2e~ **8(2,)'”. 


And thus, 





W’ (2,9,0) =W(2,q) Bron aL 2a} | (21) 


The term in square brackets can be regarded as the correction to the attenuation function as a 
result of the portion of the path from 7=0 to x. In terms of the original parameters of the 
problem, the square bracket term becomes 
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[1-(2)" = (okdy)" ‘| (22) 


This same correction emerged from the corresponding theory for the flat-earth case [Wait, 
1956b]. 


5. Solution for One Antenna at Great Heights 


When one terminal becomes elevated to a large height such that y, or y,>1, it is desirable 
[following Fock, 1945] to replace the attenuation function W by an auxiliary function V. This 
function V actually characterizes the far zone radiation pattern of a vertical antenna on the 
curved surface [Wait and Conda, 1958]. First it is noted that if |y—?¢|>1 

2 
wy (t—y) ~e~**4(y—t)-™4 exp | - 3 (y—t)*? |. (23) 
Furthermore, if also y>+¢, 
wi (t—y) cen Fy Ae 15 gto, (24) 
This latter asymptotic relation enables W(z,q) to be written in the form 


2 ae — pl /2 


Wage" 3 Ve—V Ye) a (25) 
where 
—iXt, . a 
V(X,q)=—in'*3 ee. (26) 


(t,—@)wi(t,) wi(ts) 


This suggests that the radiation pattern function V’(X,q,q:) for the mixed path be defined by 
the relation 
gi/2 





W’ (2,q,)=075" V" (eV Yor) gira (27) 
The integral formula for V’ may then be obtained directly from (8), thus 
ViX-X,gWXq) ,4 
7 sins 
V'(X,q,q) =V(X,) + _ mf ed eS lS (28) 


1/2 A 
7?) 0 (X)'2 
where 
vy reel HOV date 
ae VYa= 3 a 


It should be noted that X can be positive or negative depending on whether d is greater or less 
than (2ah,)'”. 
An alternate form for V’ is obtained from (9), thus 


La)t (X, 
po Ry 


DIX (29) 





(wi) '/2 


ad & 3 1X ¢) 4! 1 W(X 
J (Xg.q)=V(X0) + Eo | ( 
—V¥a ( 


where X,;=2,—-yy,. Since y,>1, the lower limit of this integral is effectively —@. 


Series formulas for V’(X,q,q,) can be readily obtained in the manner discussed above 
for W’(2,¢,¢:). However, these are applicable only in the present situation if X>0 since the 
series given in (26) is divergent for Y<0. A more suitable approach is to use the contour 
integral representation for V(X,q). This is given by 
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ie ea | ee | , 
\ (X,) =n p wy(t)— qu (t) w(t) dt. (30) 


The contour is chosen so that it encloses (in a clockwise sense) the poles of the integrand at 
t=t,. It is readily verified that —2z7i times the sum of the residues leads back to (26). It is 
now convenient to write 
=>) a,(q)X" (31) 
n=0 


where 
—inr/2 


wi(t—Yyo) ; 
a,(Q)= a = COE TIO) w(t) it end 


Methods for evaluating the coefficients a,(q) have been discussed by Logan [1959]. 
The above representation for V and the series form for the W given by (15) enable the 
integration in (28) for V’ to be readily carried out. The result is 


V(X, q; m)—s AX, g=*— “q >» a an(Q)An(n)em"*q ae m= =i (33) 


where 
m-1 


"dy Ci ee 
P,a=| (X—z)"(x) ° dz. (34) 
: 0 


The integral P,,™=! can be reduced by using the following formula 





_ (X—2)"(21)? nv Qn 
— m - m P-uZer (35) 
qt! gt} 





Successive applications enable P,,™-! to be expressed in terms of P m where 
2 o> 
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P “ALA m/2,7 (2)? 7? ‘ 
a=] eat (36) 


The above double-series representation for V’ converges very rapidly when both |X| and z, are 
reasonably small compared to 1. 


6. Three-Section Path 


The next most obvious extension is to a three-part medium. The situation is illustrated 
in figure 2. The approach used amounts to a successive application of the theory for the two- 
part medium. The path, between A and B, consists of three segments whose surface imped- 
ances are Z, Z, and Z;. The length of the latter two segments are d, and d,. Using the 
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Figure 2. Three-section path on a spherical earth. 
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above-mentioned method, it really follows that the expression for the resultant attenuation 
function may be written 


da 
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where W(d,Z) is the attenuation function characteristic of propagation from A to B over a 
homogeneous earth of surface impedance Z. W(d—a,Z) and W(a,Z,) are attenuation functions 
for propagation over homogeneous surfaces of surface impedances Z and Z,, respectively. The 
integration variable a can be regarded as a great circle distance measured from B. The 
W’'(a,Z,,Z:) occurring in the second integral is the appropriate attenuation function for propa- 
gation over the two-part medium from B to a point a on the middle segment. 

For numerical work it is again convenient to introduce the dimensionless parameters used 
for the two-section theory. Thus 





_ {- W(a2—7,q)W(2,q) a 
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where 


iqa= (ka/2)'*(Z2/m), t2= (ka/2)'*(d;/a), 
and the other quantities have their usual meaning. From (8) and (9) we may deduce that 


vp ph as de »\1/2 ; t—z Wi t—x’,Z Wi wy 
W (x 91,92) =W (@.n)+(=) aa) { 1W(r rier 2) 


. / pe \ U2 Pm -_ ) 
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(39) 





Of course it would be possible to obtain various series formulas for W’ for a three-section 
path, but in general these would be cumbersome. It is considered to be preferable to work 
directly with (38) and carry out the integration numerically since presumably W’ (2 ,¢1,q2) 
would be known from the two-section theory. 

There is a rather simple limiting case of (38) which is rather interesting. We take q=q, 
and 2.<z and x. To this approximation, the first integral on the right-hand side of (38) 


_vanishes and the integrand in the second integral is essentially a constant over the range 2, 


to x7, +2. Thus 


H i) 


ae “— 
W'(x,q.q,0)=We.+(2) Cg —g) In(e—n))2 W(a—21,q)W(21,9). (41) 


This simple formula has a clear physical interpretation. The second term on the right-hand 
side can be regarded as the field scattered by the strip of width z.. It is proportional to q.—gq 
(or Z,—Z) which is the surface impedance contrast and it also is proportional to the attenua- 
tion functions W(r—2,,q) and W(2,,q) which respectively account for the attenuation from 
A to the strip, and from the strip to B. Such an equation is analogous to the first Born 
approximation in scattering theory. 

In the foregoing discussion of the three-section theory the antenna heights h, and hy, of 
A and B respectively, may be arbitrary (provided that ha and h,<a). When the height of 
one terminal is large, such that y,>1, for example, it is again convenient to define a radiation 
pattern function V’. In this case, 
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Ww’ (2,9, 42,91) = e718 V"(X,9,02,0) an (42) 
where 
” , ( 1 V(X—2,q)W(z A 
(X,9,92,91) = V(X,g) + a ah 10) gi 
—q) (+1 (X—2, gw (zr, 1,92) 7 P 
+! is J (x)¥2 da (42) 


where all quantities on the right-hand side have their usual meaning. 
7. Some Extensions of the Theory 


It is of interest to apply the compensation theorem as expressed by (3) to a spherical 
earth in the case when the surface impedances are a function of distance along the great circle 
distance d between A and B. The one-dimensional form of (3) may be written in terms of the 
attenuation functions as follows 


W(r—2 ig) We, 91,9) 7A a (44) 
[?(x—z)]} 





W'(gare)=Wa2)+(4)- [“act)—a) 


where g; and q are now functions of Z (or a) the great circle distance from B. Of course, if 
q, and q are constant, (44) is identical to (7). 

It is believed that (44) could be solved directly by numerical means if q¢ is taken as a con- 
stant corresponding to the unmodified homogeneous earth and q,(7) is some specified function 
of 2 over the range from 0 to z,. Equation (44) also permits an approximate treatment of 
transition regions. To illustrate this latter point we consider a three-section path as shown 
in figure 2. Now, however, the middle section is considered to have a surface impedance which 
is a function of a. Thus, q is a function of 2 over the region from 2, to 2;+22. The attenua- 
tion function for this particular path may be expressed by 


1/2 (°ry+22 , 
W' (2,9,92(2),)=W’ (2, oa) +(=) Bi [ao(t)—q) Wt. W" (2 g1,42(2).9) pa 


iz(2—z)}* 





(45) 


where W’(z,q,q:) is the appropriate two-section attenuation function and is given explicitly 
by (8). Now on the assumption that 2, is small compared with both 2, and 2, it is permissible 
to replace W’ in the integral in (45) by W(2,,q,)._ Admittedly, this is a first-order approxima- 
tion but it should be adequate since the relative contribution of the integral is expected to be 
small in any case. To the same approximation W(r—7,q) in the integrand can be replaced 
by W(x—2,,q). Thus 


. y yo 1/2 — Do y x 
W’(21,9,92(2),qi:) &V "(2.9.0 +(~) (q2—q) iat—a)=" (x—271,qg)W (21,41) (46) 


where 
1 I1+7q a. r A 
Q=— q2(z)dz is the average value of q(z). 


8. Appendix 


In the main text of the paper the attenuation function W for a homogeneous earth con- 
sistently appears. It is of interest to see that an integral equation for W in the case h,=h,=0 
emerges directly from (44). To show this we allow x, to approach zero in (44) and then we 
imagine the unmodified surface to be a flat, perfectly conducting plane between A and B. 
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Thus, the appropriate value W(z,q) is simply exp (iz3/12). (The quantity z°/12 is the difference 
between the arc and the chord between A and B.) Furthermore, the appropriate value of the 
surface impedance along the are AB over a conducting plane is easily found to be given by 


d—a 
Z=—n 2a 





or 
q(2) ~i(a—2)/2. 
Thus 


1/2 (°z —r A A A A 
W(a,a(a)) =e4(2) "| gui SS? | ete Prong. (47) 
0 ae 


Now if g, is a constant throughout the length of the path, W’ in the above equation becomes 
by definition the attenuation function for a homogeneous path. Thus 


. 1/2 (°2 =) a ee eT A 
Wx.a)=e"+(2) i} [ o— ‘ SS] [2(a—z)]~12et#@-) 92 Wz, q:)dz. (48) 
0 2 


This integral equation has been derived by Hufford [1952] who also shows that its solution 
is the residue series 
e7 itr 


—; (49) 
t;— 1 


e rr 1/2 
We-(F) > 
v r=1,2,3... 
where ¢t, are roots of 
w, (t)—qrw, (t) =0. 


In a similar manner it may be shown that 
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These particular integral equations have also been discussed by Logan [1959]. 
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Selected Abstracts 


Two magneto-ionic phenomena permitting the observation 
of valley minima between the E and F regions in the Arctic, 
J. W. Wright, pp. 85-93 of Some ionospheric results obtained 
during the International Geophysical Year; Proc. symp. 
organized by the URSI/AGI committee, Brussels, 1959, edited 
by W. J. G. Beynon (Elsevier Publ. Co., Amsterdam, 1960). 
Two phenomena are described which result from the influence 
of electron collisions on radio wave reflections from the E 
region at high magnetic latitudes. One permits observation 
of the depth of the ‘‘valley”’ between the F and F regions; 
the other permits the measurement of electron density at a 
fixed and known height in this valley, when it exists. 


The CRPL electron density profile programme: Some features 
and early results, J. W. Wright, pp. 215-220 of Some iono- 
spheric results obtained during the International Geophysical 
Year; Proc. symp. organized by the URSI/AGI committee, 
Brussels, 1959, edited by W. J. G. Beynon (Elsevier Publ. Co., 
Amsterdam, 1960). 

Electron density data and other quantities (hmax/2, scale 
height, total electron content) provided by reduction of 
vertical soundings to N(h) profiles at CRPL, are described. 
A model of the F region above hmax/'2 is developed, corre- 
sponding to a Chapman profile of 100 km scale height, and is 
shown to be in agreement with available rocket and moon- 
echo data. Vertical cross sections of the ionosphere across 
the geomagnetic equator are illustrated and discussed. 


Pecularities of the ionosphere in the Far East: Sporadic E 
and F region scatter, E. K. Smith, Jr., and J. W. Finney, 
pp. 182-191 of Some ionospheric results obtained during the 
International Geophysical Year; Proc. symp. organized by 
the URSI/AGI committee, Brussels, 1959, edited by W. J. G. 
Beynon (Elsevier Publ. Co., Amsterdam, 1960). 

Through a study of ionosonde data from the world-wide net- 
work of stations and of miscellaneous oblique-incidence 
field-strength measurements made in Japan and the United 
States it was discovered shortly before the beginning of the 
IGY that sporadic EF appears considerably more intense in 
the Far East than in similar latitudes in the Western Hemi- 
sphere. Comparisons of sporadic E data are very difficult 
unless made under identical conditions. Therefore an experi- 
ment was designed which consisted of recording transmission 
loss over two matched 50 Me/s oblique-incidence circuits of 
approximately 800 miles in length, one in the Far East and 
the other in the Caribbean. The principal advantage of an 
oblique-incidence circuit over a vertical-incidence one at the 
equivalent frequency (about 11 Me/s in this case) is that 
although the amount of sporadic E observed would be com- 
parable, the effect of D region absorption in decibels at 
50 Me/s would be less than at 11 Me/s by a factor of about 4.5. 


Studies of scattering phenomena in the equatorial ionosphere 
based upon VHF transmissions across the magnetic equator, 
K. L. Bowles and R. Cohen, pp. 192-195 of Some ionospheric 
results obtained during the International Geophysical Year; 
Proc. symp. organised by the URSI/AGI committee, Brussels, 
1959, edited by W. J. G. Beynon (Elsevier Publ. Co., Amster- 
dam, 1960). 

Results have been attained pertaining to the phenomena of 
equatorial sporadic E, equatorial spread Ff, and propagation 
by means of ionospheric radio wave scattering at VHF near 
the magnetic equator. These studies are the results of IGY 
experimentation in South America utilizing 50 Me/s trans- 
missions crossing the magnetic equator at Huancayo, Peru, 
and over other oblique paths. Additional ionospheric experi- 
ments were also conducted in Huancayo at vertical incidence. 





It has been established that a close relationship exists between 
the magnetic manifestations of the equatorial electrojet cur- 
rent above Huancayo, the occurrence on Huancayo ionograms 
of equatorial sporadic FE, and the intensity and fading rate of 
VHF signals propagated by the ionosphere over Huancayo. 
The zone in which EF region effects of the equatorial current- 
stream are apparent extends over some 10° of latitude. Also, 
the scatter propagation in the equatorial ionosphere is 
characterized by higher signal-levels and by more rapid 
fading-rates than are observed in the ionosphere at temperate 
latitudes, and signal-strengths remain high throughout the 
night. The feasibility of utilizing ionospheric scattering for 
communications in equatorial regions seems quite promising. 
Curiously, ionospheric scatter propagation at VHF similar to 
that associated with sporadic # formations in temperate 
latitudes is observed from time to time just to the North and 
South of the magnetic equator, but its occurrence appears to 
be excluded in the immediate vicinity of that equator. 
Further a remarkable sensitivity to polarization was estab- 
lished in HF radar echoes obtained at vertical incidence from 
the diffuse region (bounded by the slant sporadic EF and 
extending to the q-type sporadic EF), comprising the typical 
magnetic-equatorial sporadic E configuration appearing on 
Huancayo ionograms. These diffuse echoes were obtained 
with an antenna-array polarized magnetic North-South, but 
were virtually eliminated upon using an East-West polarized 
antenna-array for transmitting and/or receiving. 

Ionospheric propagation by scattering from the F region was 
sought over a transequatorial path employing a 2580 km 
transmitter-receiver separation. Propagation by F-scatter 
was present only about 10 percent of the time, and this was 
during evening hours, when enhancements in signal-strength 
above the levels due to EZ region propagation were found to 
occur over this circuit, generally before midnight. A sufficient 
and nearly necessary condition for the occurrence of such 
openings was found to be the presence of equatorial spread F 
indications on the Huancayo ionograms. By the variation of 
antenna radiation patterns and by oblique incidence pulse- 
delay measurements, the height of the propagation medium 
sustaining this F-scatter was established, in general, to be 
identifiable with the lowest height of the associated equatorial 
spread F seen on Huancayo ionograms. Thus, the spread-F 
time delays are not the result of an intervening screen below 
the F region that is invisible on the ionograms, as had been 
suggested in certain theoretical explanations. Further, the 
received pulses over the oblique path were but little wider 
than those transmitted. This lack of pulse-broadening and 
the narrowing of the equatorial spread F' configuration on the 
Huancayo ionogram upon using a high-gain antenna-array 
there at vertical incidence, indicate that the equatorial spread 
F consists of a thin sheet of irregularities. From the pulse 
measurements, the thickness of this sheet is estimated at the 
order of 50 km. The not always perfect comparisons in the 
lowest heights of Huancayo spread F with heights deduced 
from pulse-delay in the F-scatter is interpreted as evidence 
for a limiting geographical extension of a given scattering 
sheet in the F region. Polarization measurements similar to 
those described above prove indirectly that the equatorial 
spread F echoes observed at Huancayo must arrive from the 
magnetic East-West direction; i.e., that orthogonality to the 
magnetic field is imposed as a requirement for obtaining them. 
This is good evidence that the irregularities are elongated in 
the North-South direction, along the earth’s magnetic lines 
of force. Thus, the time delays associated with equatorial 
spread F at Huancayo must arise from echoes to the Kast and 
West of the zenith. From the success in obtaining scattering 
at 50 Me/s up to heights of 500 km or more, it may be inferred 
that irregularities of scale size 10 metres or smaller, measured 
normal to the magnetic field lines, can exist up to these heights 
at night in the equatorial F region. 
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Standard frequencies and time signals from NBS stations 
WWV and WWVH, by Radio Standards Laboratory, NBS 
Misc. Pub. 236 (Dec. 1, 1960). 

Detailed descriptions are given of six technical services 
broadcast by National Bureau of Standards radio stations 
WWV and WWVH. The services include 1, standard radio 
frequencies; 2, standard audio frequencies; 3, standard time 
intervals; 4, standard musical pitch; 5, time signals; and 6, 
radio propagation forecasts. Other domestic and foreign 
standard frequency and time signal broadcasts are tabulated. 


Near infrared atmospheric transmission to solar radiation, 
D. “ Gates, J. Opt. Soc. Am. 50, No. 12, 1299-1304 (Dec. 
1960). 

Near infrared solar spectrum observations taken on October 
15, 1954 with a double-pass NaCl prism spectrometer have 
been analyzed for transmission coefficients for the ‘‘selective”’ 
absorption factor and for the ‘“continuum’’ factor. The 
analysis was carried out for 59 wavelength positions between 
0.872 and 2.5374. The monochromatic data fit well the law 
lnT=c,(w) 2 where w is the amount of water vapor in the 
optical path. The coefficient c; is given as a continuous 
function of the wavelength. A coefficient of extinction for 
the “continuum” factor is also given. 


Limitations of radiosonde punch-card records for radio- 
meteorological studies, B. R. Bean and B. A. Cahoon, 
J. Geophys. Research 66, No. 1, 328-331 (Jan. 1961). 
Instrumental and data reduction restrictions affect the con- 
version of radiosonde data into radio refractive index profiles. 
It is found that over half of the elevated radio ducts reported 
by the observer will be missed due solely to the standard 
criteria established for rejection of data points in the process 
of preparing punch-card records. 

An extensive collection of radiosonde data containing all 
observed levels is now available from some seventy locations 
on and near the North American Continent. 


VLF phase perturbation associated with meteor shower 
ionization, C. J. Chilton, J. Geophys. Research 66, No. 2, 
379-383 (Feb. 1961). 

Anomalies associated with the Lyrid, 6-Aquarid, and Perseid 
meteor showers were observed on the daily phase variation 
of the 16 ke/s transmission from Rugby, England, to Boulder, 
Colorado. The effective reflection heights estimated from 
these anomalies are 80, 81, and 83 km respectively. 


Microwave spectroscopy—Atomic frequency standards, J. M. 
Richardson, Encyclopedia of Spectroscopy, edited by G. L. 
Clark (Reinhold Publ. Corp., New York, N.Y., 1960). 

The topic of atomic frequency standards is surveyed in a 
general way with emphasis on principles and current results. 
Atomic transitions satisfy very well certain requirements of 
a good standard; namely, constancy, precision, renewability, 
and convenience. Technique of observation, results, and 
limitations of several types of atomic frequency standards 
are described. These are the gas absorption, atomic beam, 
maser, and optical pumping methods. he atomic resonance 
may be used as a resonator or as a frequency controlling 
element in a control mechanism. Useful applications are 
improvement in the time scale, analysis of power spectra and 
line shapes, and tests of certain postulates of special and 
general relativity. 

A theoretical study of sporadic-E structure in the light of 
ne  eneansaaeaie K. Tao, NBS TN87 (PB161588) 
($1.25). 

The theoretical aspects of the mechanisms of sporadic-E 
reflections are described from both the standpoint of a thin 
layer and a scattering model. For the thin layer model, thin 
dielectric layers which have various distributions of electron 
density are considered. It is also pointed out that the scatter- 
ing theory for which an autocorrelation function of the 
fluctuation of electron density is given by modified Bessel 
functions of the fourth through the seventh order is an 
available model for sporadic-E scatter. Moreover blobs of 
ionization which have a horizontal scale of the order of 200 m 
and a vertical scale of about 50 m are considered for sporadic-E 
scatter. The frequency and distance dependences of the 
oblique VHF propagation by means of the sporadic-F layer 
are discussed by comparing the theoretical results with 
experimental evidence. 








The NBS meteor-burst propagation project—a progress report, 
C. E. Hornback, L. D. Breyfogle, and G. R. Sugar, NBS 
TN86 (PB161587) (1960) $1.25. 

This report briefly describes a meteor-burst propagation- 
study program at NBS-Boulder Laboratories and presents 
some of the preliminary analysis results. Observations have 
been made with scaled systems over three different paths 
(Long Branch-Table Mesa, Norman-Fargo, and Barrow- 
Kenai) at frequencies of 30, 50, and 74 Mc/s. The recorded 
data is processed by a combination of manual and automatic 
methods. The preliminary results show about a 10-db diurnal 
variation ia threshold for a constant duty-cycle. Thresholds 
for a constant duty-cycle were observed to have an approxi- 
mate frequency dependence relative to 30 Me/s of 15 d 
lower for 50 Me/s and 30 db lower for 74 Me/s. There was 
no statistically-significant difference observed in the 
occurrence of meteor-bursts from a Poisson distribution. 


Oblique incidence receiving antenna array for a relative 
ionospheric opacity meter, A. Wilson, NBS TN78 
(PB161579) (Nov. 1960) 50 cents. 

Experimental measurements incidental to the design of an 
antenna for a relative ionospheric opacity meter (RIO 
Meter) are described. 

The frequency of operation is 50 Mc. The antenna require- 
ments are that the main lobe of the antenna is directed at 
23° above the horizon, the half-power beam-width in the 
vertical plane does not exceed 10°, the minimum front-to-back 
ratio is at least 13 decibels, and the side-lobe levels are at 
least 10 decibels below the maximum response in both the 
E- and H-planes. Since the antenna is for use in Alaska, it 
is to be of simple design and physically able to withstand any 
anticipated wind and ice loads. 

The final antenna design is an array of three stacked horizontal 
dipoles with two optimally spaced reflectors behind each 
dipole to obtain the required directivity. The narrow main 
lobe of the antenna directed at an angle of 23° above the 
horizon is obtained by properly spaced and phased dipoles 
above the ground. The half-power beamwidth in the vertical 
plane is computed to be 744°. The front-to-back ratio over 
the rear 180° sector is not less than 20 decibels, and the 
half-power beamwidth in the E-plane is 74°. 

Two complete receiving antenna arrays were constructed, 
adjusted, and tested. These antennas were installed in 
Alaska where they are now in use by an auroral transmission 
loss project. 


A study of auroral coruscations, W. H. Campbell and M. H. 
Rees, J. Geophys, Research 66, No. 1, 41-55 (Jan. 1961). 
Short period variations in the N2+(0, 0) auroral emission 
band represent 5 per cent of the total light in the 3914A 
region. These quasi-periodic coruscations have a dominant 
period of 6 to 10 seconds. They attain a maximum amplitude 
in the predawn hours and are closely related to magnetic 
field micropulsations and ionospheric absorption of cosmic 
noise. Spectroscopic triangulation showed that the variations 
originate in the EF region of the ionosphere. The profile of 
the electron density associated with the aurora was found to 
have a maximum value of 1.1X10® electrous per cubic 
centimeter at 98 kilometers. 


Reports on CSAGI Disciplines, Part Ai, The Airglow, by 
F. E. Roach. Annals of the International Geophysical Year, 
X 134-137 (Pergamon Press, New York, 1960). 

Hourly zenith intensities of the night airglow in rayleighs for 
IGY and IGC. Results from 28 stations. 


Propagation of electromagnetic pulses in a homogeneous 
conducting earth, J. R. Wait. Appl.Sci. Research, Section B, 
8, 213-253 (1960). 

A general analysis for the electromagnetic response of con- 
ducting media due to pulse excitation is presented. The 
treatment is based on the Laplace transform theory. First, 
a survey of the field is made and the limitations and scope of 
the previous work are pointed out. The theory of propaga- 
tion of a plane wave pulse in a conducting and homogeneous 
medium of infinite extent is then reviewed. The form of these 
results enable one to evaluate the relative importance of the 
conductivity and the dielectric constant. It is indicated, for 
sufficiently large times in the transient response, that displace- 
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ment currents may be safely neglected for sea water and for 
most geological media. Under this assumption, the wave- 
form of the electric field in a conducting medium is illustrated 
for the case where the source is an electric dipole energized by 
a step-function current. Results are also presented for 
exponential and bell-shaped source functions. The pulse 
shape of the field components is profoundly modified as they 
propagate through the medium. It is suggested that this 
property may be utilized in measuring distances in the earth’s 
crust. The more difficult problem of propagation in non- 
infinite conducting media is also considered. To account for 
the presence of the interface in a conducting half space (i.e., 
homogeneous flat ground), a rather involved analytical ex- 
pression for the transient fields is required. Certain special 
cases, such as a horizontal electric dipole at the interface, are 
illustrated by numerical results. The transient excitation of 
a wire loop lying on the surface of a homogeneous ground is 
also considered. Finally, transient coupling between pairs of 
parallel insulated wires grounded at their end points is treated 
as an extension of the earlier results. 


Ionospheric absorption at times of auroral and magnetic 
pulsations, W. H. Campbell and H. Leinbach, /. Geophys. 
Research 66, No. 1, 25-35 (Jan. 1961). 

A study in March and April 1960, showed variations in the 
auroral zone ionospheric absorption of cosmic noise to be 
closely related to magnetic field micropulsations and short 
period coruscations of \ 3914. At times of polar-cap type 
absorption, magnetic micropulsation amplitudes were dimin- 
ished. Auroral ionization in the FE region, estimated from a 
particular luminosity-height profile, accounted for 50 per 
cent, at least, of the cosmic noise absorption. 


Tests for regression coefficients when errors are correlated, 
M. M. Siddiqui. Annals of Math. Stat. 31, No. 4, 929-938 
(Dec. 1960). 

In a previous paper the covariances of least-squares estimates 
of regression coefficients and the expected value of the 
estimate of residual variance were investigated when the 
errors are assumed to be correlated. In this paper we will 
investigate the distribution of the usual test statistics for 
regression coefficients under the same assumptions. Applica- 
tions of the theory to the cases of testing a single sample mean, 
the difference between the means of two samples, the coeffi- 
cients in a linear trend and in regression on trigonometric func- 
tions will be discussed in some detail under an assumed covari- 
ance matrix for errors. 


Climatic charts and data of the radio refractive index for the 
United States and the World, B. R. Bean, J. D. Horn, and 
A. M. Ozanich, Jr., NBS Mono. 22 (Nov. 25, 1960) $2.00. 

The radio refractive index of air, n=1+ NX 10-, is a function 
of atmospheric pressure, temperature, and humidity and 
varies in a systematic fashion with climate. 

Included in this Monograph is a compilation of refractive 
index data. Data listings made up of observations from 45 
U.S. surface weather stations for 2-hour intervals over an 
8-year period are given. Mean values, maxima, minima, and 
standard deviations of the refractive index have been cal- 
culated and tabulated for these observations. Additionally, 
mean vertical profiles of the refractive index have been pre- 
pared for 43 U.S. upper air sounding stations from long-term 
means of pressure, temperature, and humidity. 

Earlier studies of refractive index climate are assimilated and 
put into perspective. One such study is an extensive analysis 
and mapping of the refractive index climate of the United 
States. A worldwide radio refractive index climatology is 
developed based upon monthly mean observations of pressure, 
temperature, and humidity. 

An important finding of these climatological investigations is 
the strong correlation of N with height. A reduced-to-sea- 
level value of the index, termed No, is used to eliminate this 
systematic height dependence. The surface value of N, N,, 
may be estimated four to five times more accurately from 
charts of No than from similar-sized charts of N, itself. 

From climatic charts of No, N, may be estimated at any given 
location in the United States throughout the day during every 
season. In addition detailed annual and diurnal cycles, as 
well as 8-year cumulative probability distributions, are given 
for 12 representative U.S. stations. 





On a worldwide basis, charts of mean No are presented for 
both summer and winter season. 


Diffractive corrections to the geometrical optics of low- 
frequency propagation, J. R. Wait. Electromagnetic Wave 
Propagation (International Conference Sponsored by the 
Postal and Telecommunications Group of the Brussels Uni- 
versal Exhibition), edited by M. Desirant and J. L. Michiels, 
pp. 87-101 (Academic Press, New York, N.Y., 1960). 

The influence of caustics in low frequency radio propagation 
is studied. Correction factors to be applied to the geometrical 
optical representations are presented in the form of graphs. 
Particular attention is paid to the treatment of the caustic 
which occurs at the geometrical horizon of the first hop sky 
wave. It is shown that the apparent infinity in the con- 
vergence coefficient can be removed when the wave nature of 
the problem is considered. Antipodal phenomena are also 
discussed. 


Notre. In this paper the publisher deleted the captions 
for Figures 7a, 7b, 8a, and 8b. These should read: 


Fig. 7a—The factor |F| for sea water, e=80e,, c=5 mhos/ 
meter, a= 4/3 < 4360 km. 


Fig. 7>—The phase, of the factor F, for sea water. 


Fig. 8a—The factor |F| for average land, «= 15e,, c=5 x 10-3 
mhos/meter, a= 4/3 < 4360 km. 


Fig. 8b>—The phase A, of the factor F, for average land. 


Furthermore, in these four figures the abscissae should be 
labelled “‘Angle below horizon (@—@,)/2.”’ Also, in equation 
(25) the factor (F) is defined by 

iP 


ep EP 
(=o e*4 
where |F| and A are the quantities plotted in Figures 7 and 8. 


Propagation of electromagnetic waves along a thin plasma 
sheet, J. R. Wait, Can. J. Phys. 38, 1586-1594 (Dec. 1960). 
It is shown that a thin ionized sheet will support a trapped 
surface wave. The effect of a constant and uniform magnetic 
field is to modify the phase velocity and polarization of the 
surface wave. The essential features are illustrated by nu- 
merical results for selected values of the electron density, 
collision frequency, and gyro frequency. The effect of locat- 
ing the plasma sheet near and parallel to a conducting plane 
is also considered. In this situation other modes of a wave- 
guide type are possible in addition to the surface wave. 


Ionospheric mapping by numerical methods, W. B. Jones and 
R. M. Gallet, Telecomm. J., Journal UIT, No. 12, 260e—264e 
(Dec. 1960). 

One of the most difficult problems encountered in forecasting 
conditions for long-distance radio propagation is the mapping 
of ionospheric characteristics. A satisfactory method has 
been developed for mapping the continuous time variations on 
a world-wide basis of any ionospheric quantity entirely by 
numerical methods and electronic computing. 

As an illustration, the method is applied to mapping the 
monthly median values of the critical frequency of the F2 
layer (foF2) for December 1957. Input to the computer con- 
sists of the data as they are tabulated at 109 available sta- 
tions. An analysis is performed by progressively fitting series 
of orthogonal functions in three coordinates: time, latitude, 
and longitude. Two main difficulties have been overcome: 
(1) the data were irregularly positioned in the two space 
dimensions, and (2) they were affected by statistical noise. 
The end product of the analysis is a numerical map, or mathe- 
matical function, defined by a table of coefficients, small in 
comparison with the number of input data. 

Many applications can be made using the table of coefficients. 
For example, a variety of graphical representations and sta- 
tistical tests for accuracy may be produced, as well as fore- 
casting ionospheric conditions using a suitable solar index. 
The numerical methods are general enough to be applied to 
any geophysical quantity much as meteorological or geomag- 
netic characteristics varying continuously with time on a 
world-wide basis. 
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Wilson and J. Mandel, Tappi 43, No. 12, 998-1004 (Dec. 
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The papers listed below will appear in early issues of Radio Propagation: 


VHF and UHF signal characteristics observed on a long knife-edge diffraction path. A. P. 
Barsis and R. S. Kirby. 


Frequency dependence of D-region scattering at VHF. J. C. Blair, R. M. Davis, Jr., and 
R. C. Kirby. 


Reflection from a sharply bounded ionosphere for VLF propagation perpendicular to the 
magnetic meridian. Douglass D. Crombie. 


Mutual interference between surface and satellite communication systems. W. J. Hartman 
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A note concerning the excitation of ELF electromagnetic waves. James R. Wait. 
Computation of whistler ray paths. Irving Yabroff. 
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The reflection of electromagnetic waves from thin ionized gaseous layers. F.H. Northover. 


On the theory of wave propagation through a concentrically stratified troposphere with a 
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The Proceedings of the Thirteenth General Assembly of the International Scientific 
Radio Union (held in London, September 5-15, 1960) is now available. This 600- 
page document may be obtained from the National Academy of Sciences, National 
Research Council, Washington 25, D.C., U.S.A. The price, including postage, is 
$5.00. 
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